Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,374,901 A→T 8.5% L11F (TTA→TTT dcuD → putative transporter DcuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,9010AT 8.5% 326.1 / 14.0 116L11F (TTA→TTTdcuDputative transporter DcuD
Reads supporting (aligned to +/- strand):  ref base A (43/63);  new base T (2/8);  total (45/71)
Fisher's exact test for biased strand distribution p-value = 3.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 128 total aligned reads displayed.

CCTACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATAT  >  NC_000913/3374752‑3375050
                                                                                                                                                     |                                                                                                                                                     
cctACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTatta                                                                                                                                                       >  1:43912/1‑150 (MQ=255)
  tACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaat                                                                                                                                                     <  1:788323/150‑1 (MQ=255)
  tACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaat                                                                                                                                                     <  2:1118060/150‑1 (MQ=255)
   gcAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaatt                                                                                                                                                    <  1:936488/149‑1 (MQ=255)
    cAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaatta                                                                                                                                                   <  2:505761/150‑1 (MQ=255)
     aaTTAAAAACACGAGATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTAc                                                                                                                                                  <  1:1180105/150‑1 (MQ=255)
     aaTTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTTTGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTTTTTATTAc                                                                                                                                                  >  2:1422586/1‑150 (MQ=255)
      aTTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACg                                                                                                                                                 >  2:558019/1‑150 (MQ=255)
      aTTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACg                                                                                                                                                 >  1:1320347/1‑150 (MQ=255)
       ttAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGa                                                                                                                                                >  1:11202/1‑150 (MQ=255)
       ttAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGa                                                                                                                                                <  1:890990/150‑1 (MQ=255)
        tAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGAt                                                                                                                                               <  2:1295534/150‑1 (MQ=255)
         aaaaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATg                                                                                                                                              <  2:323630/150‑1 (MQ=255)
         aaaaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATg                                                                                                                                              <  2:188359/150‑1 (MQ=255)
           aaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATggg                                                                                                                                            >  1:830486/1‑150 (MQ=255)
            aaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGc                                                                                                                                           <  2:123905/150‑1 (MQ=255)
            aaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGc                                                                                                                                           <  1:883780/150‑1 (MQ=255)
             acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                          <  1:560038/150‑5 (MQ=255)
             acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                          <  2:1152567/150‑5 (MQ=255)
             acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                          <  2:1478035/150‑5 (MQ=255)
             acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                          <  2:247539/150‑5 (MQ=255)
              cacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTa                                                                                                                                         >  1:55966/1‑150 (MQ=255)
              cacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacaga                                                                                                                                         <  1:1292812/150‑6 (MQ=255)
              cacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacaga                                                                                                                                         <  1:921441/150‑6 (MQ=255)
                cGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTAtt                                                                                                                                       <  2:6972/150‑1 (MQ=255)
                cGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTAtt                                                                                                                                       <  1:1200123/150‑1 (MQ=255)
                 gACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTAttt                                                                                                                                      <  1:857002/150‑1 (MQ=255)
                  aCATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTg                                                                                                                                     >  1:219429/1‑150 (MQ=255)
                  aCATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTg                                                                                                                                     <  2:231903/150‑1 (MQ=255)
                        gCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTg                                                                                                                               <  2:219429/150‑1 (MQ=255)
                          aCCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGaa                                                                                                                             >  2:830195/1‑150 (MQ=255)
                           ccATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGTCACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGaaa                                                                                                                            <  1:1171659/150‑1 (MQ=255)
                           ccATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGaaa                                                                                                                            <  1:1081383/150‑1 (MQ=255)
                            cATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGaaaa                                                                                                                           <  2:207663/150‑1 (MQ=255)
                               aaaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACt                                                                                                                        <  1:1097160/150‑1 (MQ=255)
                                aaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTa                                                                                                                       <  2:695250/150‑1 (MQ=255)
                                  tAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACa                                                                                                                     <  1:804154/150‑1 (MQ=255)
                                    gCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCATCAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACaaa                                                                                                                   >  1:19708/1‑150 (MQ=255)
                                    gCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACaaa                                                                                                                   >  1:1018978/1‑150 (MQ=255)
                                    gCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACaaa                                                                                                                   >  1:1311707/1‑150 (MQ=255)
                                    gCCTTCCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGCGATTGACCGATTTCTTTTGACACGGAGTTCAACAATTTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATATGATCCTGAAAAACAACaaa                                                                                                                   >  2:872258/1‑150 (MQ=255)
                                        tGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTc                                                                                                               >  2:1536148/1‑150 (MQ=255)
                                           aaaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAgg                                                                                                            >  1:422466/1‑150 (MQ=255)
                                           aaaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATATGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAgg                                                                                                            <  2:1517105/150‑1 (MQ=255)
                                            aaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATTGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggt                                                                                                           <  2:807338/150‑1 (MQ=255)
                                                  ataACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgct                                                                                                     >  1:1463875/1‑150 (MQ=255)
                                                   taACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctg                                                                                                    <  1:478635/150‑1 (MQ=255)
                                                    aaCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctgg                                                                                                   >  1:455426/1‑150 (MQ=255)
                                                    aaCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctgg                                                                                                   <  2:1311707/150‑1 (MQ=255)
                                                    aaCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctgg                                                                                                   <  1:1325635/150‑1 (MQ=255)
                                                    aaCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctgg                                                                                                   <  1:396303/150‑1 (MQ=255)
                                                     aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                  <  2:1326002/150‑1 (MQ=255)
                                                     aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                  <  2:55966/150‑1 (MQ=255)
                                                     aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                  <  1:323156/150‑1 (MQ=255)
                                                      cATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggct                                                                                                 <  1:1456757/150‑1 (MQ=255)
                                                            tgttTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCa                                                                                           <  2:1018978/150‑1 (MQ=255)
                                                            tgttTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCa                                                                                           <  2:1047001/150‑1 (MQ=255)
                                                             gttTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAg                                                                                          <  1:1145962/150‑1 (MQ=255)
                                                              ttttCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAgg                                                                                         >  2:933122/1‑150 (MQ=255)
                                                                  cAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATc                                                                                     <  2:11202/150‑1 (MQ=255)
                                                                  cAATATGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTAAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATc                                                                                     >  2:740181/1‑150 (MQ=255)
                                                                      cTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTcc                                                                                 <  1:558019/150‑1 (MQ=255)
                                                                       tGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCt                                                                                >  2:101279/1‑150 (MQ=255)
                                                                        gCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtg                                                                               <  2:375081/150‑1 (MQ=255)
                                                                          cGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtgat                                                                             >  2:1443399/1‑150 (MQ=255)
                                                                           gTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtgatg                                                                            <  1:969876/150‑1 (MQ=255)
                                                                            tttATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtgatga                                                                           <  2:1277430/150‑1 (MQ=255)
                                                                                 gggATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGc                                                                      <  2:207708/150‑1 (MQ=255)
                                                                                 gggATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGc                                                                      <  1:841457/150‑1 (MQ=255)
                                                                                  ggATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCg                                                                     >  1:1009226/1‑150 (MQ=255)
                                                                                  ggATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCg                                                                     <  1:663684/150‑1 (MQ=255)
                                                                                       gACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTc                                                                >  1:1425137/1‑150 (MQ=255)
                                                                                              ttCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAg                                                         >  1:1204613/1‑150 (MQ=255)
                                                                                              ttCTTTTGACACGGAGGTGAAGAATGTTCGGCATAATTTTATCTGTCATGGTATTAATTAGGATGGGCTATTTGATCCTGAAAAACTTCAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATGTTCCTGATGATGTGCGGTGTCTGGTTAg                                                         <  1:769855/150‑1 (MQ=255)
                                                                                                cTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAggg                                                       >  2:120282/1‑150 (MQ=255)
                                                                                                   ttGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTc                                                    <  2:1320347/150‑1 (MQ=255)
                                                                                                       cacGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggtg                                                >  2:394851/1‑150 (MQ=255)
                                                                                                          ggAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAATAACTACAAACCTCAGGTGGTGCTGGCTGCCTCAGGTATGTTCCTGATGATGTGCGGTGTCTGGTTAGTGTTCGgtggtg                                             <  2:661893/150‑1 (MQ=255)
                                                                                                           gAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGgtggtg                                             >  2:329462/1‑149 (MQ=255)
                                                                                                           gAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGgtggtg                                             <  1:329462/149‑1 (MQ=255)
                                                                                                                caacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCg                                       >  2:60645/1‑150 (MQ=255)
                                                                                                                caacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCg                                       >  2:477515/1‑150 (MQ=255)
                                                                                                                  acaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGAt                                     >  1:308055/1‑150 (MQ=255)
                                                                                                                     aTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATccc                                  >  2:70412/1‑150 (MQ=255)
                                                                                                                       gTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCccac                                <  2:1463875/150‑1 (MQ=255)
                                                                                                                       gTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCccac                                >  1:1148928/1‑150 (MQ=255)
                                                                                                                                aaTTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAagcag                       >  2:111149/1‑150 (MQ=255)
                                                                                                                                   tataTCTGTCATCGTATAAATTACGATGGGCTATTTGAACCCGAAAAACTACAAACCTCAGCTGGCGCTGGCTGCCGCAGGGATCTTCCGGATGATGTGCGGTGTCTGGGTAGGGTTCGGGGGTGTACTCGAGCCCACCAAAAGCAGCgg                    <  2:903079/150‑1 (MQ=255)
                                                                                                                                    atatCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGc                   <  1:969656/150‑1 (MQ=255)
                                                                                                                                         tGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACtt              >  1:275093/1‑150 (MQ=255)
                                                                                                                                         gcccATCGTAATTAATACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACtt              >  2:1002641/4‑150 (MQ=255)
                                                                                                                                           tCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACttta            <  2:275093/150‑3 (MQ=255)
                                                                                                                                              tCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATAATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcg         >  2:771086/1‑150 (MQ=255)
                                                                                                                                                gTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGTTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc       >  1:924934/1‑150 (MQ=255)
                                                                                                                                                 tATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtcg      >  2:487544/1‑150 (MQ=255)
                                                                                                                                                   ttaattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGat    <  1:1536148/150‑1 (MQ=255)
                                                                                                                                                    taattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGata   <  2:1506332/150‑1 (MQ=255)
                                                                                                                                                    taattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGata   <  1:740181/150‑1 (MQ=255)
                                                                                                                                                    taattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGata   <  1:864846/150‑1 (MQ=255)
                                                                                                                                                     aattaCGATGGGCTATTTGATCCTGAAAAACTACAAAGCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGatat  <  1:612228/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCTACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATAT  >  NC_000913/3374752‑3375050

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