Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,129,903 T→A 7.3% N326Y (AAT→TAT)  wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,129,9030TA 7.3% 296.2 / 9.8 96N326Y (AAT→TAT) wcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (43/46);  new base A (3/4);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACG  >  NC_000913/2129755‑2130052
                                                                                                                                                    |                                                                                                                                                     
gATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTg                                                                                                                                                      <  1:1619571/150‑1 (MQ=255)
  tACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTc                                                                                                                                                    <  2:496641/150‑1 (MQ=255)
  tACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACAAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTAACACAGTGATAAAAAAACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTc                                                                                                                                                    >  2:614094/1‑150 (MQ=255)
    ccGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTa                                                                                                                                                  >  2:140183/1‑150 (MQ=255)
       gAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAAGACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATgg                                                                                                                                               <  2:782414/150‑1 (MQ=255)
       gAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATgg                                                                                                                                               <  1:88990/150‑1 (MQ=255)
          tAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttt                                                                                                                                            >  1:1595000/1‑150 (MQ=255)
           aaaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtttt                                                                                                                                           >  2:1334177/1‑150 (MQ=255)
           aaaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtttt                                                                                                                                           >  1:627443/1‑150 (MQ=255)
           aaaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttg                                                                                                                                           <  1:942159/150‑2 (MQ=255)
            aaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttttt                                                                                                                                          >  1:916552/1‑150 (MQ=255)
             aaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgta                                                                                                                                         <  1:1168049/150‑4 (MQ=255)
              aaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTcc                                                                                                                                        >  1:300998/1‑150 (MQ=255)
              aaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgtat                                                                                                                                        <  2:994049/150‑5 (MQ=255)
              aaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGAGAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTcc                                                                                                                                        >  1:306880/1‑150 (MQ=255)
                catatAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgtatct                                                                                                                                      <  1:854252/149‑7 (MQ=255)
                 atatAATTGCACGCGAAAATTGCGGTTATCAACACAAGCGTTGATGATCATTTTTATCACTTACCCCATGTACCACAGTGATAAAAACAACGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCTTTGGCTATGGTTTTTCCGAc                                                                                                                                     <  2:1196083/150‑1 (MQ=255)
                          cACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                                                                                                            >  2:319763/1‑150 (MQ=255)
                            cgcgAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGt                                                                                                                          >  2:726693/1‑150 (MQ=255)
                               gAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaa                                                                                                                       <  2:1007862/150‑1 (MQ=255)
                               gAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaa                                                                                                                       >  1:383897/1‑150 (MQ=255)
                                aaaaTTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCCCTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaaa                                                                                                                      <  2:594451/150‑1 (MQ=255)
                                 aaaTTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAat                                                                                                                     <  2:1595000/150‑1 (MQ=255)
                                 aaaTTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAat                                                                                                                     <  1:1073562/150‑1 (MQ=255)
                                     tGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTc                                                                                                                 <  1:1310242/150‑1 (MQ=255)
                                      gCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                                                                                                                <  2:310726/150‑1 (MQ=255)
                                        ggTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGa                                                                                                              >  2:201001/1‑150 (MQ=255)
                                        ggTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGa                                                                                                              <  2:913703/150‑1 (MQ=255)
                                         gTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGaa                                                                                                             <  1:726693/150‑1 (MQ=255)
                                          ttATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAat                                                                                                            <  1:1334177/150‑1 (MQ=255)
                                             tCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAatgat                                                                                                         >  1:794122/1‑150 (MQ=255)
                                               aCCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                                                                                                       >  2:959423/1‑150 (MQ=255)
                                                  aaaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAc                                                                                                    >  1:135407/1‑150 (MQ=255)
                                                  aaaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAc                                                                                                    <  1:1598965/150‑1 (MQ=255)
                                                  aaaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAc                                                                                                    <  2:177875/150‑1 (MQ=255)
                                                     aGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGTAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACata                                                                                                 <  1:1542600/150‑1 (MQ=255)
                                                          tGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCAGAATGATGCGACATATTCAt                                                                                            >  1:953373/1‑150 (MQ=255)
                                                              catcatTTTTATCACTTTCCCCATTTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTTTATATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAg                                                                                        >  2:150654/1‑150 (MQ=255)
                                                              catcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACAt                                                                                                                                   >  2:948667/1‑107 (MQ=255)
                                                              catcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACAt                                                                                                                                   <  1:948667/107‑1 (MQ=255)
                                                              catcatTTTTATCACTTTCCCCATGTACAACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAAATACAAACCATTGTATATGGTTTTTCCGAAATCCGCACCGTTAAATATTCCAAATGATGCGACATATTCATaaat                                                                                        >  2:144925/1‑149 (MQ=255)
                                                               atcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt                                                                                       >  1:814987/1‑150 (MQ=255)
                                                                       tATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa                                                                               <  2:886749/150‑1 (MQ=255)
                                                                       tATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa                                                                               >  2:1051845/1‑150 (MQ=255)
                                                                         tCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATc                                                                             <  2:916552/150‑1 (MQ=255)
                                                                              ttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACt                                                                        >  1:672853/1‑150 (MQ=255)
                                                                                  ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAc                                                                                                    <  2:517413/118‑1 (MQ=255)
                                                                                  ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAc                                                                                                    >  1:517413/1‑118 (MQ=255)
                                                                                   cATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCACACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACAcc                                                                   <  1:271783/150‑1 (MQ=255)
                                                                                    aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTACGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGCCTACACCg                                                                  >  2:1363922/1‑150 (MQ=255)
                                                                                      gTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTc                                                                >  2:97242/1‑150 (MQ=255)
                                                                                       tACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCa                                                               <  2:854537/150‑1 (MQ=255)
                                                                                       tACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCa                                                               >  2:438203/1‑150 (MQ=255)
                                                                                        aCCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCaa                                                              >  2:42524/1‑150 (MQ=255)
                                                                                             aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                         <  1:1051845/150‑1 (MQ=255)
                                                                                              gTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATggg                                                        <  1:140183/150‑1 (MQ=255)
                                                                                              gTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATggg                                                        <  2:1494188/150‑1 (MQ=255)
                                                                                                gATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGaa                                                                                                             <  1:1585319/95‑1 (MQ=255)
                                                                                                gATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGaa                                                                                                             >  2:1585319/1‑95 (MQ=255)
                                                                                                gATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTc                                                      >  2:769800/1‑150 (MQ=255)
                                                                                                gATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGAGAGCCAAAACAGACTACACCGTCAataaagcac                                                      >  1:288506/1‑144 (MQ=255)
                                                                                                  tAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaa                                                    >  2:1361654/1‑150 (MQ=255)
                                                                                                  tAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaa                                                    >  1:788568/1‑150 (MQ=255)
                                                                                                    aaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAg                                                  <  2:300998/150‑1 (MQ=255)
                                                                                                    aaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAg                                                  <  2:306880/150‑1 (MQ=255)
                                                                                                        cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                              <  1:319763/150‑1 (MQ=255)
                                                                                                        cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                              <  1:334647/150‑1 (MQ=255)
                                                                                                        cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                              <  1:896495/150‑1 (MQ=255)
                                                                                                         acCGCGAACCAGGAAAAATAAATAATCAGCAGATACAACCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGGGACATATTCATAAATTGAGCCAAATCTGACTACACCGTCAATATGGGTGAATGAATaa                                             <  2:1191635/150‑1 (MQ=255)
                                                                                                          ccGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAAc                                            >  2:1521056/1‑150 (MQ=255)
                                                                                                              gAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGAc                                        <  2:814987/150‑1 (MQ=255)
                                                                                                                 ccAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCAt                                     >  2:81169/1‑150 (MQ=255)
                                                                                                                  cAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATc                                    <  2:747360/150‑1 (MQ=255)
                                                                                                                   aGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCa                                   >  1:1624530/1‑150 (MQ=255)
                                                                                                                   aGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCa                                   >  2:455609/1‑150 (MQ=255)
                                                                                                                     gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAcc                                 <  2:1283522/150‑1 (MQ=255)
                                                                                                                     gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAcc                                 <  2:618048/150‑1 (MQ=255)
                                                                                                                         aataaataATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCg                             >  1:396511/1‑150 (MQ=255)
                                                                                                                                aatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGAc                      <  2:794122/150‑1 (MQ=255)
                                                                                                                                aatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGAc                      <  2:588825/150‑1 (MQ=255)
                                                                                                                                aatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGAc                      <  2:672853/150‑1 (MQ=255)
                                                                                                                                 atCAGCAGATACAAACCATTGTCTATGTTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACa                     <  2:1176227/150‑1 (MQ=255)
                                                                                                                                  tCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACaa                    >  1:934197/1‑150 (MQ=255)
                                                                                                                                        gATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGat              <  1:614094/150‑1 (MQ=255)
                                                                                                                                        cggaaaaaCCATAGACAATGGTTTTTCAGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGat              >  1:666754/6‑150 (MQ=255)
                                                                                                                                         aTACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGata             >  1:1386179/1‑150 (MQ=255)
                                                                                                                                          tACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataa            <  2:1210889/150‑1 (MQ=255)
                                                                                                                                          gaaaaaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataa            >  2:10198/4‑150 (MQ=255)
                                                                                                                                            cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataata          >  2:1196355/1‑150 (MQ=255)
                                                                                                                                            cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataata          >  1:604738/1‑150 (MQ=255)
                                                                                                                                            cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataata          <  1:97242/150‑1 (MQ=255)
                                                                                                                                            cAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAAAATACCGAATGAAGCCACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAAAAACCGACCATCACCAGCTGACCGACAATACGataata          >  1:1489686/1‑150 (MQ=255)
                                                                                                                                            aaaaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataata          >  1:766168/2‑150 (MQ=255)
                                                                                                                                             aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAatat         >  2:808997/1‑150 (MQ=255)
                                                                                                                                              aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATatg        <  1:769800/150‑1 (MQ=255)
                                                                                                                                               aCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATatga       <  1:42524/150‑1 (MQ=255)
                                                                                                                                                    tGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACg  <  1:42277/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACG  >  NC_000913/2129755‑2130052

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