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

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

Only 100 of 111 total aligned reads displayed.

CGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAAT  >  NC_000913/2129754‑2130043
                                                                                                                                                     |                                                                                                                                            
cGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCAtt                                                                                                                                              <  2:1483229/150‑1 (MQ=255)
cGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCAAAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCAtt                                                                                                                                              >  2:151023/1‑150 (MQ=255)
 gATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTg                                                                                                                                             <  2:835764/150‑1 (MQ=255)
 gATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTg                                                                                                                                             <  1:59323/150‑1 (MQ=255)
      cGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAACACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                                                                                        >  2:697472/1‑150 (MQ=255)
      cGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                                                                                        >  2:458939/1‑150 (MQ=255)
      cGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTAt                                                                                                                                        >  2:577984/1‑150 (MQ=255)
       ggAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                                                                                                       <  2:739369/150‑1 (MQ=255)
       ggAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATg                                                                                                                                       <  2:1136805/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)
            aaaaaGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttg                                                                                                                                  <  2:724290/150‑2 (MQ=255)
            aaaaaGATATAATTGCACGAGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttg                                                                                                                                  <  2:1498319/150‑2 (MQ=255)
            aaaaaGATATAATAGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTTATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACAGCGCGTACCAGGAAAAATAAATAATCAGCAGTTACAAACCATAGACAATGGtttg                                                                                                                                  <  2:461192/150‑2 (MQ=255)
                aGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCg                                                                                                                              <  1:501736/150‑1 (MQ=255)
                aGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATAGACAATGGtttgtatc                                                                                                                              <  1:477356/150‑6 (MQ=255)
                  atatAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAc                                                                                                                            <  2:1262615/150‑1 (MQ=255)
                       aTTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCg                                                                                                                       <  2:834576/150‑1 (MQ=255)
                           cACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCAcc                                                                                                                   >  2:366013/1‑150 (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:1085511/1‑150 (MQ=255)
                                       gCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                                                                                                       <  2:930111/150‑1 (MQ=255)
                                       gCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTcc                                                                                                       <  2:1016624/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)
                                                   aaaaGCCTTGATCATCTTTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAAGCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGATAAATATTCCGAATGATGCGAc                                                                                           <  1:1529317/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:1063881/150‑1 (MQ=255)
                                                              tcatcaTTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATaaa                                                                                <  2:1592499/150‑1 (MQ=255)
                                                                atcatTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGt                                                                              <  2:593598/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:1131452/150‑1 (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)
                                                                                      tGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGt                                                        >  2:1654454/1‑150 (MQ=255)
                                                                                       gTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTc                                                       >  1:1513903/1‑150 (MQ=255)
                                                                                       gTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTc                                                       <  1:1393368/150‑1 (MQ=255)
                                                                                        tACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCa                                                      <  2:1146090/150‑1 (MQ=255)
                                                                                            acaGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAatat                                                  >  2:1007507/1‑150 (MQ=255)
                                                                                             caGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATg                                                 <  1:577984/150‑1 (MQ=255)
                                                                                             caGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATg                                                 <  1:458939/150‑1 (MQ=255)
                                                                                             caGTGAAAAAAACACCGCGAACCAGGAAAAATAAATAATCTGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAATTCTGACTACACCGTCAATATg                                                 <  1:795069/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)
                                                                                                       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:1683347/150‑1 (MQ=255)
                                                                                                              cgAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGa                                >  1:693968/1‑150 (MQ=255)
                                                                                                                 aCCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCa                             <  2:580779/150‑1 (MQ=255)
                                                                                                                  ccAGGAAGAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGAAATCCGCACCGGTAAATATTCCGAATGATGCGCCATATTCATATAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCAt                            >  1:1244908/1‑150 (MQ=255)
                                                                                                                   cAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATc                           <  1:1659344/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)
                                                                                                                                                     |                                                                                                                                            
CGATACCGGAGTAAAAAGATATAATTGCACGCGAAAATTGCGGTTATCACCAAAAGCGTTGATCATCATTTTTATCACTTTCCCCATGTACCACAGTGATAAAAACACCGCGAACCAGGAAAAATAAATAATCAGCAGATACAAACCATTGTCTATGGTTTTTCCGACATCCGCACCGTTAAATATTCCGAATGATGCGACATATTCATAAAGTGAGCCAAATCTGACTACACCGTCAATATGGGTCAAAGAATAACCGACCATCACCAGCGGACCGACAATACGATAAT  >  NC_000913/2129754‑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: