Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,967,474 G→T 100% A169A (GCG→GCT rffH → glucose‑1‑phosphate thymidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,967,4740GT97.6% 531.4 / 3.6 170A169A (GCG→GCTrffHglucose‑1‑phosphate thymidylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base T (92/74);  total (94/76)
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.60e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 178 total aligned reads displayed.

GTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGC  >  CP009273/3967329‑3967618
                                                                                                                                                 |                                                                                                                                                
gTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGt                                                                                                                                                <  1:3042003/148‑1 (MQ=255)
 tCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCtgtg                                                                                                                                               <  1:1370550/148‑1 (MQ=255)
 tCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCtgtg                                                                                                                                               >  1:335725/1‑148 (MQ=255)
 tCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCtgtg                                                                                                                                               <  1:3139367/148‑1 (MQ=255)
 tCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCtgtg                                                                                                                                               >  2:2103642/1‑148 (MQ=255)
 tCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCtgtg                                                                                                                                               >  1:2171708/1‑148 (MQ=255)
  cATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGa                                                                                                                                              >  2:2534472/1‑148 (MQ=255)
    tGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGAcc                                                                                                                                            >  1:1140154/1‑148 (MQ=255)
       tGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCggg                                                                                                                                         <  2:3130713/148‑1 (MQ=255)
         cggcgCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCt                                                                                                                                       <  1:3523755/148‑1 (MQ=255)
            cgcgcACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTa                                                                                                                                    <  2:1703598/148‑1 (MQ=255)
                     aaGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGa                                                                                                                           <  1:2486139/148‑1 (MQ=255)
                      aGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGAc                                                                                                                          >  2:1815556/1‑148 (MQ=255)
                           ctaCGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACATTTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTaa                                                                                                                     <  1:2645404/146‑1 (MQ=255)
                             aCGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAg                                                                                                                   <  2:27497/148‑1 (MQ=255)
                             aCGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAg                                                                                                                   <  1:2636925/148‑1 (MQ=255)
                              cGGTTTTTGGCTATCAGGTGATGGACCCTGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGt                                                                                                                  <  1:2673364/148‑1 (MQ=255)
                              cGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGt                                                                                                                  >  1:2745262/1‑148 (MQ=255)
                                 tttttGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTACAGTCGt                                                                                                               >  2:2689385/1‑148 (MQ=255)
                                 tttttGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTAGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGt                                                                                                               <  1:3389266/148‑1 (MQ=255)
                                  ttttGGCTATCAGGTGATGGACCCTGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTg                                                                                                              >  1:2730322/1‑148 (MQ=255)
                                   tttGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTgg                                                                                                             <  2:1139405/148‑1 (MQ=255)
                                   tttGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTgg                                                                                                             <  1:2809891/148‑1 (MQ=255)
                                   tttGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTgg                                                                                                             <  2:256910/148‑1 (MQ=255)
                                    ttGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGa                                                                                                            >  1:253797/1‑148 (MQ=255)
                                     tGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAg                                                                                                           <  1:384825/148‑1 (MQ=255)
                                         tATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACg                                                                                                       >  2:3101676/1‑148 (MQ=255)
                                         tATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACg                                                                                                       <  2:1140713/148‑1 (MQ=255)
                                          aTCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGc                                                                                                      >  1:3056230/1‑148 (MQ=255)
                                           tCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCa                                                                                                     >  1:1011939/1‑148 (MQ=255)
                                            cAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCaa                                                                                                    <  2:2029760/148‑1 (MQ=255)
                                              ggTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAg                                                                                                  >  2:1953819/1‑148 (MQ=255)
                                               gTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGc                                                                                                 >  2:2707637/1‑148 (MQ=255)
                                                  aTGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAg                                                                                                  >  2:298937/1‑144 (MQ=255)
                                                   tGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTCTTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGt                                                                                             >  2:3103372/1‑148 (MQ=255)
                                                   tGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGt                                                                                             >  2:116385/1‑148 (MQ=255)
                                                   tGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGt                                                                                             >  2:525958/1‑148 (MQ=255)
                                                    ggACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTg                                                                                            <  2:335725/148‑1 (MQ=255)
                                                    ggACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTg                                                                                            <  2:2171708/148‑1 (MQ=255)
                                                      aCCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGaa                                                                                          >  2:2281377/1‑148 (MQ=255)
                                                        ccGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGc                                                                                        <  1:2428847/148‑1 (MQ=255)
                                                         cGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGcc                                                                                       <  2:683979/148‑1 (MQ=255)
                                                         cGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGcc                                                                                       <  1:3096321/148‑1 (MQ=255)
                                                           gAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGt                                                                                     <  2:1523415/148‑1 (MQ=255)
                                                           gAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGt                                                                                     <  2:1523482/148‑1 (MQ=255)
                                                            aaCGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTAGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGAAGGTGAAGCCGTc                                                                                    >  1:1509564/1‑148 (MQ=255)
                                                                cTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAg                                                                                <  1:1603376/148‑1 (MQ=255)
                                                                   tGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCgt                                                                             >  2:2058830/1‑148 (MQ=255)
                                                                   tGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCgt                                                                             <  1:1846459/148‑1 (MQ=255)
                                                                     gCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTATAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCgtgg                                                                           <  2:3375683/148‑1 (MQ=255)
                                                                         ggtggAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGaa                                                                       <  2:2791597/148‑1 (MQ=255)
                                                                             gAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTgg                                                                   <  1:1784005/148‑1 (MQ=255)
                                                                                tttGACGACAATTTCCGCGCTATCTCGGTGGAATAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgaga                                                                <  2:3316289/148‑1 (MQ=255)
                                                                                tttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgaga                                                                >  1:1333678/1‑148 (MQ=255)
                                                                                tttGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGgaga                                                                >  1:2227876/1‑148 (MQ=255)
                                                                                   gacgacAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTa                                                             <  1:2534472/148‑1 (MQ=255)
                                                                                         aaTTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCa                                                       >  2:931960/1‑148 (MQ=255)
                                                                                          aTTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGAAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCAt                                                      >  1:652048/1‑148 (MQ=255)
                                                                                           tttCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATc                                                     >  1:1038187/1‑148 (MQ=255)
                                                                                            ttCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCa                                                    >  2:2593788/1‑148 (MQ=255)
                                                                                              ccGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAAc                                                  >  2:2497773/1‑148 (MQ=255)
                                                                                                    tATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATg                                            <  2:40771/148‑1 (MQ=255)
                                                                                                     aTCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGt                                           >  2:2307044/1‑148 (MQ=255)
                                                                                                      tctcGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTa                                          <  1:1935481/148‑1 (MQ=255)
                                                                                                       ctcGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGtcc                                         >  2:2097458/1‑146 (MQ=255)
                                                                                                       ctcGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTAc                                         >  2:452663/1‑148 (MQ=255)
                                                                                                         cGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCt                                       >  1:31242/1‑148 (MQ=255)
                                                                                                         cGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCt                                       >  2:3511435/1‑148 (MQ=255)
                                                                                                           cTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCg                                     <  2:1384963/148‑1 (MQ=255)
                                                                                                           cTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCg                                     <  1:455938/148‑1 (MQ=255)
                                                                                                           cTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCg                                     >  1:123158/1‑148 (MQ=255)
                                                                                                            tGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGa                                    <  2:1140154/148‑1 (MQ=255)
                                                                                                             ggAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAg                                   <  2:2425771/148‑1 (MQ=255)
                                                                                                              gaagaaAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAATTCGTGGAGTACTCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAgg                                  >  1:3261063/1‑148 (MQ=255)
                                                                                                              gaagaaAAGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAgg                                  <  1:1095059/148‑1 (MQ=255)
                                                                                                                  aaaaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCggg                              >  2:1992906/1‑148 (MQ=255)
                                                                                                                  aaaaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCggg                              >  2:3041302/1‑148 (MQ=255)
                                                                                                                  aaaaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCggg                              >  1:550000/1‑148 (MQ=255)
                                                                                                                  aaaaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCggg                              >  1:2080191/1‑148 (MQ=255)
                                                                                                                   aaaGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGc                             >  2:3241076/1‑148 (MQ=255)
                                                                                                                   aaaGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGc                             >  2:976335/1‑148 (MQ=255)
                                                                                                                     aGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCaa                           <  1:2185326/148‑1 (MQ=255)
                                                                                                                     aGCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCaa                           >  1:2946920/1‑148 (MQ=255)
                                                                                                                      gCCAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAAc                          <  2:864243/148‑1 (MQ=255)
                                                                                                                        cAAAACAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCt                        >  1:9523/1‑148 (MQ=255)
                                                                                                                            aCAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGAcc                    >  2:1049234/1‑148 (MQ=255)
                                                                                                                             cAGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCg                   >  2:2075430/1‑148 (MQ=255)
                                                                                                                              aGCCGAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGt                  <  2:2153678/148‑1 (MQ=255)
                                                                                                                                  gAAGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGaa              <  1:3432737/148‑1 (MQ=255)
                                                                                                                                   aaGTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAAc             >  1:2953928/1‑148 (MQ=255)
                                                                                                                                     gTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTg           <  1:2605530/148‑1 (MQ=255)
                                                                                                                                     gTCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTg           >  1:3502070/1‑148 (MQ=255)
                                                                                                                                      tCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGc          <  2:3195297/148‑1 (MQ=255)
                                                                                                                                      tCAAACTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGc          <  2:1741825/148‑1 (MQ=255)
                                                                                                                                        aaaaTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAACGTCGTGGAGTACGCAAAGCATGTGAAGCCGTCGGAGCGTGTTGAACTGTAGATTACATCCATCAACCAGATGTCCCTCGAGGCGGGCAACCTGACCGTTGAACTGCta        >  2:367287/1‑147 (MQ=255)
                                                                                                                                        aaaCTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTc        >  2:3523335/1‑148 (MQ=255)
                                                                                                                                        aaaCTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTc        <  1:2281377/148‑1 (MQ=255)
                                                                                                                                          aCTGGGCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCgg      <  2:443365/148‑1 (MQ=255)
                                                                                                                                              ggCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGgcgc  >  1:3373111/1‑148 (MQ=255)
                                                                                                                                              ggCTGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGgcgc  >  1:735801/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTCATGTTGCGGCGCGCACCGAAGGGGCGACGGTTTTTGGCTATCAGGTGATGGACCCGGAACGCTTTGGCGTGGTGGAGTTTGACGACAATTTCCGCGCTATCTCGCTGGAAGAAAAGCCAAAACAGCCGAAGTCAAACTGGGCGGTGACCGGGCTTTATTTCTACGACAGTAAAGTCGTGGAGTACGCAAAGCAGGTGAAGCCGTCGGAGCGTGGTGAACTGGAGATTACCTCCATCAACCAGATGTACCTCGAGGCGGGCAACCTGACCGTTGAACTGCTCGGGCGC  >  CP009273/3967329‑3967618

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: