Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,129,907 T→A 7.1% I324I (ATA→ATT wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,129,9070TA 7.1% 289.7 / 8.2 97I324I (ATA→ATTwcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (42/48);  new base A (3/4);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCC  >  NC_000913/2129759‑2130056
                                                                                                                                                    |                                                                                                                                                     
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                                                                                                                                               >  1:627443/1‑150 (MQ=255)
       aaaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtttt                                                                                                                                               >  2:1334177/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                                                                                                                         <  1:1073562/150‑1 (MQ=255)
                             aaaTTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAat                                                                                                                         <  2:1595000/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:1598965/150‑1 (MQ=255)
                                              aaaaGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAc                                                                                                        >  1:135407/1‑150 (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                                                                                                                                       <  1:948667/107‑1 (MQ=255)
                                                          catcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACAt                                                                                                                                       >  2:948667/1‑107 (MQ=255)
                                                          catcatTTTTATCACTTTCCCCATGTACAACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAAATACAAACCATTGTATATGGTTTTTCCGAAATCCGCACCGTTAAATATTCCAAATGATGCGACATATTCATaaat                                                                                            >  2:144925/1‑149 (MQ=255)
                                                           atcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt                                                                                           >  1:814987/1‑150 (MQ=255)
                                                                   tATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa                                                                                   >  2:1051845/1‑150 (MQ=255)
                                                                   tATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa                                                                                   <  2:886749/150‑1 (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:438203/1‑150 (MQ=255)
                                                                                   tACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCa                                                                   <  2:854537/150‑1 (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                                                        >  1:788568/1‑150 (MQ=255)
                                                                                              tAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaa                                                        >  2:1361654/1‑150 (MQ=255)
                                                                                                aaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAg                                                      <  2:306880/150‑1 (MQ=255)
                                                                                                aaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAg                                                      <  2:300998/150‑1 (MQ=255)
                                                                                                    cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                                  <  1:334647/150‑1 (MQ=255)
                                                                                                    cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                                  <  1:319763/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)
                                                                                                                                                 gTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGa     <  1:1309070/150‑1 (MQ=255)
                                                                                                                                                   cTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATc   <  2:1138956/150‑1 (MQ=255)
                                                                                                                                                   cTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATc   <  1:438203/150‑1 (MQ=255)
                                                                                                                                                    tATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAATTCTTACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATcc  <  2:722176/150‑1 (MQ=255)
                                                                                                                                                    tATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATcc  <  1:1397291/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCC  >  NC_000913/2129759‑2130056

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: