Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,729 T→A 6.5% K106* (AAA→TAA)  ssuB ← aliphatic sulfonate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,7290TA 6.5% 376.8 / 11.8 107K106* (AAA→TAA) ssuBaliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (40/60);  new base A (2/5);  total (42/65)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

Only 100 of 113 total aligned reads displayed.

AACAATAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCA  >  NC_000913/993581‑993875
                                                                                                                                                    |                                                                                                                                                  
aacaaTAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCttt                                                                                                                                                    <  1:914381/149‑1 (MQ=255)
aacaaTAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCttt                                                                                                                                                    <  1:659767/149‑1 (MQ=255)
  caaTAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTaa                                                                                                                                                  <  4:163161/149‑1 (MQ=255)
  caaTAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTaa                                                                                                                                                  <  1:395903/149‑1 (MQ=255)
  caaTAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTaa                                                                                                                                                  <  1:301038/149‑1 (MQ=255)
  caaTAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATGCTCCAGTCCTACCGCAGCCAACGCTCGACGTGGGGCATCGCGCCACTGGCCTTTaa                                                                                                                                                  >  4:644346/1‑149 (MQ=255)
   aataatCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAg                                                                                                                                                 >  1:903533/1‑149 (MQ=255)
   aataatCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAg                                                                                                                                                 >  2:587288/1‑149 (MQ=255)
    ataatCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAgg                                                                                                                                                <  3:344562/149‑1 (MQ=255)
    ataatCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAgg                                                                                                                                                <  2:473101/149‑1 (MQ=255)
     taatCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGc                                                                                                                                               <  4:339641/149‑1 (MQ=255)
     taatCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGc                                                                                                                                               <  1:198076/149‑1 (MQ=255)
       atCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCCCTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCt                                                                                                                                             <  3:890652/149‑1 (MQ=255)
         cccGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGGCTTTAAGGCCTaa                                                                                                                                           <  1:320227/149‑1 (MQ=255)
          ccGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAAc                                                                                                                                          <  4:439309/149‑1 (MQ=255)
           cGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAAcc                                                                                                                                         <  4:156702/149‑1 (MQ=255)
            ggTCGATGAATCCACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTAAAGGccagtggc                                                                                                                                        <  4:233887/149‑7 (MQ=255)
            ggTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTAAAGGccagtggc                                                                                                                                        <  4:33136/149‑7 (MQ=255)
            ggTCGATAAATCAACACTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCACAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTAAAGGccagtggc                                                                                                                                        <  2:273055/149‑7 (MQ=255)
              tCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCaa                                                                                                                                      <  4:2494/149‑1 (MQ=255)
               cGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAAc                                                                                                                                     <  2:698838/149‑1 (MQ=255)
               cGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAAc                                                                                                                                     <  4:483570/149‑1 (MQ=255)
                gATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACg                                                                                                                                    >  4:541746/1‑149 (MQ=255)
                    aaTCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTat                                                                                                                                <  3:580712/149‑1 (MQ=255)
                    aaTCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTat                                                                                                                                <  4:5531/149‑1 (MQ=255)
                     aTCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTatc                                                                                                                               <  2:814105/149‑1 (MQ=255)
                        aaCGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTatcaat                                                                                                                            >  3:712577/1‑149 (MQ=255)
                            ctctTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCAcc                                                                                                                        <  1:599787/149‑1 (MQ=255)
                             tctTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCg                                                                                                                       >  2:340504/1‑149 (MQ=255)
                                       gcCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACgg                                                                                                             <  1:551110/149‑1 (MQ=255)
                                         cACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTa                                                                                                           <  2:903533/149‑1 (MQ=255)
                                            tCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCa                                                                                                        >  2:822808/1‑149 (MQ=255)
                                            tCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCa                                                                                                        <  1:624417/149‑1 (MQ=255)
                                             cGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAg                                                                                                       <  1:838802/149‑1 (MQ=255)
                                              gCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGa                                                                                                      <  2:704417/149‑1 (MQ=255)
                                                 gTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGAcgc                                                                                                   <  2:559452/149‑1 (MQ=255)
                                                      tGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATc                                                                                              >  3:219717/1‑149 (MQ=255)
                                                      tGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATc                                                                                              <  4:824686/149‑1 (MQ=255)
                                                        cccGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCCACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCt                                                                                            <  1:557281/149‑1 (MQ=255)
                                                         ccGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCAc                                                                                                               >  3:123675/1‑129 (MQ=255)
                                                         ccGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCAc                                                                                                               <  4:123675/129‑1 (MQ=255)
                                                         ccGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTg                                                                                           <  2:669533/149‑1 (MQ=255)
                                                          cGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGa                                                                                          <  2:727590/149‑1 (MQ=255)
                                                            ccAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCTACGGTAGCAGACGCGCATCCTGaaa                                                                                        >  3:605128/1‑149 (MQ=255)
                                                                aaaGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAAcatc                                                                                    <  1:345505/149‑1 (MQ=255)
                                                                    tGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTTAAACATCATTc                                                                                <  4:652308/149‑1 (MQ=255)
                                                                    tGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTc                                                                                <  3:704659/149‑1 (MQ=255)
                                                                    tGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTc                                                                                <  3:401029/149‑1 (MQ=255)
                                                                      cGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTcgc                                                                              <  3:650871/149‑1 (MQ=255)
                                                                        gCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAATATCATTCGCgt                                                                            <  3:541746/149‑1 (MQ=255)
                                                                        gCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCgt                                                                            <  2:7020/149‑1 (MQ=255)
                                                                        gCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCgt                                                                            <  4:872690/149‑1 (MQ=255)
                                                                           ggCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTc                                                                         <  4:505346/149‑1 (MQ=255)
                                                                              cATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTc                                                                      <  4:25227/149‑1 (MQ=255)
                                                                              cATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTc                                                                      <  1:261874/149‑1 (MQ=255)
                                                                                ttCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCt                                                                    <  4:712577/149‑1 (MQ=255)
                                                                                ttCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCt                                                                    <  3:426051/149‑1 (MQ=255)
                                                                                    ccTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAAt                                                                <  3:546717/149‑1 (MQ=255)
                                                                                     cTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAAtt                                                               >  2:419940/1‑149 (MQ=255)
                                                                                      tGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAAttt                                                              >  3:733913/1‑149 (MQ=255)
                                                                                           gATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTcagcc                                                         >  1:704226/1‑149 (MQ=255)
                                                                                           gATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTcagcc                                                         >  1:448313/1‑149 (MQ=255)
                                                                                                tcCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAAt                                                                <  4:195558/137‑1 (MQ=255)
                                                                                                tcCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCgg                                                    >  2:529188/1‑149 (MQ=255)
                                                                                                tcCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCgg                                                    >  2:259872/1‑149 (MQ=255)
                                                                                                tcCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCgg                                                    >  3:832467/1‑149 (MQ=255)
                                                                                                tcCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCgg                                                    >  1:220448/1‑149 (MQ=255)
                                                                                                tcCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCgg                                                    >  1:111620/1‑149 (MQ=255)
                                                                                                   aGTCCTACCGCAGCCAACGGTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCTCCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGgtgt                                                 <  3:725287/149‑1 (MQ=255)
                                                                                                    gTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGtgtg                                                <  4:859603/149‑1 (MQ=255)
                                                                                                      ccTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTgtggt                                              <  2:782551/149‑1 (MQ=255)
                                                                                                      ccTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTgtggt                                              <  1:419940/149‑1 (MQ=255)
                                                                                                       cTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTgtggtg                                             >  1:196724/1‑149 (MQ=255)
                                                                                                          ccGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGccc                                          >  1:241103/1‑149 (MQ=255)
                                                                                                               gCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTaa                                     >  4:198212/1‑149 (MQ=255)
                                                                                                                 cAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAaca                                   >  4:896718/1‑149 (MQ=255)
                                                                                                                 cAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAaca                                   >  2:678835/1‑149 (MQ=255)
                                                                                                                  aaCGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTgtggt                                              <  2:679133/137‑1 (MQ=255)
                                                                                                                  aaCGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTgtggt                                              >  1:679133/1‑137 (MQ=255)
                                                                                                                  aaCGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAacac                                  >  4:587682/1‑149 (MQ=255)
                                                                                                                   aCGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAAcaca                                 <  4:832467/149‑1 (MQ=255)
                                                                                                                     gCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATTACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATc                               >  2:758824/1‑149 (MQ=255)
                                                                                                                         gACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCt                           >  2:851316/1‑149 (MQ=255)
                                                                                                                               cGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGgt                                                   <  3:440430/119‑1 (MQ=255)
                                                                                                                               cGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGgt                                                   >  4:440430/1‑119 (MQ=255)
                                                                                                                                ggCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTg                    <  1:423554/149‑1 (MQ=255)
                                                                                                                                  cATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGc                  <  1:151026/149‑1 (MQ=255)
                                                                                                                                   aTCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCg                 >  1:869631/1‑149 (MQ=255)
                                                                                                                                          cACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGCTTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAg          <  3:644346/149‑1 (MQ=255)
                                                                                                                                          cACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGGGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAg          >  1:265413/1‑149 (MQ=255)
                                                                                                                                           ttaggCCTTAAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTcc            >  2:381744/4‑146 (MQ=255)
                                                                                                                                           ttaggCCTTAAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTcc            <  1:381744/143‑1 (MQ=255)
                                                                                                                                           ttaggCCTTAAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGa         >  4:605364/4‑149 (MQ=255)
                                                                                                                                            cTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGAc        >  3:3308/1‑149 (MQ=255)
                                                                                                                                             tGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGAcc       >  2:575830/1‑149 (MQ=255)
                                                                                                                                             tGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGAcc       >  2:721994/1‑149 (MQ=255)
                                                                                                                                              ggCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCt      <  1:822808/149‑1 (MQ=255)
                                                                                                                                              ggCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCt      <  4:800480/149‑1 (MQ=255)
                                                                                                                                                  tttAAGGCCTAACCCAACGTTATCAATCAGCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGAcc       >  1:743623/1‑144 (MQ=255)
                                                                                                                                                  tttAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCca  <  4:223536/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AACAATAATCCCGGTCGATGAATCAACGCTCTTGCCAGCGCCACTCGCTGTTTCTGCCCGCCAGAAAGTGCGGCAGGCCATTCCCCTGCGCGATTCTCCAGTCCTACCGCAGCCAACGCTCGACGTGCGGCATCGCGCCACTGGCCTTTAAGGCCTAACCCAACGTTATCAATCACCGATTTCCACGGTAGCAGACGCGCATCCTGAAACATCATTCGCGTGTCTTCCTGAATTTCAGCCAGCGGTGTGGTGCCCGCTAACACATCGCCTGCGGTTGGCGTTTCCAGACCTGCCA  >  NC_000913/993581‑993875

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