Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,129,905 T→A 6.5% D325V (GAC→GTC)  wcaD ← putative colanic acid polymerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,129,9050TA 6.5% 303.9 / 8.2 107D325V (GAC→GTC) wcaDputative colanic acid polymerase
Reads supporting (aligned to +/- strand):  ref base T (41/59);  new base A (3/4);  total (44/63)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 107 total aligned reads displayed.

CGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAAT  >  NC_000913/2129760‑2130043
                                                                                                                                                 |                                                                                                                                          
cGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAACACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                                                                                        >  2:697472/1‑150 (MQ=255)
cGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                                                                                        >  2:458939/1‑150 (MQ=255)
cGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                                                                                        >  2:577984/1‑150 (MQ=255)
 ggAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                                                                                                       <  2:1136805/150‑1 (MQ=255)
 ggAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                                                                                                       <  2:739369/150‑1 (MQ=255)
 cgAGTAAAAAGATATAATTTCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCAACATTTTTATCACTTTGCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                                                                                                       <  2:839143/149‑1 (MQ=255)
    gTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACAGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtt                                                                                                                                    >  1:670378/1‑150 (MQ=255)
    gTAAAAAGATAAAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGtt                                                                                                                                    >  1:593598/1‑150 (MQ=255)
      aaaaaGATATAATTGCCCGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGttgt                                                                                                                                  >  2:1107684/1‑148 (MQ=255)
      aaaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttg                                                                                                                                  <  2:724290/150‑2 (MQ=255)
      aaaaaGATATAATTGCACGAGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttg                                                                                                                                  <  2:1498319/150‑2 (MQ=255)
      aaaaaGATATAATAGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTTATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACAGCGCGTACCAGGAAAAATAAATAATCAGCAGTTACAAACCATAGACAATGGtttg                                                                                                                                  <  2:461192/150‑2 (MQ=255)
         aaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTcc                                                                                                                               <  2:545063/150‑1 (MQ=255)
          aGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCg                                                                                                                              <  1:501736/150‑1 (MQ=255)
          aGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgtatc                                                                                                                              <  1:477356/150‑6 (MQ=255)
            atatAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAc                                                                                                                            <  2:1262615/150‑1 (MQ=255)
               taATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATc                                                                                                                         <  1:918400/150‑1 (MQ=255)
                 aTTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCg                                                                                                                       <  2:834576/150‑1 (MQ=255)
                    gCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAc                                                                                                                    >  1:599889/1‑150 (MQ=255)
                     cACGGGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTTATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                                                                                                   <  2:769049/150‑1 (MQ=255)
                     cACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                                                                                                   >  2:366013/1‑150 (MQ=255)
                       cgcgAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACACACCATTGTCTATGGTTTTTCCGACATCCGCACCGt                                                                                                                 <  1:1608308/150‑1 (MQ=255)
                       cgcgAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGt                                                                                                                 <  2:395210/150‑1 (MQ=255)
                        gcgAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGtt                                                                                                                >  2:952450/1‑150 (MQ=255)
                         cgAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTa                                                                                                               >  1:1146090/1‑150 (MQ=255)
                         cgAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTa                                                                                                               <  1:1078695/150‑1 (MQ=255)
                         cgAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTa                                                                                                               >  1:544217/1‑150 (MQ=255)
                          gAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTaa                                                                                                              <  1:866119/150‑1 (MQ=255)
                                 gCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                                                                                                       <  2:1016624/150‑1 (MQ=255)
                                 gCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                                                                                                       >  2:1085511/1‑150 (MQ=255)
                                 gCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                                                                                                       <  2:930111/150‑1 (MQ=255)
                                 gCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTACCCATGTACCACAGTGATAAAAAAACCGCGAACCAGGAAAAATAAATAATCATCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAAAATTcc                                                                                                       <  2:46811/150‑1 (MQ=255)
                                  cGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAATTAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCg                                                                                                      <  2:1273254/150‑1 (MQ=255)
                                  cGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCg                                                                                                      <  2:324014/150‑1 (MQ=255)
                                      tATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCACATGTAACACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCACATACAAACCATTGTCTATGGTTTTTCCGACATCCGAACCGTTAAATATTACGAatg                                                                                                  >  2:1462144/1‑150 (MQ=255)
                                          aCCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGc                                                                                              <  2:871181/150‑1 (MQ=255)
                                                  cGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATAtt                                                                                      >  1:7471/1‑150 (MQ=255)
                                                    ttGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCa                                                                                    <  1:1107684/150‑1 (MQ=255)
                                                     tGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATTTTCAt                                                                                   >  1:928970/1‑150 (MQ=255)
                                                     tGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCAt                                                                                   >  2:1135591/1‑150 (MQ=255)
                                                      gATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATa                                                                                  <  1:952450/150‑1 (MQ=255)
                                                      gATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACAGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATAATAAGAATGAAGCGACATATTCATa                                                                                  >  1:1582754/1‑150 (MQ=255)
                                                        tcatcaTTTTTATCACTTTCCCCATGTCCCACAGTGATAAAAACACCGCGAACCATGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATaaa                                                                                <  2:1058311/150‑1 (MQ=255)
                                                        tcatcaTTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATaaa                                                                                <  2:1592499/150‑1 (MQ=255)
                                                        tcatcaTTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATaaa                                                                                <  2:1063881/150‑1 (MQ=255)
                                                          atcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCAACATCCGCACCGTTAAATATTCCGAatg                                                                                                  >  2:1587795/1‑130 (MQ=255)
                                                          atcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCAACATCCGCACCGTTAAATATTCCGAatg                                                                                                  <  1:1587795/130‑1 (MQ=255)
                                                            catTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCTGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGa                                                                            >  2:1562506/1‑150 (MQ=255)
                                                             atTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAg                                                                           <  1:866013/150‑1 (MQ=255)
                                                              tttttATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGc                                                                          <  1:381052/150‑1 (MQ=255)
                                                                 ttATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATTTTCCGAATGATGCGACATATTCATAAAGTGAGCCaa                                                                       >  1:1344325/1‑150 (MQ=255)
                                                                  tATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCaaa                                                                      <  1:1085511/150‑1 (MQ=255)
                                                                   aTCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAAt                                                                     >  1:170780/1‑150 (MQ=255)
                                                                    tCACTTTCCCCATGTACCACAGTTATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATc                                                                    <  1:1691838/150‑1 (MQ=255)
                                                                       cTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGa                                                                 >  1:202251/1‑150 (MQ=255)
                                                                         ttCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACt                                                               <  2:7471/150‑1 (MQ=255)
                                                                           ccccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTac                                                             <  2:1249368/150‑1 (MQ=255)
                                                                             ccATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTacac                                                           <  1:1561098/150‑1 (MQ=255)
                                                                               aTGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCg                                                         <  1:767812/150‑1 (MQ=255)
                                                                                 ttACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTc                                                       <  2:1693954/149‑1 (MQ=255)
                                                                                 gTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTc                                                       <  1:1393368/150‑1 (MQ=255)
                                                                                 gTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTc                                                       >  1:1513903/1‑150 (MQ=255)
                                                                                  tACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCa                                                      <  2:1146090/150‑1 (MQ=255)
                                                                                      acaGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAatat                                                  >  2:1007507/1‑150 (MQ=255)
                                                                                       caGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATg                                                 <  1:458939/150‑1 (MQ=255)
                                                                                       caGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATg                                                 <  1:577984/150‑1 (MQ=255)
                                                                                        aGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATgg                                                <  2:799753/150‑1 (MQ=255)
                                                                                         gTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATggg                                               <  2:602313/150‑1 (MQ=255)
                                                                                            aTAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCa                                            <  2:202251/150‑1 (MQ=255)
                                                                                              aaaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCAATGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGAAATATTCATAAAGTGAGCCAAATCTGACAACAACGTCAATATGGGTCaaa                                          >  2:1100979/1‑150 (MQ=255)
                                                                                               aaaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACAACGTCAATATGGGTCaaat                                         >  2:476798/1‑149 (MQ=255)
                                                                                                aaaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGa                                        >  2:1277121/1‑150 (MQ=255)
                                                                                                 aaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGaa                                       <  2:112045/150‑1 (MQ=255)
                                                                                                 aaCACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGaa                                       <  1:1384203/150‑1 (MQ=255)
                                                                                                   cacCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATa                                     <  1:1333319/150‑1 (MQ=255)
                                                                                                      cgcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATTAAGTGAGCCATATCTGACTACACCGGCAATATGCGTCAAAGAATAAcc                                  >  2:1052585/1‑150 (MQ=255)
                                                                                                       gcgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCg                                 >  1:661013/1‑150 (MQ=255)
                                                                                                        cgtaCCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGa                                <  1:388995/147‑1 (MQ=255)
                                                                                                        cgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGa                                >  1:693968/1‑150 (MQ=255)
                                                                                                        cgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGa                                <  1:1683347/150‑1 (MQ=255)
                                                                                                           aCCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCa                             <  2:580779/150‑1 (MQ=255)
                                                                                                            ccAGGAAGAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAAATCCGCACCGGTAAATATTCCGAATGATGCGCCATATTCATATAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCAt                            >  1:1244908/1‑150 (MQ=255)
                                                                                                             cAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATc                           <  1:1659344/150‑1 (MQ=255)
                                                                                                               ggAAAAATAAATAATCAGCAGATACAAACCTTTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAc                         <  1:1462144/150‑1 (MQ=255)
                                                                                                               ggAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCAc                         >  1:886840/1‑150 (MQ=255)
                                                                                                                 aaaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCa                       >  2:574762/1‑150 (MQ=255)
                                                                                                                  aaaaTAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAg                      >  1:1527068/1‑150 (MQ=255)
                                                                                                                   aaataaatAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGc                     <  2:894410/150‑1 (MQ=255)
                                                                                                                    aataaataATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCg                    <  1:1464630/150‑1 (MQ=255)
                                                                                                                      taaataaTCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAATGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGa                  >  1:1582276/1‑150 (MQ=255)
                                                                                                                       aaataaTCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGAc                 <  1:1135591/150‑1 (MQ=255)
                                                                                                                        aataatCAGCAGATACAAACCATTGTCTATTGTTTTTCCGACTTCCTCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGAcc                <  1:151023/150‑1 (MQ=255)
                                                                                                                        aataatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGAcc                >  2:1626575/1‑150 (MQ=255)
                                                                                                                         ataatCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCg               <  2:170780/150‑1 (MQ=255)
                                                                                                                                gcagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACAATCACCAGCGGACAGACAATAc        >  1:224463/1‑150 (MQ=255)
                                                                                                                                 cagATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACg       >  1:569897/1‑150 (MQ=255)
                                                                                                                                  tcggaaaaaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCCGCGGACCGACAATACGa      >  1:462813/7‑150 (MQ=255)
                                                                                                                                  tcggaaaaaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGa      >  2:989227/7‑150 (MQ=255)
                                                                                                                                    ggaaaaaCCATAGACAATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGata    >  2:1187977/5‑150 (MQ=255)
                                                                                                                                      aCAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGataat  >  1:326695/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                          
CGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAAT  >  NC_000913/2129760‑2130043

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