Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,856,257 G→C 6.1% L640L (CTG→CTC fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,856,2570GC 6.1% 426.2 / 13.4 139L640L (CTG→CTCfhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base G (65/65);  new base C (6/3);  total (71/68)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

Only 100 of 145 total aligned reads displayed.

AACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGC  >  NC_000913/2856108‑2856406
                                                                                                                                                     |                                                                                                                                                     
aaCATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGCGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCGGCCGCAACGGTTGTCGCCCTg                                                                                                                                                       >  2:1206089/1‑150 (MQ=255)
aaCATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTg                                                                                                                                                       <  2:871465/150‑1 (MQ=255)
     ggAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAggg                                                                                                                                                  <  2:27029/150‑1 (MQ=255)
     ggAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAggg                                                                                                                                                  <  1:629700/150‑1 (MQ=255)
     ggAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAggg                                                                                                                                                  <  1:1178895/150‑1 (MQ=255)
      gAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGGATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGc                                                                                                                                                 <  2:1195818/150‑1 (MQ=255)
          ggCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAg                                                                                                                                             >  2:841552/1‑150 (MQ=255)
            ccGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTCCAGGGCGacaac                                                                                                                                           <  2:430020/150‑5 (MQ=255)
            ccGGGTAACGTACGCGAACTGGAAAACGGCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTCCAGGGCGacaac                                                                                                                                           <  2:1803602/150‑5 (MQ=255)
             cGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACTCCGCCTGCCGCAACGGTTGTCGCCCTGGATGGCGAAGATg                                                                                                                                          >  2:1701002/1‑150 (MQ=255)
             cGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATg                                                                                                                                          <  1:660918/150‑1 (MQ=255)
               ggTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGaa                                                                                                                                        <  1:1954874/150‑1 (MQ=255)
                gTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAat                                                                                                                                       <  2:542868/150‑1 (MQ=255)
                 tAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAata                                                                                                                                      <  2:801888/150‑1 (MQ=255)
                         gcgAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGa                                                                                                                              >  2:1508539/1‑150 (MQ=255)
                            aaCTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATtg                                                                                                                           <  1:1817344/150‑1 (MQ=255)
                              cTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATtgtg                                                                                                                         <  2:1718337/150‑1 (MQ=255)
                                ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcg                                                                                                                       >  1:1915310/1‑150 (MQ=255)
                                ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcg                                                                                                                       >  1:626798/1‑150 (MQ=255)
                                ggAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTgcc                                                                                                                       >  1:1712167/1‑149 (MQ=255)
                                      gggcATTGTGCGCGAGGTATTGCTAACACGCGGGAACGTGCTGCAGCGGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGGCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCt                                                                                                                 <  1:281766/147‑1 (MQ=255)
                                      cggcATTGAGCGCGCGGTATTGCTACCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATTAATATCAGTTGATTGTGCGCGTGCt                                                                                                                 <  1:1320022/147‑1 (MQ=255)
                                      cGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCGGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCt                                                                                                                 >  1:1274767/1‑150 (MQ=255)
                                       gTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTg                                                                                                                <  2:408792/150‑1 (MQ=255)
                                         cATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTgaa                                                                                                              >  1:673230/1‑150 (MQ=255)
                                               gcgcgcGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAAc                                                                                                        <  2:503211/150‑1 (MQ=255)
                                                cgggcgGTATTGCTAACACGCGGAAACGTGCTGCAGCTGACATTGCCAGATATTGTTTTACCGGAACCTGCAAAGCCGCCTGCAGCAACGGATGTCGCCCTGGAGGGCGAAGATGAATATCAGTTCATTGTGAGCGTGCTGAAAGAAAcc                                                                                                       >  1:109678/4‑150 (MQ=255)
                                                cgcgcgGGATTGCTAACACGCGGTAACGTGCGGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAAcc                                                                                                       <  1:274268/150‑1 (MQ=255)
                                                 gagctgTATTGCTAAGACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAATGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCa                                                                                                      <  2:1293714/145‑1 (MQ=255)
                                                      gTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGc                                                                                                 <  2:725714/150‑1 (MQ=255)
                                                        aTTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGCCGt                                                                                               >  2:1965118/1‑150 (MQ=255)
                                                          tGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGTTTGTGCGCGTGCTGAATGAAACCAACGGCGTgg                                                                                             <  1:1206089/150‑1 (MQ=255)
                                                          tGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTgg                                                                                             >  1:1055177/1‑150 (MQ=255)
                                                          tGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTgg                                                                                             <  2:680301/150‑1 (MQ=255)
                                                          tGCTAACACGCGGTAACGTGCTGCAGCGGGCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTgg                                                                                             <  2:1818744/150‑1 (MQ=255)
                                                           gCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGt                                                                                            >  2:1575618/1‑150 (MQ=255)
                                                           gCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGt                                                                                            >  1:1489934/1‑150 (MQ=255)
                                                            cTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGtt                                                                                           >  1:1506525/1‑150 (MQ=255)
                                                              aaCACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGc                                                                                         >  2:794271/1‑150 (MQ=255)
                                                                cacGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCg                                                                                       >  1:722473/1‑150 (MQ=255)
                                                                  cgcgGTAACGTGCTGCAGCTGTCATTGCCAGATCTTGTTTTACCGGACCCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCggg                                                                                     <  1:1587804/150‑1 (MQ=255)
                                                                  cgcgGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCggg                                                                                     >  2:676059/1‑150 (MQ=255)
                                                                  cgcgGTAACGTGCTGCAGCGGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGGCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCggg                                                                                     <  1:1928303/150‑1 (MQ=255)
                                                                     ggTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCt                                                                                  <  2:936141/150‑1 (MQ=255)
                                                                      gTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTa                                                                                 >  1:1834088/1‑150 (MQ=255)
                                                                          cGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAgg                                                                             >  2:1388760/1‑150 (MQ=255)
                                                                                gcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCATTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgc                                                                       >  2:1745482/1‑150 (MQ=255)
                                                                                gcAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgc                                                                       >  1:45553/1‑150 (MQ=255)
                                                                                gcAGCTGTCAGTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTgc                                                                       >  2:1195768/1‑150 (MQ=255)
                                                                                   gCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCa                                                                    >  1:1256238/1‑150 (MQ=255)
                                                                                    cTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGGGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCaa                                                                   >  1:1901562/1‑150 (MQ=255)
                                                                                     tGTCATTGCCAGATATTGTTGTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGAGCGCGTGCTGAAAGAAACCAACGGCGTGGGTGCCGGGCCTAAAGGCGCTGCGCAAc                                                                  <  1:107768/150‑1 (MQ=255)
                                                                                      gTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACg                                                                 <  1:841552/150‑1 (MQ=255)
                                                                                      gTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACg                                                                 >  1:701655/1‑150 (MQ=255)
                                                                                       tCATTGACAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTGGTCGCCCTGGAGGGCGAAGATGAATCTCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGt                                                                <  1:113689/150‑1 (MQ=255)
                                                                                        cATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTc                                                               >  2:1413374/1‑150 (MQ=255)
                                                                                        cATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGCTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTc                                                               >  2:1130812/1‑150 (MQ=255)
                                                                                            gCCAGATATTGTTTTCCCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTggg                                                           <  1:190814/150‑1 (MQ=255)
                                                                                             ccAGATATTGTGTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGAGTGTGCGCGTGCTGAAAGACAACAACGCCGTGGTTGCCGCGCCTAAAGGCGATGAGCAACGTCTgggg                                                          >  2:246584/1‑150 (MQ=255)
                                                                                              cAGATATTGTTTTACCGGACCCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGc                                                         <  1:1321767/150‑1 (MQ=255)
                                                                                               aGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGAGCCTAAAGGCGCTGCGCAACGTCTGGGGCt                                                        >  2:1320576/1‑150 (MQ=255)
                                                                                                gATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCTCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTg                                                       >  1:350552/1‑150 (MQ=255)
                                                                                                gATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTg                                                       >  1:1509204/1‑150 (MQ=255)
                                                                                                 atatTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGa                                                      <  1:1136782/150‑1 (MQ=255)
                                                                                                      gTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGc                                                 >  2:1307493/1‑150 (MQ=255)
                                                                                                        tttACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCAc                                               >  1:1028044/1‑150 (MQ=255)
                                                                                                         ttACCGGAACCTGAAACGCCGCCTGCCGCATCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGACAGAAACCACCGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACTTCTGGGGCTGAAACGCACg                                              <  2:1035700/150‑1 (MQ=255)
                                                                                                           aCCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAc                                            >  1:511657/1‑150 (MQ=255)
                                                                                                             gggAACCTGAAACGCCGCCTGCCTCAACGGTTCTCGCCCTGGAGGGCGATGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGAccc                                          <  1:1668636/149‑1 (MQ=255)
                                                                                                              ggAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGCAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCAAGACCct                                         >  2:1249062/1‑150 (MQ=255)
                                                                                                                aaCCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCATTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgc                                       <  1:560605/150‑1 (MQ=255)
                                                                                                                  ccTGAACCGCCGCATGCCGCAACCGTACTCGCCCTGCCGGGCGAAGATGAATATCAGTTGCTTGGGCGCGTGCTGAAACCAACCAACGGCCTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctg                                     <  2:315591/150‑1 (MQ=255)
                                                                                                                  ccTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctg                                     >  2:642368/1‑150 (MQ=255)
                                                                                                                   cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCctgctga                                    >  1:485238/1‑149 (MQ=255)
                                                                                                                   cTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGt                                    >  2:1637805/1‑150 (MQ=255)
                                                                                                                    tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTc                                   <  1:1411446/150‑1 (MQ=255)
                                                                                                                    tGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTc                                   <  2:209998/150‑1 (MQ=255)
                                                                                                                       aaCGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACg                                <  2:573596/150‑1 (MQ=255)
                                                                                                                       aaCGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGACGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACg                                >  1:542088/1‑150 (MQ=255)
                                                                                                                        aCGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACgg                               <  2:1274767/150‑1 (MQ=255)
                                                                                                                        aCGCCGCCGGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACgg                               <  1:1387524/150‑1 (MQ=255)
                                                                                                                            cgccTGCCGCAACGGGTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGa                           <  2:548985/150‑1 (MQ=255)
                                                                                                                                    gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTg                   >  1:1845969/1‑150 (MQ=255)
                                                                                                                                    gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTg                   >  1:33154/1‑150 (MQ=255)
                                                                                                                                    gCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTg                   >  2:71961/1‑150 (MQ=255)
                                                                                                                                     caccGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTgg                  <  1:1114926/147‑1 (MQ=255)
                                                                                                                                     cAACGGTTGTCGCCCTGGAGGGCGATGATGAATATCAGTTGATTGTGCTCGTGCTGAATGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGAGCAACGTCTGGGGCTGAAACTCACGACCCTGCTGTCACGGATGAAGCGGCTgg                  <  2:109678/150‑1 (MQ=255)
                                                                                                                                       aCGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGa                <  2:701655/150‑1 (MQ=255)
                                                                                                                                        cGGGTGTCGCCCTGGAGGGCGATGATGAATATCAGTGGATTGTGCGCGGGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGaa               <  1:202977/150‑1 (MQ=255)
                                                                                                                                         gggtGTCGCCCTGGAGGGCGAAGATGAATAGCAGTTGATTGTGCGCGAGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAt              <  1:1743030/147‑1 (MQ=255)
                                                                                                                                         catcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGAGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAAt              >  2:303027/6‑150 (MQ=255)
                                                                                                                                          atcttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGGGTGCTGAAAGAAACCAACGTCGTGTTTGCCGGGCGTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCGAGCCCCTGCTGTCAAGGATGAAGCGGCTGGGAAtt             >  1:1249229/5‑150 (MQ=255)
                                                                                                                                           ttGGCGCCCTGGATGGCGAAGATGTATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTCAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGACGCGGCTGGGAATTg            <  2:162826/150‑1 (MQ=255)
                                                                                                                                            ggTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGGTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa           <  1:1910442/149‑1 (MQ=255)
                                                                                                                                            cttCGCCCTCCAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGa           >  2:406310/3‑150 (MQ=255)
                                                                                                                                             gTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAACCCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGAt          >  2:1286838/1‑150 (MQ=255)
                                                                                                                                                   cTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCg    >  2:844955/1‑150 (MQ=255)
                                                                                                                                                    tGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAACTCgg   >  1:785942/1‑150 (MQ=255)
                                                                                                                                                    tGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGAGGCTGGGAATTGATAAATCgg   <  2:496437/150‑1 (MQ=255)
                                                                                                                                                     ggAGGGGGAAGATGTATTTCAGTTGATTTTGCGCGTGCTGAAATAAACCAACGGCGTGGTTGCCGGGGCTATAGGGGCTGCGCAAAGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGc  <  1:1425081/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACATGGAGTGGCCGGGTAACGTACGCGAACTGGAAAACGTCATTGAGCGCGCGGTATTGCTAACACGCGGTAACGTGCTGCAGCTGTCATTGCCAGATATTGTTTTACCGGAACCTGAAACGCCGCCTGCCGCAACGGTTGTCGCCCTGGAGGGCGAAGATGAATATCAGTTGATTGTGCGCGTGCTGAAAGAAACCAACGGCGTGGTTGCCGGGCCTAAAGGCGCTGCGCAACGTCTGGGGCTGAAACGCACGACCCTGCTGTCACGGATGAAGCGGCTGGGAATTGATAAATCGGC  >  NC_000913/2856108‑2856406

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