Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,598,924 A→T 15.7% K154N (AAA→AAT)  yjjG → dUMP phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,598,9240AT15.7% 874.5 / 143.5 305K154N (AAA→AAT) yjjGdUMP phosphatase
Reads supporting (aligned to +/- strand):  ref base A (135/122);  new base T (29/19);  total (164/141)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 308 total aligned reads displayed.

GCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGA  >  CP009273/4598778‑4599073
                                                                                                                                                  |                                                                                                                                                     
gCCATTCGAGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAAtt                                                                                                                                                    >  2:553407/1‑150 (MQ=255)
 ccATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATAttt                                                                                                                                                   <  1:2641223/150‑1 (MQ=255)
 ccATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAAttt                                                                                                                                                   <  1:1987154/150‑1 (MQ=255)
    ttCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCg                                                                                                                                                <  1:1587316/150‑1 (MQ=255)
    ttCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCg                                                                                                                                                <  1:3662694/150‑1 (MQ=255)
    ttCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCg                                                                                                                                                <  1:284464/150‑1 (MQ=255)
        tGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTa                                                                                                                                            >  1:4206644/1‑150 (MQ=255)
          gCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATg                                                                                                                                          >  1:1486553/1‑150 (MQ=255)
           cAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTGTCGATTATgc                                                                                                                                         >  2:3558712/1‑150 (MQ=255)
            aaCGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATgcg                                                                                                                                        >  1:1967437/1‑150 (MQ=255)
            aaCGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATgcg                                                                                                                                        >  2:772259/1‑150 (MQ=255)
               gCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTg                                                                                                                                     <  1:2611311/150‑1 (MQ=255)
                       cGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGc                                                                                                                             >  2:1583137/1‑150 (MQ=255)
                          catcatCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCggg                                                                                                                          >  1:4125713/1‑150 (MQ=255)
                           atcatcACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGc                                                                                                                         <  1:1188351/150‑1 (MQ=255)
                            tcatcaCCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCa                                                                                                                        >  2:915467/1‑150 (MQ=255)
                              atcaCCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAAt                                                                                                                      >  1:3159577/1‑150 (MQ=255)
                                 aCCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACAGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCt                                                                                                                   <  1:613301/150‑1 (MQ=255)
                                  ccAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTg                                                                                                                  <  1:2266310/150‑1 (MQ=255)
                                     aCGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGAcc                                                                                                               >  1:3108419/1‑150 (MQ=255)
                                      cGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCg                                                                                                              >  2:2090852/1‑150 (MQ=255)
                                      cGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCg                                                                                                              >  2:4366645/1‑150 (MQ=255)
                                        gCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGtt                                                                                                            >  2:3788707/1‑150 (MQ=255)
                                          tttAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCa                                                                                                          >  1:3476515/1‑150 (MQ=255)
                                          tttAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCa                                                                                                          >  1:2699134/1‑150 (MQ=255)
                                           ttAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCAc                                                                                                         <  2:1102042/150‑1 (MQ=255)
                                               tGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGt                                                                                                     <  1:4383938/150‑1 (MQ=255)
                                               tGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGt                                                                                                     >  1:3814850/1‑150 (MQ=255)
                                                  cTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCt                                                                                                  <  2:341817/150‑1 (MQ=255)
                                                  cTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCt                                                                                                  <  2:2424072/150‑1 (MQ=255)
                                                   ttGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGATg                                                                                                 <  1:1830769/150‑1 (MQ=255)
                                                     gCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAt                                                                                               <  2:1205240/150‑1 (MQ=255)
                                                     gCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCAAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGAt                                                                                               >  2:138585/1‑150 (MQ=255)
                                                       aaCAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATCAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATgg                                                                                             >  2:3625082/1‑150 (MQ=255)
                                                       aaCAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATgg                                                                                             <  2:3601760/150‑1 (MQ=255)
                                                        aCAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGt                                                                                            <  1:3412049/150‑1 (MQ=255)
                                                         cAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGAGCGTTCACGCGTGCTGATGGtt                                                                                           <  2:4313439/150‑1 (MQ=255)
                                                            gTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGc                                                                                        <  2:1869936/150‑1 (MQ=255)
                                                             tGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCg                                                                                       <  2:2390184/150‑1 (MQ=255)
                                                                gTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGaca                                                                                    <  1:4498000/150‑1 (MQ=255)
                                                                  cTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACt                                                                                  <  2:1950939/150‑1 (MQ=255)
                                                                  cTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACt                                                                                  <  2:3900911/150‑1 (MQ=255)
                                                                  cTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACt                                                                                  <  2:3036735/150‑1 (MQ=255)
                                                                   tGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTg                                                                                 >  1:1002499/1‑150 (MQ=255)
                                                                   tGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTg                                                                                 >  1:831381/1‑150 (MQ=255)
                                                                     gAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGcc                                                                               >  2:1438938/1‑150 (MQ=255)
                                                                      aaCGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCg                                                                              <  2:2160873/150‑1 (MQ=255)
                                                                         gCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGt                                                                           <  1:1653948/150‑1 (MQ=255)
                                                                         gCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCNNGt                                                                           <  2:978576/150‑1 (MQ=255)
                                                                          cACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTc                                                                          >  2:1164321/1‑150 (MQ=255)
                                                                          cACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTc                                                                          <  2:2266353/150‑1 (MQ=255)
                                                                            cGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCg                                                                        >  1:3922480/1‑150 (MQ=255)
                                                                              ggCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGat                                                                      <  1:1346315/150‑1 (MQ=255)
                                                                                ccTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGatat                                                                    <  2:3161473/150‑1 (MQ=255)
                                                                                  tGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATtc                                                                  <  1:4187359/150‑1 (MQ=255)
                                                                                     gTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCg                                                               <  2:3893007/150‑1 (MQ=255)
                                                                                      tGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCgg                                                              <  1:1219209/150‑1 (MQ=255)
                                                                                       gATTACTTCGATTTGCTGGTGATTTCCGTAGAAGTTGGCGTTGCCAAACCGTATAAGAAAATTTTCGATTATGCGCTGGAACTGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGt                                                             <  2:2220278/150‑1 (MQ=255)
                                                                                       gATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGt                                                             <  1:3850774/150‑1 (MQ=255)
                                                                                        aTTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTg                                                            >  2:1655119/1‑150 (MQ=255)
                                                                                        aTTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTg                                                            >  2:2353229/1‑150 (MQ=255)
                                                                                          tACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGc                                                          >  2:4187126/1‑150 (MQ=255)
                                                                                          tACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGc                                                          <  1:3261248/150‑1 (MQ=255)
                                                                                            cTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCAt                                                        <  1:2555767/150‑1 (MQ=255)
                                                                                              tCGATTTGCTGGTGATTTCCGAAGATGTTGGAGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACATGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCTATTTTCTCGGTGGCAtaa                                                      <  2:2877880/150‑3 (MQ=255)
                                                                                              tCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCa                                                      <  1:1753673/150‑1 (MQ=255)
                                                                                               cGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCaa                                                     <  2:3483351/150‑1 (MQ=255)
                                                                                               cGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCaa                                                     >  2:2742501/1‑150 (MQ=255)
                                                                                               cGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCaa                                                     >  2:1185932/1‑150 (MQ=255)
                                                                                                 aTTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACg                                                   <  2:845436/150‑1 (MQ=255)
                                                                                                    tGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCg                                                <  1:3882701/150‑1 (MQ=255)
                                                                                                      cTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCggg                                              >  2:642374/1‑150 (MQ=255)
                                                                                                        ggTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCt                                            >  1:410994/1‑150 (MQ=255)
                                                                                                         gTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCtt                                           >  1:543316/1‑150 (MQ=255)
                                                                                                            aTTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCg                                        >  2:4021721/1‑150 (MQ=255)
                                                                                                             tttCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGa                                       >  1:2990325/1‑150 (MQ=255)
                                                                                                              ttCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGAc                                      >  1:1497063/1‑150 (MQ=255)
                                                                                                                ccGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACct                                    <  2:2795160/150‑1 (MQ=255)
                                                                                                                  gaagaaGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACctgc                                  >  2:1041332/1‑150 (MQ=255)
                                                                                                                  gaagaaGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACctgc                                  >  1:1322519/1‑150 (MQ=255)
                                                                                                                  gaagaaGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACctgc                                  >  2:2969651/1‑150 (MQ=255)
                                                                                                                      aagTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTgctggc                              >  1:4293262/1‑150 (MQ=255)
                                                                                                                       agTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTgctggct                             >  2:805954/1‑150 (MQ=255)
                                                                                                                        gTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTgctggctg                            <  2:2468950/150‑1 (MQ=255)
                                                                                                                         ttGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGa                           <  1:4211700/150‑1 (MQ=255)
                                                                                                                         gtGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTCTGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGa                           <  1:2103415/149‑1 (MQ=255)
                                                                                                                           ggCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAAt                         <  2:1262858/150‑1 (MQ=255)
                                                                                                                             cGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGc                       >  2:493578/1‑150 (MQ=255)
                                                                                                                                tGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGcaca                    >  1:1098799/1‑150 (MQ=255)
                                                                                                                                  ccAAACCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACAcc                  >  1:3738047/1‑150 (MQ=255)
                                                                                                                                   cAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCa                 >  2:897576/1‑150 (MQ=255)
                                                                                                                                    aaaCCGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCAt                >  1:143278/1‑150 (MQ=255)
                                                                                                                                    aaaCCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCAt                >  1:2531392/1‑150 (MQ=255)
                                                                                                                                       ccGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATcgc             >  1:2376064/1‑150 (MQ=255)
                                                                                                                                        cGAATAAGAATATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATcgcg            <  2:2685953/150‑1 (MQ=255)
                                                                                                                                        cGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATcgcg            >  2:1384534/1‑150 (MQ=255)
                                                                                                                                          aaTAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAg          >  1:2512673/1‑150 (MQ=255)
                                                                                                                                          aaTAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAg          >  1:4425645/1‑150 (MQ=255)
                                                                                                                                           aTAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGc         >  2:1720836/1‑150 (MQ=255)
                                                                                                                                              aGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAAc      >  1:2298168/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCCATTCGTGGCAACGCCAAAATCGGCATCATCACCAACGGCTTTAGTGCCTTGCAACAGGTGCGTCTGGAACGCACGGGCCTGCGTGATTACTTCGATTTGCTGGTGATTTCCGAAGAAGTTGGCGTTGCCAAACCGAATAAGAAAATTTTCGATTATGCGCTGGAACAGGCGGGCAATCCTGACCGTTCACGCGTGCTGATGGTTGGCGACACTGCCGAGTCCGATATTCTCGGTGGCATCAACGCCGGGCTTGCGACCTGCTGGCTGAATGCACACCATCGCGAGCAACCAGA  >  CP009273/4598778‑4599073

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: ‑