Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 992,767 C→A 89.4% D68Y (GAT→TAT)  ssuE ← NAD(P)H‑dependent FMN reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273992,7670CA89.4% 442.9 / 48.3 177D68Y (GAT→TAT) ssuENAD(P)H‑dependent FMN reductase
Reads supporting (aligned to +/- strand):  ref base C (10/9);  new base A (83/75);  total (93/84)
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.71e-01

Only 100 of 177 total aligned reads displayed.

CGCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCC  >  CP009273/992621‑992909
                                                                                                                                                  |                                                                                                                                              
cGCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCg                                                                                                                                               >  2:1815127/1‑148 (MQ=255)
  ccAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGc                                                                                                                                             <  2:1084322/148‑1 (MQ=255)
   cagcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGcc                                                                                                                                            >  1:1939537/1‑148 (MQ=255)
    agcagATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACACTCAGCCCATAGGCCt                                                                                                                                           >  1:1750240/1‑148 (MQ=255)
       agATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGtt                                                                                                                                        <  2:1290599/148‑1 (MQ=255)
         aTGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGc                                                                                                                                      >  1:402803/1‑148 (MQ=255)
                aCGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGtt                                                                                                                               >  1:2585452/1‑148 (MQ=255)
                 cGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTc                                                                                                                              >  2:360493/1‑148 (MQ=255)
                 cGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTc                                                                                                                              >  1:2207066/1‑148 (MQ=255)
                         cccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAgg                                                                                                                      >  2:470158/1‑148 (MQ=255)
                         cccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAgg                                                                                                                      <  1:2028722/148‑1 (MQ=255)
                         cccGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAgg                                                                                                                      <  1:962128/148‑1 (MQ=255)
                             tCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCtt                                                                                                                  >  2:1211219/1‑148 (MQ=255)
                              cGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTg                                                                                                                 <  1:2406552/148‑1 (MQ=255)
                               gCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGa                                                                                                                <  1:2980570/148‑1 (MQ=255)
                                ccAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAg                                                                                                               <  2:2268775/148‑1 (MQ=255)
                                  aGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTGCGGTGAAGGTCTTGAGTg                                                                                                             >  2:716615/1‑148 (MQ=255)
                                    cGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGcc                                                                                                           <  2:1233199/148‑1 (MQ=255)
                                         gCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACt                                                                                                      >  1:1046894/1‑148 (MQ=255)
                                          caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTa                                                                                                     <  1:2012935/148‑1 (MQ=255)
                                          caccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTa                                                                                                     <  1:2234365/148‑1 (MQ=255)
                                           accacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTAt                                                                                                    >  1:3118238/1‑148 (MQ=255)
                                            ccacTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATc                                                                                                   >  2:2591908/1‑148 (MQ=255)
                                              acTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGa                                                                                                 >  2:77565/1‑148 (MQ=255)
                                               cTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGaa                                                                                                >  1:3260602/1‑148 (MQ=255)
                                               cTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGaa                                                                                                <  1:1263591/148‑1 (MQ=255)
                                                 ttGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAAc                                                                                              >  2:33825/1‑148 (MQ=255)
                                                  tGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACg                                                                                             <  2:1736576/148‑1 (MQ=255)
                                                     cTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGc                                                                                          <  1:486274/148‑1 (MQ=255)
                                                      ttGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCa                                                                                         >  2:887826/1‑148 (MQ=255)
                                                          aaaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaaa                                                                                     <  2:529003/148‑1 (MQ=255)
                                                          aaaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaaa                                                                                     <  2:585242/148‑1 (MQ=255)
                                                          aaaGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATaaa                                                                                     <  2:585853/148‑1 (MQ=255)
                                                            aGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGt                                                                                   >  1:1304231/1‑148 (MQ=255)
                                                            aGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGt                                                                                   >  2:2975012/1‑148 (MQ=255)
                                                            aGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGt                                                                                   >  2:2975008/1‑148 (MQ=255)
                                                                cgTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGAt                                                                               >  2:3256576/1‑148 (MQ=255)
                                                                  ttCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCt                                                                             >  1:1973884/1‑148 (MQ=255)
                                                                  ttCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCt                                                                             <  2:2387670/148‑1 (MQ=255)
                                                                   tCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCtt                                                                            >  2:1264364/1‑148 (MQ=255)
                                                                   tCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCtt                                                                            >  2:1264439/1‑148 (MQ=255)
                                                                    cTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTc                                                                           <  2:127336/148‑1 (MQ=255)
                                                                       gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                        <  1:1646742/148‑1 (MQ=255)
                                                                       gcagcaGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCATGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCgg                                                                        >  2:2104972/1‑148 (MQ=255)
                                                                          gcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGc                                                                     >  2:1843964/1‑148 (MQ=255)
                                                                          gcagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGc                                                                     >  2:1843920/1‑148 (MQ=255)
                                                                           cagGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCg                                                                    <  1:2884562/148‑1 (MQ=255)
                                                                             ggTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGaa                                                                  <  1:3295668/148‑1 (MQ=255)
                                                                              gTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAg                                                                 >  1:178071/1‑148 (MQ=255)
                                                                              gTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAg                                                                 >  1:1848918/1‑148 (MQ=255)
                                                                              gTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAg                                                                 >  2:1445140/1‑148 (MQ=255)
                                                                               tCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGt                                                                <  2:1133235/148‑1 (MQ=255)
                                                                                cGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtt                                                               <  1:686101/148‑1 (MQ=255)
                                                                                cGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGtt                                                               <  1:2344713/148‑1 (MQ=255)
                                                                                   gCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTg                                                            <  1:779182/148‑1 (MQ=255)
                                                                                    cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGc                                                           <  2:1379554/148‑1 (MQ=255)
                                                                                    cAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGc                                                           >  1:162217/1‑148 (MQ=255)
                                                                                     aGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCa                                                          >  2:1048799/1‑148 (MQ=255)
                                                                                      gggTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                         >  2:1916940/1‑148 (MQ=255)
                                                                                      gggTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAg                                                         <  2:63501/148‑1 (MQ=255)
                                                                                         ttttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGAtt                                                      <  2:204690/148‑1 (MQ=255)
                                                                                         ttttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGAtt                                                      >  1:2278923/1‑148 (MQ=255)
                                                                                          tttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTc                                                     >  1:665218/1‑148 (MQ=255)
                                                                                          tttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTc                                                     <  2:3005235/148‑1 (MQ=255)
                                                                                           ttCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTcc                                                    <  1:1054625/148‑1 (MQ=255)
                                                                                             cAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAg                                                  <  1:1937844/148‑1 (MQ=255)
                                                                                             cAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAg                                                  <  2:1552547/148‑1 (MQ=255)
                                                                                             cAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAg                                                  <  2:1004953/148‑1 (MQ=255)
                                                                                               aCGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTg                                                <  1:738432/148‑1 (MQ=255)
                                                                                                cGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGa                                               <  1:3287851/148‑1 (MQ=255)
                                                                                                 gCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGAt                                              <  2:1608745/148‑1 (MQ=255)
                                                                                                  cACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATa                                             >  1:862982/1‑148 (MQ=255)
                                                                                                    cctgAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTTAGTGCCGGACTATCTAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaaa                                           <  1:2403248/145‑1 (MQ=255)
                                                                                                    ccGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATaaa                                           >  2:2425787/1‑148 (MQ=255)
                                                                                                       gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACct                                        <  2:1939537/148‑1 (MQ=255)
                                                                                                       gAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACct                                        <  1:687107/148‑1 (MQ=255)
                                                                                                        aaTAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACctc                                       >  2:2644030/1‑148 (MQ=255)
                                                                                                          tAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTa                                     >  2:2173631/1‑148 (MQ=255)
                                                                                                           aGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTAc                                    >  2:651764/1‑148 (MQ=255)
                                                                                                             gcggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACAt                                  <  1:778991/148‑1 (MQ=255)
                                                                                                               ggcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATcc                                <  1:1471363/148‑1 (MQ=255)
                                                                                                                gcTTTATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCa                               <  1:33825/148‑1 (MQ=255)
                                                                                                                   ttATACACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGc                            <  2:833256/148‑1 (MQ=255)
                                                                                                                      taCACAGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCAt                         <  1:1170740/148‑1 (MQ=255)
                                                                                                                         acaGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTa                      >  2:1292167/1‑148 (MQ=255)
                                                                                                                          caGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAg                     <  1:269851/148‑1 (MQ=255)
                                                                                                                          caGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAg                     <  1:269793/148‑1 (MQ=255)
                                                                                                                          caGGCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAg                     >  1:1296965/1‑148 (MQ=255)
                                                                                                                            ggCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGtt                   >  1:1244328/1‑148 (MQ=255)
                                                                                                                            ggCGTGGAGACAATCAGCCCATAGGCCTTTTGCAGCTGTTCGGTGAAGGTCTTGAGTGACTGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGAGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGtt                   <  2:2416606/148‑1 (MQ=255)
                                                                                                                             gCGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGttt                  >  2:1661819/1‑148 (MQ=255)
                                                                                                                              cGTGGCGACAATCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGtttt                 <  2:3097916/148‑1 (MQ=255)
                                                                                                                                        aTCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCat       >  1:1382915/1‑148 (MQ=255)
                                                                                                                                        aTCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCat       <  2:402803/148‑1 (MQ=255)
                                                                                                                                         tCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata      >  1:2025998/1‑148 (MQ=255)
                                                                                                                                         tCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata      >  1:238226/1‑148 (MQ=255)
                                                                                                                                         tCAGCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCata      >  1:25198/1‑148 (MQ=255)
                                                                                                                                            gCCCATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTc   >  2:1449160/1‑148 (MQ=255)
                                                                                                                                             cccATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTcc  >  1:1627071/1‑148 (MQ=255)
                                                                                                                                             cccATAGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTcc  >  1:3224511/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
CGCCAGCAGATGGGCCACGGTACCGCCCGTCGCCAGCGGTAGCACCACTTTGCCTTGCAAAGCGCGTTCTGGCAGCAGGTCGAGCAGGGTTTTCAACGCACCGGAATAGGCGGCTTTATACACAGGCGTGGCGACAATCAGCCCATCGGCCTGTTGCAGCTGTTCGGTGAAGGTCTTGAGTGCCGGACTATCGAAACGAGCATAAAGTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCC  >  CP009273/992621‑992909

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