Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,678,604 G→T A60S (GCA→TCA)  ydgH → DUF1471 domain‑containing protein YdgH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,678,6040GT100.0% 401.8 / NA 114A60S (GCA→TCA) ydgHDUF1471 domain‑containing protein YdgH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (54/60);  total (54/60)

Only 100 of 116 total aligned reads displayed.

CTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAA  >  NC_000913/1678467‑1678740
                                                                                                                                         |                                                                                                                                        
cTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCa                                                                                                                                        <  1:1227557/140‑1 (MQ=255)
  gCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGGAGCGGTGAAAGCCGTTTCTCGTCGCTCAGa                                                                                                                                      >  2:55672/1‑140 (MQ=255)
      tGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATaaa                                                                                                                                  >  1:620209/1‑140 (MQ=255)
      tGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATaaa                                                                                                                                  >  1:1667113/1‑140 (MQ=255)
        gCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAaaga                                                                                                                                <  2:947979/140‑1 (MQ=255)
           ttCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAgg                                                                                                                             >  1:1147162/1‑140 (MQ=255)
           ttCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAgg                                                                                                                             >  2:358569/1‑140 (MQ=255)
           ttCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAgg                                                                                                                             >  2:20054/1‑140 (MQ=255)
           ttCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAgg                                                                                                                             >  1:426834/1‑140 (MQ=255)
              tcCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgc                                                                                                                          <  2:561632/140‑1 (MQ=255)
                 gTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgc                                                                                                                       <  2:559586/140‑1 (MQ=255)
                 gTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgc                                                                                                                       <  1:1172100/140‑1 (MQ=255)
                 gTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgc                                                                                                                       <  2:1036869/140‑1 (MQ=255)
                 gTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgc                                                                                                                       <  1:1367462/140‑1 (MQ=255)
                 gTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgc                                                                                                                       <  1:1684466/140‑1 (MQ=255)
                  ttAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGTGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgcc                                                                                                                      <  1:920918/140‑1 (MQ=255)
                  ttAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgcc                                                                                                                      <  1:1189262/140‑1 (MQ=255)
                  ttAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgcc                                                                                                                      <  1:724673/140‑1 (MQ=255)
                  ttAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgcc                                                                                                                      <  1:432878/140‑1 (MQ=255)
                   tAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCct                                                                                                                     >  1:267006/1‑140 (MQ=255)
                        cagcaACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCtttt                                                                                                                >  1:531627/1‑140 (MQ=255)
                        cagcaACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCtttt                                                                                                                >  2:1738362/1‑140 (MQ=255)
                        cagcaACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCtttt                                                                                                                >  1:1331884/1‑140 (MQ=255)
                          gcaACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTa                                                                                                              >  1:163648/1‑140 (MQ=255)
                              cAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTAtgtt                                                                                                          <  1:1529595/140‑1 (MQ=255)
                                   cTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGa                                                                                                     <  1:569393/140‑1 (MQ=255)
                                   cTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGa                                                                                                     <  1:20054/140‑1 (MQ=255)
                                        acCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACtt                                                                                                >  2:1317720/1‑140 (MQ=255)
                                        acCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACtt                                                                                                >  2:1001806/1‑140 (MQ=255)
                                        acCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACtt                                                                                                >  1:375397/1‑140 (MQ=255)
                                        acCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACtt                                                                                                <  2:669900/140‑1 (MQ=255)
                                           ggAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTg                                                                                             <  1:958083/140‑1 (MQ=255)
                                           ggAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTg                                                                                             <  2:1331884/140‑1 (MQ=255)
                                            gAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGa                                                                                            >  1:818397/1‑140 (MQ=255)
                                             aGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGGAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGAt                                                                                           >  2:1789226/1‑140 (MQ=255)
                                             aGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGAt                                                                                           <  1:1246787/140‑1 (MQ=255)
                                              gCAAGCGGCATCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGAtt                                                                                          >  1:1507943/1‑140 (MQ=255)
                                                     gCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTaa                                                                                   <  1:1462400/140‑1 (MQ=255)
                                                     gCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTaa                                                                                   <  2:1147162/140‑1 (MQ=255)
                                                     gCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTaa                                                                                   <  1:1001806/140‑1 (MQ=255)
                                                       aGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACa                                                                                 >  2:159931/1‑140 (MQ=255)
                                                       aGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACa                                                                                 >  2:1677137/1‑140 (MQ=255)
                                                       aGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACa                                                                                 <  1:1738362/140‑1 (MQ=255)
                                                         cGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGc                                                                               <  2:880198/140‑1 (MQ=255)
                                                            ttAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGt                                                                            <  2:1757823/140‑1 (MQ=255)
                                                                            tgtAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCg                                                            >  1:567818/1‑140 (MQ=255)
                                                                            tgtAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCg                                                            >  1:483122/1‑140 (MQ=255)
                                                                            tgtAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCg                                                            <  2:267006/140‑1 (MQ=255)
                                                                              tAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  2:1447967/140‑1 (MQ=255)
                                                                              tAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  2:163648/140‑1 (MQ=255)
                                                                              tAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  1:1266964/140‑1 (MQ=255)
                                                                              tAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  1:876814/140‑1 (MQ=255)
                                                                              tAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  2:655627/140‑1 (MQ=255)
                                                                                 tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  1:1458661/140‑1 (MQ=255)
                                                                                 tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  1:375704/140‑1 (MQ=255)
                                                                                 tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  2:1031366/140‑1 (MQ=255)
                                                                                 tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  2:483122/140‑1 (MQ=255)
                                                                                 tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  1:34078/140‑1 (MQ=255)
                                                                                 tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  1:195744/140‑1 (MQ=255)
                                                                                   gTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACct                                                     >  1:1145500/1‑140 (MQ=255)
                                                                                     tACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACctct                                                   <  2:378970/140‑1 (MQ=255)
                                                                                      aCCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGTGTGTGGTCGCTGACCTCta                                                  >  2:847157/1‑140 (MQ=255)
                                                                                      aCCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCta                                                  >  2:514194/1‑140 (MQ=255)
                                                                                      aCCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCta                                                  >  1:1023347/1‑140 (MQ=255)
                                                                                        cGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCtata                                                >  1:1053609/1‑140 (MQ=255)
                                                                                        cGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCtata                                                >  1:38932/1‑140 (MQ=255)
                                                                                          gTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATaaa                                              >  2:1459039/1‑140 (MQ=255)
                                                                                          gTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATaaa                                              >  2:912655/1‑140 (MQ=255)
                                                                                          gTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATaaa                                              >  1:955381/1‑140 (MQ=255)
                                                                                           tCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAg                                             >  2:1584363/1‑140 (MQ=255)
                                                                                               tttAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGa                                         <  2:1347109/140‑1 (MQ=255)
                                                                                               tttAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGa                                         <  2:375397/140‑1 (MQ=255)
                                                                                                     gCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGa                                   <  1:358569/140‑1 (MQ=255)
                                                                                                        aTTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaaa                                >  1:625657/1‑140 (MQ=255)
                                                                                                         ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaaaa                               >  1:559204/1‑140 (MQ=255)
                                                                                                         ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaaaa                               >  2:800332/1‑140 (MQ=255)
                                                                                                         ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaaaa                               >  2:720525/1‑140 (MQ=255)
                                                                                                         ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaaaa                               >  2:166359/1‑140 (MQ=255)
                                                                                                         ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaa                                 <  1:826472/138‑1 (MQ=255)
                                                                                                         ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaa                                 >  2:826472/1‑138 (MQ=255)
                                                                                                            gCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGTGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCa                            <  1:847157/140‑1 (MQ=255)
                                                                                                             cGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCag                           <  2:1147829/140‑1 (MQ=255)
                                                                                                                aGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCagaaa                        >  1:1701344/1‑139 (MQ=255)
                                                                                                                 gCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCagaaga                       >  2:89684/1‑140 (MQ=255)
                                                                                                                     tGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACa                   >  1:421490/1‑140 (MQ=255)
                                                                                                                        aaGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGt                >  1:459797/1‑140 (MQ=255)
                                                                                                                         aGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTa               >  2:112961/1‑140 (MQ=255)
                                                                                                                          gCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTaa              <  1:1317720/140‑1 (MQ=255)
                                                                                                                            cGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATc            >  2:228237/1‑140 (MQ=255)
                                                                                                                            cGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcg           <  2:567818/140‑1 (MQ=255)
                                                                                                                             gTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcg           <  2:625657/140‑1 (MQ=255)
                                                                                                                               ttCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTCTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATcgcg         <  1:1613743/140‑1 (MQ=255)
                                                                                                                                tctcGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGt        <  2:1053609/140‑1 (MQ=255)
                                                                                                                                 ctcGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTa       <  1:1459039/140‑1 (MQ=255)
                                                                                                                                   cgtcgCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAAt     >  1:12461/1‑140 (MQ=255)
                                                                                                                                    gtcgCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAAtt    <  1:800332/140‑1 (MQ=255)
                                                                                                                                     tcgCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTa   >  2:74220/1‑140 (MQ=255)
                                                                                                                                     tcgCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTa   >  2:1648401/1‑140 (MQ=255)
                                                                                                                                      cgctCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTaa  <  2:1725032/140‑1 (MQ=255)
                                                                                                                                      cgctCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTaa  <  2:38932/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                        
CTGCTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAA  >  NC_000913/1678467‑1678740

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: