Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,088,193 C→A intergenic (+333/‑91) metK → / → galP S‑adenosylmethionine synthetase/D‑galactose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,088,1930CA98.5% 414.1 / ‑3.0 132intergenic (+333/‑91)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (67/62);  total (67/65)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

Only 100 of 132 total aligned reads displayed.

TTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCT  >  NC_000913/3088050‑3088338
                                                                                                                                               |                                                                                                                                                 
    tattattatTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACt                                                                                                                                          <  2:83031/149‑1 (MQ=255)
     attattattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTg                                                                                                                                         >  1:1167467/1‑149 (MQ=255)
      ttattattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGa                                                                                                                                        >  1:232332/1‑149 (MQ=255)
       tattattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATATATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGAt                                                                                                                                       <  1:807850/149‑1 (MQ=255)
        attattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAATCTCTTACACtaaca                                                                                                                                      <  2:499438/149‑5 (MQ=255)
        attattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAATCTCTTACACtaaca                                                                                                                                      <  2:386600/149‑5 (MQ=255)
        attattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGAtg                                                                                                                                      >  1:1412023/1‑149 (MQ=255)
        attattCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGATTAATCTATGATGATATAACTCAATTATTTTCATTCACTTAAATCATAACTAATATAAATGTTAGTGTAATATATTACACTGAtt                                                                                                                                      >  1:1283679/1‑148 (MQ=255)
          tattCGGATTTCACAGTTTTTACATTTCTTTTCAGTAAAGTATTAATTGCATATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGAtgtg                                                                                                                                    <  2:1758035/149‑1 (MQ=255)
            ttaggTTTTCACAGTTTTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCTCTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTTAAATTATGTGAt                                                                                                                                  >  2:1573095/4‑149 (MQ=255)
            ttCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGAt                                                                                                                                  >  2:680123/1‑149 (MQ=255)
            ttCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGAt                                                                                                                                  >  1:58783/1‑149 (MQ=255)
             tCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGAtt                                                                                                                                 <  1:899456/149‑1 (MQ=255)
              cGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATATCAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGAttt                                                                                                                                <  2:1697645/149‑1 (MQ=255)
                 ttttCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATTACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGaagcaaa                                                                                                                             >  1:264647/1‑143 (MQ=255)
                 ttttCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCt                                                                                                                             <  2:413546/149‑1 (MQ=255)
                  tttCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTTCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACATAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGaagcaaat                                                                                                                            <  1:945129/149‑8 (MQ=255)
                   ttCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTc                                                                                                                           <  2:212382/149‑1 (MQ=255)
                     cacaGGTGGTACATTTCTTTTCAGTACAGTCTTACTTCCAGCTAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTcac                                                                                                                         <  1:1342195/149‑1 (MQ=255)
                      acaGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATTAATGTTAGTGTAATAGATTACACTGATGTGATTTGCTTcaca                                                                                                                        >  2:1546273/1‑149 (MQ=255)
                      acaGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAATAGATAACACTGATGTGATTTGCTTcaca                                                                                                                        >  2:978702/1‑149 (MQ=255)
                      acaGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTcaca                                                                                                                        >  2:1541368/1‑149 (MQ=255)
                      acaGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTcaca                                                                                                                        >  2:1411779/1‑149 (MQ=255)
                           tgttACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCttt                                                                                                                   <  1:1647705/149‑1 (MQ=255)
                             ttACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGGACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCttttt                                                                                                                 <  1:1426797/149‑1 (MQ=255)
                             ttACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCttttt                                                                                                                 >  2:1331599/1‑149 (MQ=255)
                              tACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTGCACATCTTTTTa                                                                                                                <  2:814021/149‑1 (MQ=255)
                                 aTTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTAcgt                                                                                                             >  1:562510/1‑149 (MQ=255)
                                 aTTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTAcgt                                                                                                             >  1:1059252/1‑149 (MQ=255)
                                  tttCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTAcgtc                                                                                                            <  1:1025880/149‑1 (MQ=255)
                                   ttCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTAcgtcg                                                                                                           <  2:1206289/149‑1 (MQ=255)
                                      ttttCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTTATTTGCTTCACATCTTTTTACGTCTTAc                                                                                                        >  1:1212146/1‑149 (MQ=255)
                                      ttttCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACACAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTAc                                                                                                        >  2:15990/1‑149 (MQ=255)
                                      ttttCAGTAAAGTCTTAATTGAAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCTCTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTTATTTGCTTCAAATCTTTTTACGTCGTAc                                                                                                        >  1:1186658/1‑149 (MQ=255)
                                       tttCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACt                                                                                                       >  2:556240/1‑149 (MQ=255)
                                         tCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCa                                                                                                     <  1:1090394/149‑1 (MQ=255)
                                           aGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCAcc                                                                                                   >  2:648024/1‑149 (MQ=255)
                                            gTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCt                                                                                                  <  1:673220/149‑1 (MQ=255)
                                            gTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCAAGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCt                                                                                                  >  2:304249/1‑149 (MQ=255)
                                             tAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTa                                                                                                 >  1:178485/1‑149 (MQ=255)
                                                aGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCt                                                                                              <  1:878860/149‑1 (MQ=255)
                                                aGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCt                                                                                              <  1:761169/149‑1 (MQ=255)
                                                  tCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTa                                                                                            >  2:1121658/1‑149 (MQ=255)
                                                  tCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTa                                                                                            >  1:1532527/1‑149 (MQ=255)
                                                   cTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTaa                                                                                           <  1:209723/149‑1 (MQ=255)
                                                    ttAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAAt                                                                                          >  1:485878/1‑149 (MQ=255)
                                                     tAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAAtt                                                                                         <  2:13090/149‑1 (MQ=255)
                                                       aTTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTATGTCGTACTCACCTATCTTAATTca                                                                                       <  1:761324/149‑1 (MQ=255)
                                                            aGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTCACGTCGTACTCACATATCTTAATTCACAATa                                                                                  >  1:547376/1‑149 (MQ=255)
                                                                   aGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATa                                                                           >  2:1559596/1‑149 (MQ=255)
                                                                    tcGTTTAATCTATTATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTATTGTAAGAGATTAAACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATaa                                                                          <  1:325954/148‑1 (MQ=255)
                                                                    gCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATaa                                                                          <  1:599074/149‑1 (MQ=255)
                                                                     cGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAAc                                                                         <  1:1331599/149‑1 (MQ=255)
                                                                       tttAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTTGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACTTAGTACTCACCTATCTTAATTCACAATATAAAATAACCa                                                                       >  1:1008560/1‑149 (MQ=255)
                                                                       tttAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCa                                                                       >  1:861102/1‑149 (MQ=255)
                                                                           aTCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATAtt                                                                   >  2:1020758/1‑149 (MQ=255)
                                                                           aTCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATAtt                                                                   >  2:343575/1‑149 (MQ=255)
                                                                              tATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGa                                                                >  2:972936/1‑149 (MQ=255)
                                                                              tATGATGATATAACACAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGa                                                                <  1:15990/149‑1 (MQ=255)
                                                                               atgatgATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAg                                                               <  2:477280/149‑1 (MQ=255)
                                                                                tgatgaTATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAgg                                                              <  2:39684/149‑1 (MQ=255)
                                                                                 gatgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAggg                                                             >  2:1046476/1‑149 (MQ=255)
                                                                                 gatgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAggg                                                             >  1:1362577/1‑149 (MQ=255)
                                                                                 gatgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAggg                                                             >  1:382374/1‑149 (MQ=255)
                                                                                 gatgatATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAggg                                                             >  1:1449751/1‑149 (MQ=255)
                                                                                     atatAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGcatc                                                         >  2:351691/1‑149 (MQ=255)
                                                                                      tataACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGcatcat                                                       >  1:32452/1‑149 (MQ=255)
                                                                                        taACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATg                                                      >  2:181629/1‑149 (MQ=255)
                                                                                         aaCTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGc                                                     >  2:679168/1‑149 (MQ=255)
                                                                                           atCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCt                                                   <  2:1441262/148‑1 (MQ=255)
                                                                                            tCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGATTCACATCTTTTTACGTCGTACTAACATATCTTAATTCACAATAAAAAATAACCATATTTGAGGGCATAATGCCTg                                                  >  1:1040560/1‑149 (MQ=255)
                                                                                              aaTTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGAc                                                <  2:1710164/149‑1 (MQ=255)
                                                                                              aaTTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGAc                                                >  1:44751/1‑149 (MQ=255)
                                                                                                  attTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATTAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTa                                            <  1:624380/149‑1 (MQ=255)
                                                                                                    tttCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTaaa                                          <  1:679168/149‑1 (MQ=255)
                                                                                                        aTGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAAc                                      >  2:1402033/1‑149 (MQ=255)
                                                                                                         tGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACa                                     <  2:572621/149‑1 (MQ=255)
                                                                                                         tGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATAATGCCTGACGCGAATAAAcg                                     >  2:1649456/1‑148 (MQ=255)
                                                                                                           cACTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTAACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAgg                                   >  1:126728/1‑149 (MQ=255)
                                                                                                            aCTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAggg                                  >  1:1462742/1‑149 (MQ=255)
                                                                                                             cTTAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATACACAATAAAAAATAACCATATTGGAGGGCATCATGCCAAACGATAAAAAACAggtg                                 >  1:995788/1‑147 (MQ=255)
                                                                                                              ttAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAg                                                               <  2:1103776/118‑1 (MQ=255)
                                                                                                              ttAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAg                                                               >  1:1103776/1‑118 (MQ=255)
                                                                                                               tAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCg                               <  2:1495874/149‑1 (MQ=255)
                                                                                                               tAAATCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCGGACGCTAAAAAACAGGGGCg                               <  2:851851/149‑1 (MQ=255)
                                                                                                                   tCATAACTAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCa                           >  2:905430/1‑149 (MQ=255)
                                                                                                                      aaaCTAAGATAAATTTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAAc                        <  2:1212146/148‑1 (MQ=255)
                                                                                                                          tAAGATAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAgg                    <  2:485878/149‑1 (MQ=255)
                                                                                                                              aTAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAAt                <  1:1020758/149‑1 (MQ=255)
                                                                                                                              aTAAATGTTAGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAAt                <  1:343575/149‑1 (MQ=255)
                                                                                                                                      tAGTTTAAGAGATGACACTGATGTGATTTGCATCACATCTTTTTACGTCGTACTCACCTATCTTCATTCACCATAAACAATCCCCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGGCAAACAAGGCAATGACGtttt        <  1:89679/149‑1 (MQ=255)
                                                                                                                                       aGTGTAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGttttt       <  2:1532527/149‑1 (MQ=255)
                                                                                                                                        gtgtAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTc      <  1:1411779/149‑1 (MQ=255)
                                                                                                                                        gtgtAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTc      <  2:562510/149‑1 (MQ=255)
                                                                                                                                        gtgtAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTc      <  2:126728/149‑1 (MQ=255)
                                                                                                                                        gtgtAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTc      <  2:178485/149‑1 (MQ=255)
                                                                                                                                         tgtAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCg     <  1:1559596/149‑1 (MQ=255)
                                                                                                                                           tAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTc   >  1:308540/1‑149 (MQ=255)
                                                                                                                                           tAAGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTc   >  1:1199910/1‑149 (MQ=255)
                                                                                                                                            aaGAGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCt  <  1:680123/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
TTCATATTATTATTCGGTTTTCACAGTTGTTACATTTCTTTTCAGTAAAGTCTTAATTGCAGATAACAGCGTTTAATCTATGATGATATAACTCAATTATTTTCATGCACTTAAATCATAACTAAGATAAATGTTAGTGTAAGCGATTACACTGATGTGATTTGCTTCACATCTTTTTACGTCGTACTCACCTATCTTAATTCACAATAAAAAATAACCATATTGGAGGGCATCATGCCTGACGCTAAAAAACAGGGGCGGTCAAACAAGGCAATGACGTTTTTCGTCT  >  NC_000913/3088050‑3088338

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: