Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,097,262 C→T 100% F297F (TTC→TTT mutY → adenine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,097,2620CT95.9% 470.5 / 8.0 147F297F (TTC→TTTmutYadenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base C (0/6);  new base T (62/79);  total (62/85)
Fisher's exact test for biased strand distribution p-value = 3.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 148 total aligned reads displayed.

TCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGC  >  CP009273/3097118‑3097406
                                                                                                                                                |                                                                                                                                                
tCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCCTTTTCAc                                                                                                                                               >  1:2735250/1‑148 (MQ=255)
 ccgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACt                                                                                                                                              >  2:765671/1‑148 (MQ=255)
 ccgccgAGCGGATTGTGGGGCGGTTTATACTGTTACCCGCAGTTTGCCGACGAAGAAAGTTTGCTGCAGTGCCTGGAGCAACGGCAGATTGCTGCCGATAACATGACGCAACTGACCGCGTTTCGGCATACATTCAGCCATTTTCACt                                                                                                                                              >  2:443364/1‑148 (MQ=255)
  cgccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACtt                                                                                                                                             >  2:132343/1‑148 (MQ=255)
   gccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTa                                                                                                                                            <  2:1277160/148‑1 (MQ=255)
    ccgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAg                                                                                                                                           <  1:1443496/148‑1 (MQ=255)
     cgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGa                                                                                                                                          <  2:2956773/148‑1 (MQ=255)
     cgAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTCTCGGCATACCTTCAGCCATTTTCACTTAGa                                                                                                                                          <  2:2888800/148‑1 (MQ=255)
      gAGCGGATTGTGGTGCGGTTTATACTGTTTCGCGCAGTTTGCCGACGAAGAAAGTTTGCGTCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGat                                                                                                                                         <  2:2889666/148‑1 (MQ=255)
      gAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGat                                                                                                                                         <  2:2291267/148‑1 (MQ=255)
        gCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGatat                                                                                                                                       <  1:2145983/148‑1 (MQ=255)
            aTTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATtgtg                                                                                                                                   >  2:1922808/1‑148 (MQ=255)
              tgtgGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGcc                                                                                                                                 <  2:745087/148‑1 (MQ=255)
                 ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAt                                                                                                                              <  1:1821635/148‑1 (MQ=255)
                 ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAt                                                                                                                              <  1:1858294/148‑1 (MQ=255)
                 ggggCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTAt                                                                                                                              <  2:1444855/148‑1 (MQ=255)
                    gCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATGTTCACTTAGATAGTGTGCCTAggtg                                                                                                                           >  2:1497353/1‑148 (MQ=255)
                     cGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTgg                                                                                                                          <  1:2843698/148‑1 (MQ=255)
                       gTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCt                                                                                                                        <  2:209400/148‑1 (MQ=255)
                       gTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCt                                                                                                                        <  2:1155558/148‑1 (MQ=255)
                         ttATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTc                                                                                                                      <  1:2156785/148‑1 (MQ=255)
                         ttATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTc                                                                                                                      >  2:2424059/1‑148 (MQ=255)
                         ttATACTGTGTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTc                                                                                                                      <  1:3394983/148‑1 (MQ=255)
                              cTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCgtg                                                                                                                 <  1:1325841/148‑1 (MQ=255)
                              cTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCgtg                                                                                                                 >  2:194781/1‑148 (MQ=255)
                              cTGTGTCCCGCAGTTTGCCGTCGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTGTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCGTCCCgtg                                                                                                                 <  1:2304970/148‑1 (MQ=255)
                               tGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCgtgt                                                                                                                <  1:398972/148‑1 (MQ=255)
                                gTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtc                                                                                                               >  1:2571555/1‑148 (MQ=255)
                                gTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtc                                                                                                               >  1:2571567/1‑148 (MQ=255)
                                   tCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTGTCGGCATACCTTCAGCCAGTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTgtcgtc                                                                                                            >  2:2978601/1‑148 (MQ=255)
                                    cccGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCa                                                                                                           <  2:242890/148‑1 (MQ=255)
                                      cGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCAtt                                                                                                         >  1:358105/1‑148 (MQ=255)
                                      cGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCAtt                                                                                                         <  1:1450412/148‑1 (MQ=255)
                                        cAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGACAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCa                                                                                                       >  2:1129654/1‑148 (MQ=255)
                                            ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCgg                                                                                                   <  1:483967/148‑1 (MQ=255)
                                            ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCgg                                                                                                   <  1:2998235/148‑1 (MQ=255)
                                            ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCgg                                                                                                   >  1:1112506/1‑148 (MQ=255)
                                            ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCgg                                                                                                   <  1:978189/148‑1 (MQ=255)
                                            ttGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCgg                                                                                                   <  1:2630507/148‑1 (MQ=255)
                                                  acgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCAt                                                                                             <  2:2080152/148‑1 (MQ=255)
                                                   cgaAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATg                                                                                            >  1:2748583/1‑148 (MQ=255)
                                                    gaagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATgg                                                                                           >  1:639966/1‑148 (MQ=255)
                                                    gaagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATCGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATgg                                                                                           <  2:2596797/148‑1 (MQ=255)
                                                     aagatAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGa                                                                                          <  1:2664375/148‑1 (MQ=255)
                                                     aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGa                                                                                          <  2:2705954/148‑1 (MQ=255)
                                                     aagaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGa                                                                                          <  1:1239115/148‑1 (MQ=255)
                                                      agaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGAt                                                                                         >  1:3063287/1‑148 (MQ=255)
                                                       gaaAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATg                                                                                        <  2:3043007/148‑1 (MQ=255)
                                                        aaaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGa                                                                                       >  2:2833009/1‑148 (MQ=255)
                                                         aaGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGaa                                                                                      >  1:2643147/1‑148 (MQ=255)
                                                          aGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAg                                                                                     >  2:2269985/1‑148 (MQ=255)
                                                          aGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAg                                                                                     >  2:144340/1‑148 (MQ=255)
                                                           gTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAgg                                                                                    <  2:3373057/148‑1 (MQ=255)
                                                             ttGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCa                                                                                  >  1:1877430/1‑148 (MQ=255)
                                                              tGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCaa                                                                                 <  2:979828/148‑1 (MQ=255)
                                                              tGCGGCAGTGGCTGGCGAAACGGCAGATTGCTGCCGATCACCTGAAGCAACTGACCGCGTTTCGGCTTACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTTTCGTCATTCACCGGCTGCATGGATGAAGGCaa                                                                                 <  2:266746/148‑1 (MQ=255)
                                                               gCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAAt                                                                                <  2:679270/148‑1 (MQ=255)
                                                                   cAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgcgc                                                                            >  2:1751400/1‑148 (MQ=255)
                                                                   cAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATgcgc                                                                            >  2:1783583/1‑148 (MQ=255)
                                                                         ctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGt                                                                      <  1:2390819/148‑1 (MQ=255)
                                                                         ctggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGt                                                                      >  1:694967/1‑148 (MQ=255)
                                                                           ggcGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGtat                                                                    >  1:2801058/1‑148 (MQ=255)
                                                                               cAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACt                                                                <  1:759105/148‑1 (MQ=255)
                                                                               cAACGGCAGATTGCTGCCGAGAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACt                                                                <  2:2279797/148‑1 (MQ=255)
                                                                                aaCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                                               <  1:971323/148‑1 (MQ=255)
                                                                                aaCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                                               <  2:3187665/148‑1 (MQ=255)
                                                                                aaCGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACtt                                                               <  1:2496060/148‑1 (MQ=255)
                                                                                    gCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgcg                                                           >  2:3440792/1‑148 (MQ=255)
                                                                                    gCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAgcg                                                           >  2:1672591/1‑148 (MQ=255)
                                                                                        aTTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAc                                                       >  2:2295389/1‑148 (MQ=255)
                                                                                         ttGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAcc                                                      >  1:1274420/1‑148 (MQ=255)
                                                                                           gctgcCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAAccgc                                                    <  2:1367447/148‑1 (MQ=255)
                                                                                                   tAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                                            >  2:2325615/1‑148 (MQ=255)
                                                                                                   tAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGtt                                            >  1:2884266/1‑148 (MQ=255)
                                                                                                       cTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGcc                                        <  2:3200143/148‑1 (MQ=255)
                                                                                                          aCGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAg                                     >  2:1078466/1‑148 (MQ=255)
                                                                                                             cAACTGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCgg                                  >  1:1356599/1‑148 (MQ=255)
                                                                                                                cTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTc                               <  1:861145/148‑1 (MQ=255)
                                                                                                                 tGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTcc                              <  1:2111191/148‑1 (MQ=255)
                                                                                                                 tGACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTcc                              <  2:3027361/148‑1 (MQ=255)
                                                                                                                  gACCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTccc                             >  1:1930880/1‑148 (MQ=255)
                                                                                                                   aCCGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCg                            >  2:2437009/1‑148 (MQ=255)
                                                                                                                    ccGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGt                           <  1:3087029/148‑1 (MQ=255)
                                                                                                                    ccGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGt                           <  1:3087166/148‑1 (MQ=255)
                                                                                                                    ccGCGTTTCGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGt                           >  1:2100146/1‑148 (MQ=255)
                                                                                                                            cGGCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGtt                   <  2:2801058/148‑1 (MQ=255)
                                                                                                                              gCATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTg                 <  2:609552/148‑1 (MQ=255)
                                                                                                                               cATACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGt                >  2:150913/1‑148 (MQ=255)
                                                                                                                                aTACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGtt               <  1:2841565/148‑1 (MQ=255)
                                                                                                                                 tACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTa              >  2:1642975/1‑148 (MQ=255)
                                                                                                                                 tACCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTa              >  1:671956/1‑148 (MQ=255)
                                                                                                                                  aCCTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAc             >  1:1378757/1‑148 (MQ=255)
                                                                                                                                   ccTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAca            >  1:1420496/1‑148 (MQ=255)
                                                                                                                                    cTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAcag           >  1:3461595/1‑148 (MQ=255)
                                                                                                                                    cTTCAGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTAcag           >  1:3461563/1‑148 (MQ=255)
                                                                                                                                        aGCCATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGt       <  2:464629/148‑1 (MQ=255)
                                                                                                                                           cATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTAc    >  1:246239/1‑148 (MQ=255)
                                                                                                                                           cATTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTAc    >  1:3439884/1‑148 (MQ=255)
                                                                                                                                            aTTTTCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACg   >  2:1064933/1‑148 (MQ=255)
                                                                                                                                             ttttCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGc  <  2:905139/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TCCGCCGAGCGGATTGTGGGGCGGTTTATACTGTTTCCCGCAGTTTGCCGACGAAGAAAGTTTGCGGCAGTGGCTGGCGCAACGGCAGATTGCTGCCGATAACCTGACGCAACTGACCGCGTTTCGGCATACCTTCAGCCATTTCCACTTAGATATTGTGCCTATGTGGCTTCCCGTGTCGTCATTCACCGGCTGCATGGATGAAGGCAATGCGCTCTGGTATAACTTAGCGCAACCGCCGTCAGTTGGCCTAGCGGCTCCCGTGGAGCGTTTGTTACAGCAGTTACGC  >  CP009273/3097118‑3097406

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: