Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,900,327 A→C 7.3% intergenic (‑54/‑265) ygcW ← / → yqcE putative SDR family oxidoreductase/putative MFS transporter, inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,900,3270AC 7.3% 990.7 / 52.8 286intergenic (‑54/‑265)ygcW/yqcEputative SDR family oxidoreductase/putative MFS transporter, inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (110/155);  new base C (6/15);  total (116/170)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 299 total aligned reads displayed.

CTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAACAC  >  NC_000913/2900178‑2900476
                                                                                                                                                      |                                                                                                                                                     
 tGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgag                                                                                                                                                      >  1:3754732/1‑150 (MQ=255)
  ggCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgaga                                                                                                                                                     <  1:625463/150‑1 (MQ=255)
  ggCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCgaga                                                                                                                                                     <  2:4054489/150‑1 (MQ=255)
   gCCTAAACCGCTCTTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGagag                                                                                                                                                    <  2:4587477/150‑1 (MQ=255)
   gCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGagag                                                                                                                                                    >  2:1071286/1‑150 (MQ=255)
      tAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATGTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGa                                                                                                                                                 <  2:5178380/150‑1 (MQ=255)
       aaaCCGCTATTCCCACCGGTAACAATTGCGTTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                                <  1:756166/150‑1 (MQ=255)
       aaaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGaa                                                                                                                                                <  2:1433572/150‑1 (MQ=255)
       aaaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGaa                                                                                                                                                <  1:5753880/150‑1 (MQ=255)
        aaCCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAg                                                                                                                                               <  2:3783890/150‑1 (MQ=255)
           cGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGTTa                                                                                                                                            <  2:4972329/150‑1 (MQ=255)
           cGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTa                                                                                                                                            >  1:5149193/1‑150 (MQ=255)
            gCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                           >  2:381425/1‑150 (MQ=255)
            gCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                           <  2:2359362/150‑1 (MQ=255)
              tagtCCCACCGGTAACAATTGCGGTTTTACCTCTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGTTACCt                                                                                                                                         <  1:4671490/147‑1 (MQ=39)
                 tCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtc                                                                                                                                      <  2:3979755/150‑1 (MQ=255)
                  cccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCTCTCGCGGGGGAAGTTACCTTtct                                                                                                                                     <  1:4957104/150‑1 (MQ=255)
                   ccACCGGTAACTATTGCGTTTTTACCTTTCATGTATAAAAAATCCATTGATAACGCATTGAGAGATTCGATTGACATAGTTTGGTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtctc                                                                                                                                    <  1:4551964/150‑1 (MQ=255)
                   ccACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTtctc                                                                                                                                    >  1:3083718/1‑150 (MQ=255)
                            aCAAGTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGACCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATctct                                                                                                                           <  1:3974748/150‑1 (MQ=255)
                              aaTTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTgg                                                                                                                         <  2:587268/150‑1 (MQ=255)
                              aaTTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTgg                                                                                                                         <  2:3662960/150‑1 (MQ=255)
                               aTTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGt                                                                                                                        <  2:1695413/150‑1 (MQ=255)
                                 tGCGGTTTTACCTTTCAGGTATAAAAATTCAATTGTGAACGATTTTTTATTTTCGATTTACATTTTTACCTCCATTTTCCGTAAGGAAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTaa                                                                                                                      <  2:711635/150‑1 (MQ=255)
                                    ggTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTg                                                                                                                   >  1:5700361/1‑150 (MQ=255)
                                     gTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTgg                                                                                                                  >  1:1187881/1‑150 (MQ=255)
                                     gTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTgg                                                                                                                  <  2:1816469/150‑1 (MQ=255)
                                      ttttACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGga                                                                                                                 >  2:5307033/1‑150 (MQ=255)
                                        ttACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGgaga                                                                                                               <  1:2657021/150‑1 (MQ=255)
                                         tACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGAt                                                                                                              >  2:4786191/1‑150 (MQ=255)
                                              ttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAttt                                                                                                         >  1:1623159/1‑150 (MQ=255)
                                              ttCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTGTCTCATTATCTCTGGTAACTGGAGATAAttt                                                                                                         >  1:3029264/1‑150 (MQ=255)
                                                cAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAAttttt                                                                                                       <  1:5373292/150‑1 (MQ=255)
                                                  gggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAg                                                                                                     <  1:3846047/150‑1 (MQ=255)
                                                   ggAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGc                                                                                                    >  1:3328476/1‑150 (MQ=255)
                                                       aaaaaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAGAAAGGTAACTTCCCACGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAg                                                                                               >  2:3442514/1‑150 (MQ=255)
                                                         aaaaTCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAgg                                                                                              >  1:4840806/1‑150 (MQ=255)
                                                               cATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCCTaa                                                                                        >  2:4188566/1‑150 (MQ=255)
                                                                 ttGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACt                                                                                      >  1:5003142/1‑150 (MQ=255)
                                                                  tGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTg                                                                                     <  2:2447739/150‑1 (MQ=255)
                                                                  tGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTg                                                                                     <  2:3887182/150‑1 (MQ=255)
                                                                  tGAGAACGCATTGAGAGATTCGATTGACAGAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTg                                                                                     <  2:549953/150‑1 (MQ=255)
                                                                    agaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCt                                                                                   >  2:301541/1‑150 (MQ=255)
                                                                     gaACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTa                                                                                  >  2:1305640/1‑150 (MQ=255)
                                                                       aCGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTTACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAt                                                                                <  1:3310525/150‑1 (MQ=255)
                                                                         gCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAAtgt                                                                              <  1:1407002/150‑1 (MQ=255)
                                                                            ttGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGAt                                                                           >  1:2784266/1‑150 (MQ=255)
                                                                             tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                          >  1:1507219/1‑150 (MQ=255)
                                                                             tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                          <  1:1071286/150‑1 (MQ=255)
                                                                             tGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATc                                                                          <  1:4450936/150‑1 (MQ=255)
                                                                                gagaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTa                                                                       >  2:246544/1‑150 (MQ=255)
                                                                                 agaTTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTaa                                                                      <  1:944707/150‑1 (MQ=255)
                                                                                    ttCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATca                                                                   <  2:5616410/150‑1 (MQ=255)
                                                                                       gATTGACATAGTTAGCTACAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTCACTGGAGATAATTTTTAGAATAAGGGCATAACTGCTAATGTGATCGTAATCACAg                                                                >  2:1713872/1‑150 (MQ=255)
                                                                                        aTTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgt                                                               <  2:491358/150‑1 (MQ=255)
                                                                                          tGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAgtgt                                                             <  1:5226290/150‑1 (MQ=255)
                                                                                            aCATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGa                                                           >  1:591825/1‑150 (MQ=255)
                                                                                              aTAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGata                                                         >  2:1055174/1‑150 (MQ=255)
                                                                                                aGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATAtt                                                       <  2:2916862/150‑1 (MQ=255)
                                                                                                 gTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTc                                                      <  1:5284241/150‑1 (MQ=255)
                                                                                                   tAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGt                                                    <  1:2532714/150‑1 (MQ=255)
                                                                                                   tAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGt                                                    <  2:1789759/150‑1 (MQ=255)
                                                                                                     tcccccaTTTACGTAAAGCCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTa                                                   <  1:4196032/144‑1 (MQ=255)
                                                                                                      cTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTAtt                                                 >  1:5293136/1‑150 (MQ=255)
                                                                                                      cTCCAATTTCCGTAAGGCAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTg                                                >  2:3199427/1‑150 (MQ=255)
                                                                                                       tCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACGGCTAATGTGATCGTAATCACAGTGTGATATTCGTCTTg                                                >  2:4731706/1‑150 (MQ=255)
                                                                                                       tCCAATTTCCGTAAGGCAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGt                                               <  1:5432304/150‑1 (MQ=255)
                                                                                                         cAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTc                                              <  2:3332025/150‑1 (MQ=255)
                                                                                                         cAATGTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTc                                              <  2:5355763/150‑1 (MQ=255)
                                                                                                           aTTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCGGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCaa                                            >  2:450686/1‑150 (MQ=255)
                                                                                                           aTTTCCCTAAGTCAAAAAAAAACAAACCTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCaa                                            <  1:3536663/150‑1 (MQ=255)
                                                                                                            tttCCGTAAGCCAAAAAAAAAGAAAGTTAACTTCCCCCTCGAGAGTGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATATCTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCaaa                                           <  1:3172910/150‑1 (MQ=255)
                                                                                                              tCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGGGATATTCGTATTGTCAaata                                         >  2:961400/1‑150 (MQ=255)
                                                                                                                cGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAaatgat                                       >  1:3449137/1‑150 (MQ=255)
                                                                                                                 gTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtt                                      <  1:129352/150‑1 (MQ=255)
                                                                                                                 gTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtt                                      <  1:2435520/150‑1 (MQ=255)
                                                                                                                 gTAAGGCAAAAAAAAAGAAAGGGAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAAtt                                      <  2:1356137/150‑1 (MQ=255)
                                                                                                                     aGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                           >  1:5133865/1‑46 (MQ=255)
                                                                                                                     aGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTAc                                                                                                                                           <  2:5133865/46‑1 (MQ=255)
                                                                                                                      ggCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaa                                  <  1:1941608/150‑1 (MQ=255)
                                                                                                                        cAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaa                                >  1:3855624/1‑150 (MQ=255)
                                                                                                                        cAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaa                                <  2:2779076/150‑1 (MQ=255)
                                                                                                                        cAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTAATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGGAATCACAGTGTGATATTCGTATTGTCAAATAATGTTaaaa                                >  2:3821214/1‑150 (MQ=255)
                                                                                                                         aaaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaaa                               <  2:572273/150‑1 (MQ=255)
                                                                                                                         aaaaaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTaaaaa                               >  1:5874804/1‑150 (MQ=255)
                                                                                                                            aaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCa                            >  1:4791892/1‑150 (MQ=255)
                                                                                                                            aaaaaaGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCa                            <  1:2947142/150‑1 (MQ=255)
                                                                                                                                    aagGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTg                    <  2:852069/150‑1 (MQ=255)
                                                                                                                                      ggTAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                  >  2:4396509/1‑150 (MQ=255)
                                                                                                                                      ggTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGAt                  <  2:3931748/150‑1 (MQ=255)
                                                                                                                                        tAACTTCCCTCTCGCGGGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTa                >  2:2036574/1‑150 (MQ=255)
                                                                                                                                        tAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTa                >  1:5645751/1‑150 (MQ=255)
                                                                                                                                         aaCTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTAc               >  1:2843028/1‑150 (MQ=255)
                                                                                                                                          aCTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTAct              <  1:4351570/150‑1 (MQ=255)
                                                                                                                                           cTTCCCTCTCGCGGGGGAAGTTACCTTTCTTATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActc             >  2:777219/1‑150 (MQ=255)
                                                                                                                                           cTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTActc             <  1:3510925/150‑1 (MQ=255)
                                                                                                                                               ccccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcac         >  2:5246575/1‑150 (MQ=255)
                                                                                                                                                cccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACGGGAGATAATTTTTAGCATAAGGGCCTAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAATCAGAAGCGTGATTACTCTcaca        >  1:2081624/1‑150 (MQ=255)
                                                                                                                                                 ccGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTcacac       >  1:4819238/1‑150 (MQ=255)
                                                                                                                                                    cgAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACacaac    >  2:4522564/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGGCCTAAACCGCTATTCCCACCGGTAACAATTGCGGTTTTACCTTTCAGGGAGAAAAAATCCATTGAGAACGCATTGAGAGATTCGATTGACATAGTTAGCTCCAATTTCCGTAAGGCAAAAAAAAAGAAAGGTAACTTCCCCCGCGAGAGGGAAGTTACCTTTCTCATTATCTCTGGTAACTGGAGATAATTTTTAGCATAAGGGCATAACTGCTAATGTGATCGTAATCACAGTGTGATATTCGTATTGTCAAATAATTTTAAAAAGCAGAAGCGTGATTACTCTCACACAACAC  >  NC_000913/2900178‑2900476

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: