Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,129,974 A→G 5.8% F302S (TTT→TCT)  wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,129,9740AG 5.8% 340.7 / 8.0 104F302S (TTT→TCT) wcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base A (32/66);  new base G (3/3);  total (35/69)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 104 total aligned reads displayed.

ATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACACCAACCAGG  >  NC_000913/2129827‑2130120
                                                                                                                                                   |                                                                                                                                                  
aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAAt                                                                                                                                                  <  1:1149445/150‑1 (MQ=255)
aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAAt                                                                                                                                                  <  1:1151389/150‑1 (MQ=255)
aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAAt                                                                                                                                                  <  1:1845641/150‑1 (MQ=255)
aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCAATGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAAt                                                                                                                                                  <  1:809845/150‑1 (MQ=255)
 tCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCt                                                                                                                                                <  1:422258/150‑1 (MQ=255)
  cACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCt                                                                                                                                                <  1:708724/150‑1 (MQ=255)
    cTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                                                                                              <  1:1083830/150‑1 (MQ=255)
    cTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                                                                                              <  1:858715/150‑1 (MQ=255)
    cTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                                                                                              <  1:65599/150‑1 (MQ=255)
    cTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                                                                                              <  1:581335/150‑1 (MQ=255)
    cTTTCCCCATGGACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                                                                                              <  1:65594/150‑1 (MQ=255)
     tttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGAc                                                                                                                                             >  1:2332495/1‑150 (MQ=255)
     tttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGAc                                                                                                                                             >  1:2205926/1‑150 (MQ=255)
     tttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGAc                                                                                                                                             >  1:1802198/1‑150 (MQ=255)
          ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTacac                                                                                                                                        >  1:1119186/1‑150 (MQ=255)
             tGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGt                                                                                                                                     >  2:1049514/1‑150 (MQ=255)
                     aGTGATAAAAACACCGCTAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                                                                                             <  2:2085942/150‑1 (MQ=255)
                     aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                                                                                             <  1:1334373/150‑1 (MQ=255)
                     aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                                                                                             <  1:345826/150‑1 (MQ=255)
                     aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                                                                                             <  1:345854/150‑1 (MQ=255)
                     aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                                                                                             <  1:480973/150‑1 (MQ=255)
                     aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                                                                                             <  2:1356110/150‑1 (MQ=255)
                     aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                                                                                             <  2:1051717/150‑1 (MQ=255)
                       tGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGt                                                                                                                           >  2:2378242/1‑150 (MQ=255)
                           aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaaa                                                                                                                       >  2:2004104/1‑150 (MQ=255)
                           aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaaa                                                                                                                       <  1:558846/150‑1 (MQ=255)
                           aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaaa                                                                                                                       <  1:160409/150‑1 (MQ=255)
                           aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaaa                                                                                                                       <  1:2286773/150‑1 (MQ=255)
                           aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCaaa                                                                                                                       <  1:2286923/150‑1 (MQ=255)
                                 acCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATaa                                                                                                                 <  1:1336415/150‑1 (MQ=255)
                                  ccGCGAACCAGGAAAAATAAATTATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAAc                                                                                                                <  2:1972164/150‑1 (MQ=255)
                                   cgcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAAcc                                                                                                               <  2:1682465/150‑1 (MQ=255)
                                     cgAACCAGTAAAAATAAATAATCACCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATACATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGa                                                                                                             <  1:1520434/150‑1 (MQ=255)
                                             gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAcc                                                                                                     >  2:1725019/1‑150 (MQ=255)
                                             gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAcc                                                                                                     >  1:1725517/1‑150 (MQ=255)
                                             gAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAcc                                                                                                     >  2:434335/1‑150 (MQ=255)
                                                    aaataaTCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATTTGGGTCAAAGAATAACCGACCATCACCAGCGGAc                                                                                              <  1:593146/150‑1 (MQ=255)
                                                    aaataaTCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCCAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGAc                                                                                              <  2:1343857/150‑1 (MQ=255)
                                                     aataatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGAcc                                                                                             >  2:739893/1‑150 (MQ=255)
                                                       taatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGa                                                                                           <  2:546601/150‑1 (MQ=255)
                                                       taatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGa                                                                                           <  2:485903/150‑1 (MQ=255)
                                                        aatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGAc                                                                                          <  2:1119186/150‑1 (MQ=255)
                                                         atCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAGATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACa                                                                                         >  1:25502/1‑150 (MQ=255)
                                                         atCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAGATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACa                                                                                         >  1:1939718/1‑150 (MQ=255)
                                                            agcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATa                                                                                      <  2:1239237/150‑1 (MQ=255)
                                                            agcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATa                                                                                      <  2:1239254/150‑1 (MQ=255)
                                                            agcagATACAAACCATTGTCTATGGTTTTTCCGAAATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACaaga                                                                                      <  2:2275753/150‑3 (MQ=255)
                                                                     aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAatat                                                                             >  1:2190943/1‑150 (MQ=255)
                                                                     aaaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAatat                                                                             >  2:1384194/1‑150 (MQ=255)
                                                                      aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATatg                                                                            >  1:501731/1‑150 (MQ=255)
                                                                      aaCCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATatg                                                                            >  1:120333/1‑150 (MQ=255)
                                                                           ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATTACCAGCGGACCGACAATACGGTAATATGATGAc                                                                       <  2:313350/150‑1 (MQ=255)
                                                                           ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGAc                                                                       <  2:2120774/150‑1 (MQ=255)
                                                                           ttGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGAc                                                                       <  2:2361852/150‑1 (MQ=255)
                                                                              tCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGAt                                                                    <  2:2099163/150‑1 (MQ=255)
                                                                              tCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGAt                                                                    <  2:587331/150‑1 (MQ=255)
                                                                              tCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGAt                                                                    <  2:2180431/150‑1 (MQ=255)
                                                                               cTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATc                                                                   <  1:433797/150‑1 (MQ=255)
                                                                               cTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATc                                                                   <  2:1991629/150‑1 (MQ=255)
                                                                               cTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATc                                                                   <  2:2227293/150‑1 (MQ=255)
                                                                               cTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATc                                                                   <  1:1645846/150‑1 (MQ=255)
                                                                                tATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATcc                                                                  >  2:1634139/1‑150 (MQ=255)
                                                                                 aTGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCt                                                                 >  2:856731/1‑150 (MQ=255)
                                                                                   ggTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTc                                                               >  1:1238697/1‑150 (MQ=255)
                                                                                    gTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGTGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGAAAATACGATAATATGATGACGATCCTTcc                                                              <  2:1652589/150‑1 (MQ=255)
                                                                                    gTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTcc                                                              >  2:2116004/1‑150 (MQ=255)
                                                                                    gTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTcc                                                              <  2:95659/150‑1 (MQ=255)
                                                                                    gTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTcc                                                              <  2:95641/150‑1 (MQ=255)
                                                                                               aTCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATc                                                   <  1:1202582/150‑1 (MQ=255)
                                                                                               aTCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATc                                                   <  1:1407981/150‑1 (MQ=255)
                                                                                                   gCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACTATACGATAATATGATGACGATCCTTCCGTACCTAAATCGccc                                               <  1:174307/150‑1 (MQ=255)
                                                                                                    cACCGTTAAAAATTCCGAATGATTCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGTGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCa                                              >  1:1836336/1‑150 (MQ=255)
                                                                                                           aaaTATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGt                                       >  1:1901761/1‑150 (MQ=255)
                                                                                                           aaaTATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGt                                       <  1:763995/150‑1 (MQ=255)
                                                                                                           aaaTATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGACAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGt                                       <  1:1154035/150‑1 (MQ=255)
                                                                                                                ttCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGaaa                                  >  2:1158254/1‑150 (MQ=255)
                                                                                                                ttCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGaaa                                  >  2:2353029/1‑150 (MQ=255)
                                                                                                                    gAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAAtata                              <  1:2004104/150‑1 (MQ=255)
                                                                                                                     aaTGATGCGACATATTCATAAAGTGAGCCAAATCTGCCTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATaa                             <  1:1871658/150‑1 (MQ=255)
                                                                                                                     aaTGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATaa                             <  1:604217/150‑1 (MQ=255)
                                                                                                                      atgatgCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAg                            <  1:1443801/150‑1 (MQ=255)
                                                                                                                        gatgCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGa                          >  2:1107304/1‑150 (MQ=255)
                                                                                                                          tgCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGaaa                        >  1:296899/1‑150 (MQ=255)
                                                                                                                             gACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAAcgc                     <  2:1238697/150‑1 (MQ=255)
                                                                                                                                   ttCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTAc               <  2:1036173/150‑1 (MQ=255)
                                                                                                                                   ttCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTAc               <  2:1352614/150‑1 (MQ=255)
                                                                                                                                     cATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACca             <  1:2378242/150‑1 (MQ=255)
                                                                                                                                     cATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACca             <  1:1984761/150‑1 (MQ=255)
                                                                                                                                      aTAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACcac            >  2:531500/1‑150 (MQ=255)
                                                                                                                                       tAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACcaca           <  1:2380980/150‑1 (MQ=255)
                                                                                                                                          aGTGAGCCAGATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACacca        <  2:25502/150‑1 (MQ=255)
                                                                                                                                          aGTGAGCCAGATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACacca        <  2:1939718/150‑1 (MQ=255)
                                                                                                                                          aGTGAGCCAGATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACacca        <  2:1939644/150‑1 (MQ=255)
                                                                                                                                          aGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACacca        >  1:2094799/1‑150 (MQ=255)
                                                                                                                                          aGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACacca        >  1:1851675/1‑150 (MQ=255)
                                                                                                                                          aGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACacca        >  1:1851644/1‑150 (MQ=255)
                                                                                                                                               gCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACACCAACCAg   <  2:2190943/150‑1 (MQ=255)
                                                                                                                                               gCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACACCAACCAg   <  1:1384194/150‑1 (MQ=255)
                                                                                                                                                ccAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACACCAACCAgg  >  2:2027356/1‑150 (MQ=255)
                                                                                                                                                ccAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACACCAACCAgg  <  2:19051/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAATATGATGACGATCCTTCCGTACCTAAATCGCCCAGACGGGTGGAAATATAAGGAAACGCGATTACCACACCAACCAGG  >  NC_000913/2129827‑2130120

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