Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 838,985 A→T 9.9% L75M (TTG→ATG)  ybiX ← Fe(II)‑dependent oxygenase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913838,9850AT 9.9% 308.9 / 18.3 120L75M (TTG→ATG) ybiXFe(II)‑dependent oxygenase superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (60/47);  new base T (2/10);  total (62/58)
Fisher's exact test for biased strand distribution p-value = 1.32e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.48e-01

Only 100 of 120 total aligned reads displayed.

TGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCA  >  NC_000913/838839‑839132
                                                                                                                                                  |                                                                                                                                                   
tgGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGc                                                                                                                                                  >  2:707816/1‑150 (MQ=255)
 ggATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGTGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGcc                                                                                                                                                 >  2:1082378/1‑150 (MQ=255)
  gATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCcg                                                                                                                                                >  1:234777/1‑150 (MQ=255)
     cGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCgg                                                                                                                                             <  1:1218353/150‑1 (MQ=255)
      gCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCTTCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgta                                                                                                                                            <  2:1719651/150‑3 (MQ=255)
      gCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGc                                                                                                                                            <  1:1679851/150‑1 (MQ=255)
      gCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGc                                                                                                                                            <  1:842446/150‑1 (MQ=255)
       cTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCa                                                                                                                                           >  2:62566/1‑150 (MQ=255)
       cTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCa                                                                                                                                           <  1:1083926/150‑1 (MQ=255)
       cTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCa                                                                                                                                           >  2:448381/1‑150 (MQ=255)
        ttAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtacc                                                                                                                                          <  2:1381454/150‑5 (MQ=255)
         tAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtaccc                                                                                                                                         <  2:762715/150‑6 (MQ=255)
         tAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtaccc                                                                                                                                         <  1:1486849/150‑6 (MQ=255)
         tAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtaccc                                                                                                                                         <  2:1004048/150‑6 (MQ=255)
         tAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtaccc                                                                                                                                         <  1:1603244/150‑6 (MQ=255)
          aaaaaCAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtaccct                                                                                                                                        <  2:191377/150‑7 (MQ=255)
          aaaaaCAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGGGGAAAGGGTACGCGGCAAGGCCGCGGCAAAg                                                                                                                                        >  2:1128508/1‑150 (MQ=255)
           aaaaCAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtaccctt                                                                                                                                       <  1:402162/150‑8 (MQ=255)
               cAGCGTTGCAGAAACGTCAGTACGCATGCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGCAAACCTTAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTTCGCGGCAAGGCCGCGGCAAAGAATaa                                                                                                                                   <  2:353936/150‑1 (MQ=255)
                 gCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAg                                                                                                                                 >  2:744084/1‑150 (MQ=255)
                     ggCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATg                                                                                                                             >  2:1660200/1‑150 (MQ=255)
                        aGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTg                                                                                                                          >  2:1528208/1‑150 (MQ=255)
                           aaGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGtt                                                                                                                       >  2:1393302/1‑150 (MQ=255)
                             ggTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGATAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTaa                                                                                                                     >  2:131572/1‑150 (MQ=255)
                             ggTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTaa                                                                                                                     >  1:1663746/1‑150 (MQ=255)
                              gTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAAc                                                                                                                    <  1:1393302/150‑1 (MQ=255)
                                cAGTACACATCCAGCCGTTTTGTGGATGGCTGCGTACAGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACcg                                                                                                                  <  1:1530218/150‑1 (MQ=255)
                                  gTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGATTCATTGTTCTGAGAGCGATTAAACAGCGGCGGGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACcgcg                                                                                                                >  1:71262/1‑150 (MQ=255)
                                     cGCATCCAGCCGTTTTGTGTATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACATCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTc                                                                                                             <  1:843150/150‑1 (MQ=255)
                                       cATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATAAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAg                                                                                                           >  2:1022340/1‑150 (MQ=255)
                                        aTCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGc                                                                                                          <  1:8583/150‑1 (MQ=255)
                                          cctgCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGTTTTAACCGCGGCGTGGAAAGGGTACGCGGCATGGCCGCGGCAAAGAATAAAGCATGTTGGTTTACCGCGATCAGCAc                                                                                                        >  2:102651/4‑150 (MQ=255)
                                           cAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCAcc                                                                                                       >  2:321115/1‑150 (MQ=255)
                                           cAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGAAcc                                                                                                       >  2:329635/1‑150 (MQ=255)
                                            aGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTACCCGCGTTCAGCAccg                                                                                                      >  2:56214/1‑149 (MQ=255)
                                             gCCGTTTTGTGGATGGCTGCGTACTGCGCCGTCCACATGGAAACCATAGGTTTCATTGTTCTGATTGCGATTAAACAGCGGCGTGGAAAGGGTCCGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTc                                                                                                     >  1:688903/1‑150 (MQ=255)
                                             gCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTc                                                                                                     >  1:574984/1‑150 (MQ=255)
                                                 ccatgtgGATGGCGCAGTACGCAGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCAttt                                                                                                 >  1:1343290/4‑150 (MQ=255)
                                                 ccatgtgGATGGCGCAGTACGCAGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCAttt                                                                                                 >  1:1517420/4‑150 (MQ=255)
                                                  tttGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCAtttt                                                                                                <  2:75144/150‑1 (MQ=255)
                                                   ttGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTg                                                                                               >  1:1015559/1‑150 (MQ=255)
                                                   atgtgGATGGCGCAGTACGCAGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTg                                                                                               >  2:530012/2‑150 (MQ=255)
                                                     gtgGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCa                                                                                             >  2:1390152/1‑150 (MQ=255)
                                                     gtgGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCa                                                                                             >  1:140027/1‑150 (MQ=255)
                                                     gtgGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGGTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCa                                                                                             <  2:850604/150‑1 (MQ=255)
                                                      tgGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCCTTGCcgcgtaccc                                                                                                                                         <  1:1576899/105‑6 (MQ=255)
                                                       ggATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAgg                                                                                           >  2:460610/1‑150 (MQ=255)
                                                           ggCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGggcgg                                                                                       >  2:303109/1‑150 (MQ=255)
                                                             cTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGgcggcg                                                                                     >  1:1454488/1‑150 (MQ=255)
                                                             cTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAAGAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGgcggcg                                                                                     >  2:1104999/1‑150 (MQ=255)
                                                                 gTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGAGGCGTGGAAAGGGTACGCTGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGtata                                                                                 <  2:1279395/150‑1 (MQ=255)
                                                                   aCTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAAc                                                                               <  1:1276748/150‑1 (MQ=255)
                                                                    cTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACg                                                                              <  2:1305251/150‑1 (MQ=255)
                                                                    cTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACg                                                                              <  1:434996/150‑1 (MQ=255)
                                                                      gcgcCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTTAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtg                                                                            >  1:278172/1‑150 (MQ=255)
                                                                        gcCATCCACATGGAACCCATAGGTATCATTGTTCTGATCGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtgct                                                                          <  2:198637/150‑1 (MQ=255)
                                                                          cATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGtgctgc                                                                        <  2:1373345/150‑1 (MQ=255)
                                                                            tCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGCTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTAGAACGTGCTGCgg                                                                      >  2:179979/1‑150 (MQ=255)
                                                                              cacaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGTTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCGGCGCgt                                                                    >  2:842615/1‑150 (MQ=255)
                                                                              cacaTGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGgt                                                                    >  2:880390/1‑150 (MQ=255)
                                                                                  ttgacaCCATAGGTTTCATTGTTCTTATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTTGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCg                                                                <  2:136309/145‑1 (MQ=255)
                                                                                  tGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCg                                                                >  1:625064/1‑150 (MQ=255)
                                                                                  tGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCg                                                                >  2:1491178/1‑150 (MQ=255)
                                                                                  tGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCg                                                                <  1:538082/150‑1 (MQ=255)
                                                                                  tGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCg                                                                >  2:1683530/1‑150 (MQ=255)
                                                                                       aCCATAGGTTTCATTGGTCTGAAAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCACGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTg                                                           <  1:277641/150‑1 (MQ=255)
                                                                                         cATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGGTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGtt                                                         <  1:1042398/150‑1 (MQ=255)
                                                                                            aGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGAt                                                      <  1:131572/150‑1 (MQ=255)
                                                                                            aGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGAt                                                      >  1:889127/1‑150 (MQ=255)
                                                                                                   aTTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCtt                                               <  2:690771/150‑1 (MQ=255)
                                                                                                      gTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAAc                                            >  2:935324/1‑150 (MQ=255)
                                                                                                      gTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAAc                                            >  1:499101/1‑150 (MQ=255)
                                                                                                         cTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTg                                         <  2:520980/150‑1 (MQ=255)
                                                                                                         cTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTg                                         <  2:478135/150‑1 (MQ=255)
                                                                                                           gATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTgcg                                       <  2:1663746/150‑1 (MQ=255)
                                                                                                               gCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCAcc                                   >  2:1530005/1‑150 (MQ=255)
                                                                                                                  aTTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTtg                                >  1:126794/1‑150 (MQ=255)
                                                                                                                    tAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTTACCGCGTTCAGCACCTCATTTTGCTGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTtgtg                              <  1:1070310/150‑1 (MQ=255)
                                                                                                                    tAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTtgtg                              <  2:814189/150‑1 (MQ=255)
                                                                                                                    tAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTGCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTtgtg                              <  2:946093/150‑1 (MQ=255)
                                                                                                                     aaaCAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTgc                             <  1:460610/150‑1 (MQ=255)
                                                                                                                           cggcgTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCggt                       <  2:574984/150‑1 (MQ=255)
                                                                                                                              cgTGGCAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCAGGTTGGTTAACCGCGTTCAGCACCTCAATTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCggtggt                    <  2:237357/150‑1 (MQ=255)
                                                                                                                              cgTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCggtggt                    <  2:234777/150‑1 (MQ=255)
                                                                                                                                 ggAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGCTTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGAc                 >  2:1060952/1‑150 (MQ=255)
                                                                                                                                  gAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGAcg                >  2:1540271/1‑150 (MQ=255)
                                                                                                                                   aaaGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGAcgc               >  1:1468215/1‑150 (MQ=255)
                                                                                                                                     aGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGt             >  2:1310531/1‑150 (MQ=255)
                                                                                                                                      gggTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGtc            <  2:773304/150‑1 (MQ=255)
                                                                                                                                         tACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTGGACCTGTTGATTGTGCTTAACTTGCGCACCTTGTGCGCCGGTGGTGCCGCGtccat         >  1:829441/1‑150 (MQ=255)
                                                                                                                                            gcgGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGtccatcca      <  1:448381/150‑1 (MQ=255)
                                                                                                                                            gcgGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTCCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTCTGCGCCGGTGGTGACGCGtccctccc      >  1:71298/1‑149 (MQ=255)
                                                                                                                                             cgGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTc                                                                 >  1:1630520/1‑90 (MQ=255)
                                                                                                                                             cgGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTc                                                                 <  2:1630520/90‑1 (MQ=255)
                                                                                                                                              ggCCAGGCCGCGGCAAAGAATAAAGCATGTTGGTTATCCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGGGCTGCGGGTGGCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCAcc    <  2:356179/150‑1 (MQ=255)
                                                                                                                                              ggCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCAcc    >  1:596677/1‑150 (MQ=255)
                                                                                                                                              ggCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCAcc    <  1:492138/150‑1 (MQ=255)
                                                                                                                                               gCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCAccc   <  2:150896/150‑1 (MQ=255)
                                                                                                                                               gCAAGCCCGCGGAAAAGAATAAAGCATGTTGGTTAACCGCGTACAGCCCCTCATTTTGCAGGGCGGCGTATAACGTCCTGCGGGTGTCGACCTGTAGATTGTTCTTAACTTGCGCACCTTGAGCGCCTGTGCTCCCGCGTCCATCCAccc   >  1:369793/1‑150 (MQ=255)
                                                                                                                                                cAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCa  <  1:1082378/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TGGATCGCTTAAAAACAGCGTGGCAGAAAGGTCAGTACGCATCCAGCCGTTTTGTGGATGGCTGCGTACTGCGCCATCCACATGGAAACCATAGGTTTCATTGTTCTGATAGCGATTAAACAGCGGCGTGGAAAGGGTACGCGGCAAGGCCGCGGCAAAGAATAAAGCATGTTGGTTAACCGCGTTCAGCACCTCATTTTGCAGGGCGGCGTATAACGTGCTGCGGGTGTCGACCTGTTGATTGTTCTTAACTTGCGCACCTTGTGCGCCGGTGGTGACGCGTCCATCCACCCA  >  NC_000913/838839‑839132

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: ‑