Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,019 G→T 5.7% intergenic (+35/‑30) lomR → / → stfR pseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,0190GT 5.7% 331.6 / 6.5 113intergenic (+35/‑30)lomR/stfRpseudogene, Rac prophage lom homolog;Phage or Prophage Related; interrupted by IS5 and N‑ter deletion/Rac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (55/51);  new base T (2/5);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 2.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 121 total aligned reads displayed.

GGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTC  >  NC_000913/1428871‑1429167
                                                                                                                                                    |                                                                                                                                                     
      cAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGGAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAgg                                                                                                                                                >  1:1008967/1‑150 (MQ=255)
       aGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGGGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGc                                                                                                                                               <  1:1837858/150‑1 (MQ=255)
          ttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCggg                                                                                                                                            >  2:963726/1‑150 (MQ=255)
          ttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCggg                                                                                                                                            <  2:1509311/150‑1 (MQ=255)
          ttAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCggg                                                                                                                                            >  2:957204/1‑150 (MQ=255)
           tAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCTGAACCGGCGGGc                                                                                                                                           <  1:1032166/150‑1 (MQ=255)
             aCCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCTGAACCGGCGGGC‑tg                                                                                                                                        <  1:1229548/150‑2 (MQ=255)
             aCCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑tt                                                                                                                                        >  1:1657042/1‑150 (MQ=255)
               ccgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑tttt                                                                                                                                      >  2:492809/1‑150 (MQ=255)
                 gaccgaATCCGGGGACATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑tttttt                                                                                                                                    <  2:1397722/150‑1 (MQ=255)
                 caccgaATCCGTGGCCATTGCTATTGCTCATGTAGCCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGGCGGTTATAAGATCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCTGAACCGGCGGGC‑tgtcat                                                                                                                                    <  2:164482/149‑6 (MQ=255)
                  accgaCTCCGGGGCCATTGATATTGCTTATGAAGGCCCCGGCAGGGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGGTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTtg                                                                                                                                   <  2:790980/150‑1 (MQ=255)
                     gaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTTATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTgg                                                                                                                                <  2:1444513/150‑1 (MQ=255)
                      aaaccGTGGCCATTGATATTGCTTCTGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGGTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGGTCAGGCGGGC‑TTTTTTGTggg                                                                                                                               <  2:912277/147‑1 (MQ=255)
                      aaTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTggg                                                                                                                               <  2:787218/150‑1 (MQ=255)
                             ggCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAata                                                                                                                        <  1:502867/150‑1 (MQ=255)
                                  ttGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGAGGGGTGAATATGGCa                                                                                                                   <  1:1357632/150‑1 (MQ=255)
                                   tGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAg                                                                                                                  <  1:1239879/150‑1 (MQ=255)
                                   tGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAg                                                                                                                  <  1:1331181/150‑1 (MQ=255)
                                    gATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGt                                                                                                                 >  2:280130/1‑150 (MQ=255)
                                        ttGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATTGCAGTAAAg                                                                                                             >  2:468265/1‑150 (MQ=255)
                                           cTTATGAAGGCCCCGGCAGTGGTGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGAtt                                                                                                          >  1:1759839/1‑150 (MQ=255)
                                           cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGAtt                                                                                                          >  2:513789/1‑150 (MQ=255)
                                           cTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGAtt                                                                                                          >  1:194664/1‑150 (MQ=255)
                                              aTGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCa                                                                                                       >  1:1852679/1‑150 (MQ=255)
                                                taaGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAgg                                                                                                     <  2:432168/149‑1 (MQ=255)
                                                  aGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGgtg                                                                                                   <  2:1016951/150‑1 (MQ=255)
                                                  aGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGgtg                                                                                                   >  1:1813394/1‑150 (MQ=255)
                                                     ccccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAgtacacc                                                                                                <  1:1492679/150‑7 (MQ=255)
                                                     ccccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTAc                                                                                                >  1:164826/1‑150 (MQ=255)
                                                     ccccGGCAGTGGCGACTGGCGCACTGACGGGTTCATCGTGGGTGTCGGTTATAAGATCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTAc                                                                                                <  1:1470262/150‑1 (MQ=255)
                                                      cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACt                                                                                               <  2:1364628/150‑1 (MQ=255)
                                                      cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACt                                                                                               >  1:968778/1‑150 (MQ=255)
                                                      cccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACt                                                                                               <  2:1246028/150‑1 (MQ=255)
                                                       ccGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTg                                                                                              <  2:676790/150‑1 (MQ=255)
                                                         ggCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGaa                                                                                            <  2:866713/150‑1 (MQ=255)
                                                           cAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAg                                                                                          <  1:1729716/150‑1 (MQ=255)
                                                                 cGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGca                                                                                    <  1:115820/150‑1 (MQ=255)
                                                                  gACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGAGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGcac                                                                                   <  1:1307244/150‑1 (MQ=255)
                                                                           aCTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAAcc                                                                          <  1:455667/150‑1 (MQ=255)
                                                                             ggACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCgg                                                                        <  1:890254/149‑1 (MQ=255)
                                                                              gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGt                                                                       <  2:1297346/150‑1 (MQ=255)
                                                                              gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGt                                                                       <  2:1609632/150‑1 (MQ=255)
                                                                                cGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTAc                                                                     >  2:954835/1‑150 (MQ=255)
                                                                                cGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGAATAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTa                                                                      <  2:812401/150‑1 (MQ=255)
                                                                                      cctCGTGGGTGGCGGTTATAAGTTCTGAGTTGCCAGGTACCACAGTGGTATGACCGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGATAACCGGTACAGAACt                                                               <  2:394202/148‑1 (MQ=255)
                                                                                        tCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGc                                                             <  1:1314153/150‑1 (MQ=255)
                                                                                           tGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGcaca                                                          >  1:1531740/1‑150 (MQ=255)
                                                                                              gtgtCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATc                                                       >  2:281809/1‑150 (MQ=255)
                                                                                               tgtCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATcc                                                      >  2:19206/1‑150 (MQ=255)
                                                                                               ggtCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGTGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATcc                                                      <  2:510894/149‑1 (MQ=255)
                                                                                               ggtCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATcc                                                      <  1:654932/149‑1 (MQ=255)
                                                                                                  cGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGc                                                   <  2:1852679/150‑1 (MQ=255)
                                                                                                    gTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTg                                                 <  2:1366693/150‑1 (MQ=255)
                                                                                                    gTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTg                                                 >  1:187769/1‑150 (MQ=255)
                                                                                                    gTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTg                                                 >  1:1618621/1‑150 (MQ=255)
                                                                                                    gTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTg                                                 <  2:1792103/150‑1 (MQ=255)
                                                                                                         aaGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGc                                            >  1:874338/1‑150 (MQ=255)
                                                                                                          aGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCa                                           >  1:611791/1‑150 (MQ=255)
                                                                                                           gTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCaa                                          >  1:821857/1‑150 (MQ=255)
                                                                                                           gTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCaa                                          >  1:1503158/1‑150 (MQ=255)
                                                                                                             tCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCTCAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCaaaa                                        >  2:739101/1‑150 (MQ=255)
                                                                                                                gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGt                                     >  1:1727181/1‑150 (MQ=255)
                                                                                                                 aTTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAACCCGGTACAGACCTGCACAATCCAGCTGAAAGCAAAACGTa                                    >  1:585848/1‑150 (MQ=255)
                                                                                                                 aTTAGCCAGGTAACACAGTGATATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTa                                    <  1:877768/150‑1 (MQ=255)
                                                                                                                   tAGCCTGTTAACACAGTGTTATGACAGCCCGCAGGTTCAGGCGGGC‑TTTTTTGTGGGGTTAATATGTCAGTAAAGATTTCAGTTGTACTGAAAGAATGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAATCAAAACGTAAc                                  <  1:301639/150‑1 (MQ=255)
                                                                                                                   tAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAAc                                  <  2:1657042/150‑1 (MQ=255)
                                                                                                                   tAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAAc                                  <  1:1064725/150‑1 (MQ=255)
                                                                                                                    aGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACa                                 <  1:1581192/150‑1 (MQ=255)
                                                                                                                    aGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACa                                 >  1:74738/1‑150 (MQ=255)
                                                                                                                     gCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAg                                >  1:1525528/1‑150 (MQ=255)
                                                                                                                      ccAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGc                               <  2:1503158/150‑1 (MQ=255)
                                                                                                                         ggTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGcacc                            >  2:366148/1‑150 (MQ=255)
                                                                                                                            aaCACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACg                         >  2:975692/1‑150 (MQ=255)
                                                                                                                            aaCACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACg                         >  2:118614/1‑150 (MQ=255)
                                                                                                                             acacAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGATTATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACgg                        >  1:315535/1‑150 (MQ=255)
                                                                                                                             acacAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACgg                        <  1:1726650/150‑1 (MQ=255)
                                                                                                                              cacaGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggt                       >  2:167360/1‑150 (MQ=255)
                                                                                                                               acaGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtg                      >  1:67428/1‑150 (MQ=255)
                                                                                                                                caGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGACCTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtgg                     >  1:1158879/1‑150 (MQ=255)
                                                                                                                                caGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtgg                     <  2:1024647/150‑1 (MQ=255)
                                                                                                                                 aGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggt                    >  2:1532637/1‑150 (MQ=255)
                                                                                                                                 aGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggt                    >  1:1106542/1‑150 (MQ=255)
                                                                                                                                  gtgtTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtg                   >  1:881633/1‑150 (MQ=255)
                                                                                                                                    gtTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt                 >  1:1409242/1‑150 (MQ=255)
                                                                                                                                     ttATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGgtggtggtg                <  2:1379877/150‑1 (MQ=255)
                                                                                                                                     ttATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGgtggtggtg                <  1:534127/150‑1 (MQ=255)
                                                                                                                                      tATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGa               >  1:1025356/1‑150 (MQ=255)
                                                                                                                                      tATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTCAATATGGCAGAAACGATTTCAGGTGTACTGAAAGACGGCACAGGAAACCCGGTACAGAACTCCACAATCCCGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGa               <  2:1325467/150‑1 (MQ=255)
                                                                                                                                      caaaaaaGCCCGCCTGCACCGGCGGGC‑TATTTTGTGGGGTGAATATGGCAGTAAAGCTTTCAGGTGTACTGAAAGACGGCCCAGGACAACCGGGACAGAACTGTACAAGCCAGCTGAAAGCACAACGTCGCAGCACCACGgtgggggtgg               >  1:86674/7‑149 (MQ=255)
                                                                                                                                        tGACAGCCCGCCTGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAac             <  1:393763/150‑1 (MQ=255)
                                                                                                                                          aaaGCCCGCCTGAACCGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAacac           >  1:1761169/3‑150 (MQ=255)
                                                                                                                                          aCAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAACCCGGTACAGAACTGCACAATCCCGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAaccc           >  1:1250134/1‑148 (MQ=255)
                                                                                                                                             gCCCGCCGGTTCAGGCGGGC‑TTTTTTGTTGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTTAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAATCGTAACAGCACCACGGTGGTGGTGAACACGCt        >  2:1128078/1‑150 (MQ=255)
                                                                                                                                              cccGCCGGTTCAGGCGGGC‑TTATTAGTGGGGTGAATATGGCAGTAAAGAAATCAGGAGTACTGAAAGACTGCCCAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACCTAACAGCACCACGGTGGTGGTTCACACGcag       >  1:938852/1‑148 (MQ=255)
                                                                                                                                                cgccgGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGc     >  1:529484/1‑150 (MQ=255)
                                                                                                                                                 gccgGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGcc    <  1:724862/150‑1 (MQ=255)
                                                                                                                                                 gccgGTTCAGGCGGGC‑TTTTTTGTGGGGGGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTAGTGAACACGCTGGcc    <  2:1230272/150‑1 (MQ=255)
                                                                                                                                                  ccgGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCTCGGTGGTGGTGAACACGCTGGCCt   >  2:228221/1‑150 (MQ=255)
                                                                                                                                                   cgGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTc  <  2:621761/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGC‑TTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTC  >  NC_000913/1428871‑1429167

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑