Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 252,594 G→C 5.2% W98C (TGG→TGC yafO → mRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913252,5940GC 5.2% 376.5 / 9.9 115W98C (TGG→TGCyafOmRNA interferase toxin of the YafO‑YafN toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (52/57);  new base C (4/2);  total (56/59)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

Only 100 of 117 total aligned reads displayed.

CCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTG  >  NC_000913/252446‑252742
                                                                                                                                                    |                                                                                                                                                    
ccTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATTATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                     >  2:132902/1‑150 (MQ=255)
ccTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                     >  1:1008623/1‑150 (MQ=255)
ccTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                     >  2:659978/1‑150 (MQ=255)
ccTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGt                                                                                                                                                     <  1:1817349/150‑1 (MQ=255)
  tGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCTCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTg                                                                                                                                                   >  1:253665/1‑150 (MQ=255)
     ccATTATTCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTc                                                                                                                                                >  1:952704/1‑150 (MQ=255)
      cATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGCCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                               >  1:1099113/1‑150 (MQ=255)
      cATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCa                                                                                                                                               <  2:483555/150‑1 (MQ=255)
       aTTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTTGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAt                                                                                                                                              <  2:1619563/150‑1 (MQ=255)
        ttAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAACCATGCTTGCTCAtc                                                                                                                                             <  2:1819969/150‑2 (MQ=255)
          aaTCAAATTTGAGCGAGCTCCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATACACCAGAACGAATGACGAAGCGCATTTGGTATATTGACAGGGGGCGTGTGATGACCAAGCAGGGTTGCTCATTGc                                                                                                                                           <  2:291558/150‑1 (MQ=255)
            tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCa                                                                                                                                         <  2:703888/150‑1 (MQ=255)
            tCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCa                                                                                                                                         <  1:134097/150‑1 (MQ=255)
                aTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATTTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                                     >  1:639352/1‑150 (MQ=255)
                aTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                                     >  2:1059771/1‑150 (MQ=255)
                aTTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCt                                                                                                                                     >  1:1455386/1‑150 (MQ=255)
                          gTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTAAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGCGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGaa                                                                                                                           >  1:803523/1‑150 (MQ=255)
                            tGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAAcc                                                                                                                         <  1:1804274/150‑1 (MQ=255)
                             gCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCt                                                                                                                        >  1:341848/1‑150 (MQ=255)
                              cTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTc                                                                                                                       >  2:440572/1‑150 (MQ=255)
                                cATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCAt                                                                                                                     <  1:1732355/150‑1 (MQ=255)
                                 atatTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                                                                                    >  1:1351471/1‑150 (MQ=255)
                                 atatTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                                                                                    <  1:373767/150‑1 (MQ=255)
                                 atatTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATa                                                                                                                    >  2:311787/1‑150 (MQ=255)
                                   atTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaa                                                                                                                  >  2:1452411/1‑150 (MQ=255)
                                    ttCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATaaaa                                                                                                                 >  2:1752253/1‑149 (MQ=255)
                                     tCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACt                                                                                                                <  1:249849/150‑1 (MQ=255)
                                       aTCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTgg                                                                                                              <  1:1208894/150‑1 (MQ=255)
                                        tCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGc                                                                                                             <  2:163494/150‑1 (MQ=255)
                                        tCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGCGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAATCTGGc                                                                                                             >  2:309367/1‑150 (MQ=255)
                                           ggCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAACCCTGAACCTCATAAACTGGCTCg                                                                                                          <  2:1520064/150‑1 (MQ=255)
                                           ggCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                                                                          >  1:919880/1‑150 (MQ=255)
                                           ggCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCg                                                                                                          <  1:583974/150‑1 (MQ=255)
                                             cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGACGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                                                                        <  1:365792/150‑1 (MQ=255)
                                             cAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCgag                                                                                                        >  1:207036/1‑150 (MQ=255)
                                                     aataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATaacaac                                                                                                >  2:839039/1‑150 (MQ=255)
                                                      atattCCATTTCCGCCAGAGTTTCGCAAATTCAGAAGTACGAATGACGAAGCGCATTTGGTTTATTGTCTGGGGGCTTTTTATGAGCTAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATATTCTGGCTCGAGTTAACTAcc                                                                                               <  2:720511/150‑1 (MQ=255)
                                                      ataatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAAcc                                                                                               <  2:1258539/150‑1 (MQ=255)
                                                       taatCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTAAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCa                                                                                              <  2:1601992/150‑1 (MQ=255)
                                                         atCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGCACCTCATAAACTGGCTCGAGATAGCAACCaaa                                                                                            >  1:1507760/1‑150 (MQ=255)
                                                          tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAt                                                                                           >  2:641175/1‑150 (MQ=255)
                                                          tCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAt                                                                                           >  1:1225375/1‑150 (MQ=255)
                                                          tCCATGTCCGCCACAGTTGCGCGAATTCAGCAGAACGAATGTCGAAGAGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAAt                                                                                           <  1:74885/150‑1 (MQ=255)
                                                               ttCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATa                                                                                      <  2:555423/150‑1 (MQ=255)
                                                                    ccACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAAtt                                                                                 >  1:1689010/1‑150 (MQ=255)
                                                                       caGTTGCGCCAATTCAGCAGCACGAATGACGATGCGCTGTTGGTTTATTGTCAGGGGGCGTTTGATGAGGAAGCATGGTTGCTCATTGACATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCTTATAATTggg                                                                              >  2:1438785/1‑150 (MQ=255)
                                                                       aaGTTGCGCCAATTAAGCAGAACGAATGACGAAGCGCCTTAGGTATATTGTCAGGTGGCGTTTGATGAGCTAGCATTGTTGCTCATTGCCCTTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTggg                                                                              <  2:290821/149‑1 (MQ=255)
                                                                               ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGTGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATgac                                                                      <  1:1443246/150‑3 (MQ=255)
                                                                               ccAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGc                                                                      <  1:1167128/150‑1 (MQ=255)
                                                                                cAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCa                                                                     <  1:996844/150‑1 (MQ=255)
                                                                                   ttCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa                                                                  <  2:1285771/150‑1 (MQ=255)
                                                                                   ttCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGaa                                                                  <  1:487589/150‑1 (MQ=255)
                                                                                     cagcagAACGAATGCCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGc                                                                <  1:311787/150‑1 (MQ=255)
                                                                                      agcagAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCg                                                               <  1:836320/150‑1 (MQ=255)
                                                                                          gAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTc                                                           >  2:325565/1‑150 (MQ=255)
                                                                                             cGACTGACGAAGCGCATTTGGTATATTTTCAGGTGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa                                                        <  2:502788/150‑1 (MQ=255)
                                                                                             cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa                                                        <  2:1490884/150‑1 (MQ=255)
                                                                                             cGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCa                                                        >  2:1079770/1‑150 (MQ=255)
                                                                                             cGAATGACGAAGCGCATTTGCTATATTGTAAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATTAACGGGCTCGAGATAAACACCAAATGCATAAAAGTTGGAAAATGGCAGAAGCGTTTCGCa                                                        >  1:1395897/1‑150 (MQ=255)
                                                                                              gAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCAt                                                       >  1:993596/1‑150 (MQ=255)
                                                                                               aaTGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATg                                                      <  1:1753007/150‑1 (MQ=255)
                                                                                                   aCGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGtc                                                  >  2:1127407/1‑149 (MQ=255)
                                                                                                     gAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGATCCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGTGAAAATGGCCCAAGCGTTTCGCATGCGtttt                                                >  2:871087/1‑150 (MQ=255)
                                                                                                       cgcgcAGTTGGACTATTGTCAGGGGGCGTATGATGAGCAACCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAACTGCATAAAATTGGGAAAATGGCAGAAGCGGTTCGCATGCGTTTTTg                                              <  2:792534/149‑1 (MQ=255)
                                                                                                       aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTg                                              <  1:1731327/150‑1 (MQ=255)
                                                                                                       aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTg                                              <  2:12377/150‑1 (MQ=255)
                                                                                                       aGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTg                                              >  1:1823279/1‑150 (MQ=255)
                                                                                                          gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGCAGCGTTTCGCATGCGTTTTTGAAt                                           >  1:314784/1‑150 (MQ=255)
                                                                                                          gcATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAAt                                           <  1:1277498/150‑1 (MQ=255)
                                                                                                                ggTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTatat                                     >  2:1568415/1‑150 (MQ=255)
                                                                                                                  tataTTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtatta                                   <  1:1254436/150‑1 (MQ=255)
                                                                                                                   atatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtattat                                  >  1:985764/1‑150 (MQ=255)
                                                                                                                    tatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGATCCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTAtattatt                                 >  1:496090/1‑149 (MQ=255)
                                                                                                                    tatTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAATTGCATAATATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATg                                 <  1:440572/150‑1 (MQ=255)
                                                                                                                      ttGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGaa                               <  2:207036/150‑1 (MQ=255)
                                                                                                                        gTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATa                             <  2:1689010/150‑1 (MQ=255)
                                                                                                                           aGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACa                          >  2:607383/1‑150 (MQ=255)
                                                                                                                            gggggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACAt                         <  1:659978/150‑1 (MQ=255)
                                                                                                                              gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAACCTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATAc                       <  2:527072/150‑1 (MQ=255)
                                                                                                                              gggCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATAc                       <  1:1322563/150‑1 (MQ=255)
                                                                                                                                gCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaa                     >  2:1520622/1‑150 (MQ=255)
                                                                                                                                gCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaa                     <  2:279002/150‑1 (MQ=255)
                                                                                                                                 cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa                    <  2:1351471/150‑1 (MQ=255)
                                                                                                                                 cGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACaaa                    >  2:644754/1‑150 (MQ=255)
                                                                                                                                   ggcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATa                  >  2:597143/5‑150 (MQ=255)
                                                                                                                                    ttGAAGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAACTGCATAAAAATGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATaa                 <  1:1492174/150‑1 (MQ=255)
                                                                                                                                    gcaatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATaa                 >  2:1795253/4‑150 (MQ=255)
                                                                                                                                     tgatgaGCAAGCATGGTTGCTCATTGTCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg                <  1:1770382/150‑1 (MQ=255)
                                                                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAg                >  1:329944/3‑150 (MQ=255)
                                                                                                                                     caatgaGCAACCATGCTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATAAAAATAAg                >  1:148185/3‑150 (MQ=255)
                                                                                                                                       atgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGaa              <  2:817718/150‑1 (MQ=255)
                                                                                                                                        tgaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTTTGAATAACATACAAATAAGaaa             >  1:950536/1‑150 (MQ=255)
                                                                                                                                         gaGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAAc            <  2:1506251/150‑1 (MQ=255)
                                                                                                                                          agcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACt           >  1:1318709/1‑150 (MQ=255)
                                                                                                                                           gcaagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTa          >  1:1667564/1‑150 (MQ=255)
                                                                                                                                            caagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTAt         >  1:108987/1‑150 (MQ=255)
                                                                                                                                             aagcaTGGTTGCTCATTGCCATTCTGACACCTGAACATCATAAACTGGCTCGAGATAACAACCACATGCATAAAATTGGGTAAATGGCAGAACCTTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAAGAAGAAACAATc        >  2:876628/1‑150 (MQ=255)
                                                                                                                                             aagcaTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATc        <  2:1225375/150‑1 (MQ=255)
                                                                                                                                                caTGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGc     <  2:1455386/150‑1 (MQ=255)
                                                                                                                                                   ggTTGCTCATTGCAATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTg  <  2:1135715/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTGGCCATTAATCAAATTTGAGCGAGTTGCTCATATTCATCTGGCAAATGAGAATAATCCATTTCCGCCACAGTTGCGCCAATTCAGCAGAACGAATGACGAAGCGCATTTGGTATATTGTCAGGGGGCGTTTGATGAGCAAGCATGGTTGCTCATTGCCATTCTGAAACCTGAACCTCATAAACTGGCTCGAGATAACAACCAAATGCATAAAATTGGGAAAATGGCAGAAGCGTTTCGCATGCGTTTTTGAATTTATATTATGAATAACATACAAATAAGAAACTATCAGCCTG  >  NC_000913/252446‑252742

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: