Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,643,955 C→A 39.6% A30A (GCC→GCA dicB → Qin prophage; cell division inhibition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,643,9550CA39.6% 237.1 / 198.5 199A30A (GCC→GCAdicBQin prophage; cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base C (61/59);  new base A (39/39);  total (100/99)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 200 total aligned reads displayed.

TCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAG  >  CP009273/1643812‑1644101
                                                                                                                                               |                                                                                                                                                  
tCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTa                                                                                                                                                <  2:1654688/148‑1 (MQ=255)
tCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTa                                                                                                                                                >  2:2855625/1‑148 (MQ=255)
 ctccTGGAGGCTCGAGGGACTGCATCACAAAATTCTTTGTTGAGGTCGCGATAATGAAAACGTTATTAGCAAAGGTTAATACGTCTGTAGGTTGTTTTGAAATTGGTGTCACTATCAGTATCCCAGTATTTACTGATGATGCTATTat                                                                                                                                               <  2:95872/146‑2 (MQ=255)
 cACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTaa                                                                                                                                               >  1:733598/1‑148 (MQ=255)
 cACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTaa                                                                                                                                               >  2:2885280/1‑148 (MQ=255)
  aCCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAAc                                                                                                                                              >  2:1450245/1‑148 (MQ=255)
      ggAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAaga                                                                                                                                          >  2:286442/1‑148 (MQ=255)
             ccAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACa                                                                                                                                   >  2:757829/1‑148 (MQ=255)
             ccAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACa                                                                                                                                   <  1:579193/148‑1 (MQ=255)
               aGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAg                                                                                                                                 <  2:2075664/148‑1 (MQ=255)
                 gCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGaa                                                                                                                               >  2:1225987/1‑148 (MQ=255)
                 gCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGaa                                                                                                                               <  2:2292030/148‑1 (MQ=255)
                 gCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGaa                                                                                                                               <  1:225350/148‑1 (MQ=255)
                    cTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACgg                                                                                                                            <  1:2237956/148‑1 (MQ=255)
                      gCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACggta                                                                                                                          >  2:1141655/1‑146 (MQ=255)
                          cacaAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTa                                                                                                                      <  1:635713/148‑1 (MQ=255)
                             aaaaTTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATta                                                                                                                   <  2:2053234/148‑1 (MQ=255)
                                  tCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATtaaataa                                                                                                              <  2:449847/148‑1 (MQ=255)
                                   cATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATtaaataaa                                                                                                             <  1:2064139/148‑1 (MQ=255)
                                   cATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATtaaataaa                                                                                                             >  1:56782/1‑148 (MQ=255)
                                     ttgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaat                                                                                                           <  2:1923939/148‑1 (MQ=255)
                                     ttgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaat                                                                                                           <  1:179193/148‑1 (MQ=255)
                                      tgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTaaataaaata                                                                                                          <  1:2414015/148‑1 (MQ=255)
                                        ttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCGGGTACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATg                                                                                                        >  2:3155120/1‑148 (MQ=255)
                                         tgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGc                                                                                                       <  2:1373039/148‑1 (MQ=255)
                                            ggACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCAtt                                                                                                    >  2:1061796/1‑148 (MQ=255)
                                             gACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTg                                                                                                   <  1:378971/148‑1 (MQ=255)
                                                 cgATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTc                                                                                               >  1:2802165/1‑148 (MQ=255)
                                                  gATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCa                                                                                              <  1:892613/148‑1 (MQ=255)
                                                  gATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCa                                                                                              <  1:257456/148‑1 (MQ=255)
                                                     aaTGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATg                                                                                           <  2:333667/148‑1 (MQ=255)
                                                      aTGAAAACTTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGc                                                                                          <  1:3170182/148‑1 (MQ=255)
                                                       tGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCt                                                                                         <  2:848017/148‑1 (MQ=255)
                                                       tGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCt                                                                                         >  2:1415884/1‑148 (MQ=255)
                                                       tGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCt                                                                                         <  2:2209550/148‑1 (MQ=255)
                                                         aaaaCGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTgg                                                                                       <  1:2983514/148‑1 (MQ=255)
                                                         aaaaCGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTgg                                                                                       <  2:535359/148‑1 (MQ=255)
                                                         aaaaCGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATCAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTgg                                                                                       <  2:1336065/148‑1 (MQ=255)
                                                           aaCGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCt                                                                                     >  1:2701220/1‑148 (MQ=255)
                                                             cGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCg                                                                                   >  1:1549445/1‑148 (MQ=255)
                                                             cGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCg                                                                                   <  1:2116831/148‑1 (MQ=255)
                                                               ttattaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGtt                                                                                 >  1:525125/1‑148 (MQ=255)
                                                               ttattaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGtt                                                                                 <  1:287625/148‑1 (MQ=255)
                                                                 attaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTa                                                                               <  2:2717276/148‑1 (MQ=255)
                                                                  ttaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCAGGGTTTCAAGGCTGGCTCGTTGAc                                                                              >  1:3220565/1‑148 (MQ=255)
                                                                   taCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACg                                                                             >  2:2955458/1‑148 (MQ=255)
                                                                   taCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACg                                                                             >  2:712236/1‑148 (MQ=255)
                                                                    aCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGt                                                                            <  1:2526587/148‑1 (MQ=255)
                                                                     ccAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTc                                                                           >  1:614866/1‑148 (MQ=255)
                                                                      cAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCt                                                                          >  1:3139187/1‑148 (MQ=255)
                                                                      cAAACGTTAATAACACTGAAGGTTGTTTTGAAATTGGTGACACTATCAGTAAAACACTATTTACTAAAGATTCCATTAACAATAGAAAACAAGAACGGGAGCTATTAAAAAAAATATGAATTGTTTAAATGATGGCTCGTTTAAGTCt                                                                          >  2:484859/1‑148 (MQ=255)
                                                                         aCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGAt                                                                       >  2:1251764/1‑148 (MQ=255)
                                                                         aCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGAt                                                                       >  1:1976747/1‑148 (MQ=255)
                                                                         aCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGAt                                                                       >  1:2349942/1‑148 (MQ=255)
                                                                            ttAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGcc                                                                    >  2:1122827/1‑148 (MQ=255)
                                                                             tAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACTGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCa                                                                   >  2:2919832/1‑148 (MQ=255)
                                                                             tAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCa                                                                   <  1:2805468/148‑1 (MQ=255)
                                                                             tAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCa                                                                   >  2:2236692/1‑148 (MQ=255)
                                                                                tACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCcaaa                                                                >  2:125262/1‑148 (MQ=255)
                                                                                 aCGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAg                                                               <  2:1725766/148‑1 (MQ=255)
                                                                                 aCGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAg                                                               >  2:2308261/1‑148 (MQ=255)
                                                                                 aCGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAg                                                               >  1:1618919/1‑148 (MQ=255)
                                                                                 aCGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAg                                                               <  1:2942599/148‑1 (MQ=255)
                                                                                       gAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGAtgtg                                                         >  2:2103181/1‑148 (MQ=255)
                                                                                       gAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGAtgtg                                                         >  1:625044/1‑148 (MQ=255)
                                                                                         aGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGca                                                       >  1:391430/1‑148 (MQ=255)
                                                                                          ggTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGcac                                                      >  2:6588/1‑148 (MQ=255)
                                                                                           gttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGcaca                                                     <  2:1022640/148‑1 (MQ=255)
                                                                                             tgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAAt                                                   <  1:1673381/148‑1 (MQ=255)
                                                                                             tgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAAt                                                   <  1:3310509/148‑1 (MQ=255)
                                                                                               ttttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTGAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGa                                                 >  2:943184/1‑148 (MQ=255)
                                                                                               ttttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGa                                                 <  1:857823/148‑1 (MQ=255)
                                                                                                   gAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                             <  1:1029810/148‑1 (MQ=255)
                                                                                                   gAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                             <  1:2885280/148‑1 (MQ=255)
                                                                                                     aaTTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAg                                           >  1:2128110/1‑148 (MQ=255)
                                                                                                     aaTTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAg                                           <  1:2855625/148‑1 (MQ=255)
                                                                                                       ttGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCa                                         >  1:2025331/1‑148 (MQ=255)
                                                                                                          gtgtCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCAttt                                      <  1:2796077/148‑1 (MQ=255)
                                                                                                           tgtCACTATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTtg                                     >  2:2016858/1‑148 (MQ=255)
                                                                                                               aCTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCt                                 <  2:1387076/148‑1 (MQ=255)
                                                                                                               aCTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAACGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCt                                 <  1:184624/148‑1 (MQ=255)
                                                                                                                 tATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTg                               >  2:931863/1‑148 (MQ=255)
                                                                                                                 tATCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTg                               <  1:1148289/148‑1 (MQ=255)
                                                                                                                   tCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGtt                             >  2:1178442/1‑148 (MQ=255)
                                                                                                                   tCAGTAACCCAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGtt                             >  1:505755/1‑148 (MQ=255)
                                                                                                                         aCCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACa                       >  2:3108967/1‑148 (MQ=255)
                                                                                                                          cccAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAg                      <  2:2644534/148‑1 (MQ=255)
                                                                                                                           ccAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGt                     <  2:2198683/148‑1 (MQ=255)
                                                                                                                           ccAGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGt                     <  2:2399695/148‑1 (MQ=255)
                                                                                                                             aGTATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGttt                   >  1:2306410/1‑148 (MQ=255)
                                                                                                                               tATTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGttttt                 >  1:1965524/1‑148 (MQ=255)
                                                                                                                                aTTTACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTc                >  2:1129801/1‑148 (MQ=255)
                                                                                                                                 tttACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCt               >  2:1643228/1‑148 (MQ=255)
                                                                                                                                  ttACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCtt              <  1:703195/148‑1 (MQ=255)
                                                                                                                                  ttACTGAAGATGCAATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCtt              >  2:971425/1‑148 (MQ=255)
                                                                                                                                         aGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCg       >  2:1949460/1‑148 (MQ=255)
                                                                                                                                          gATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCgg      >  1:2991577/1‑148 (MQ=255)
                                                                                                                                            tGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGgag    >  1:800899/1‑148 (MQ=255)
                                                                                                                                             gCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGgaga   >  2:3244999/1‑148 (MQ=255)
                                                                                                                                              ccATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGagag  >  1:127020/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
TCACCTGGAGGCACCAGGCACTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAG  >  CP009273/1643812‑1644101

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: