Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,386,672 C→A 12.6% P982P (CCG→CCT yhdP ← DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,386,6720CA12.6% 1058.5 / 84.2 318P982P (CCG→CCTyhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base C (130/148);  new base A (22/18);  total (152/166)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 319 total aligned reads displayed.

GTATCAATCAGTCCATTGGTCAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCAACCAGGGTGC  >  CP009273/3386523‑3386820
                                                                                                                                                     |                                                                                                                                                    
   tcaatcaGTCCATTGGTCAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCggg                                                                                                                                                   <  2:2302713/150‑1 (MQ=255)
         aGTCCATTGGTCAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGt                                                                                                                                             >  1:4166957/1‑150 (MQ=255)
               ttGGTCAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGAc                                                                                                                                       <  1:2587384/150‑1 (MQ=255)
               ttGGTCAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGAc                                                                                                                                       <  1:2586319/150‑1 (MQ=255)
                 ggTCAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACgg                                                                                                                                     <  2:1385319/150‑1 (MQ=255)
                    cAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGt                                                                                                                                  <  2:1433041/150‑1 (MQ=255)
                     aGCGTTAACGTATCGCCAGACATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGtt                                                                                                                                 <  1:3934333/150‑1 (MQ=255)
                     aGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGtt                                                                                                                                 >  1:2252166/1‑150 (MQ=255)
                       cGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGa                                                                                                                               >  2:604940/1‑150 (MQ=255)
                        gTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGaa                                                                                                                              >  2:1676853/1‑150 (MQ=255)
                         ttAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAAt                                                                                                                             <  2:2882045/150‑1 (MQ=255)
                          tAACGTATCGCCAGAAATGGTGATATCACTTTCATTGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATc                                                                                                                            <  1:2506925/150‑1 (MQ=255)
                          tAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATc                                                                                                                            >  2:221646/1‑150 (MQ=255)
                          tAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATc                                                                                                                            >  2:2044444/1‑150 (MQ=255)
                           aaCGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATca                                                                                                                           <  1:606893/150‑1 (MQ=255)
                             cGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATcacc                                                                                                                         <  1:521469/150‑1 (MQ=255)
                              gTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATcacca                                                                                                                        >  1:2078558/1‑150 (MQ=255)
                                 tCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGa                                                                                                                     >  2:1433963/1‑150 (MQ=255)
                                 tCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGa                                                                                                                     >  1:2783837/1‑150 (MQ=255)
                                  cGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAg                                                                                                                    <  1:604486/150‑1 (MQ=255)
                                    ccAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGtt                                                                                                                  >  2:1596485/1‑150 (MQ=255)
                                    ccAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGtt                                                                                                                  <  1:4359275/150‑1 (MQ=255)
                                     cAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGttt                                                                                                                 <  1:1255289/150‑1 (MQ=255)
                                     cAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGttt                                                                                                                 <  1:3528942/150‑1 (MQ=255)
                                     cAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGttt                                                                                                                 >  1:4280115/1‑150 (MQ=255)
                                     cAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGCTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGttt                                                                                                                 >  1:637158/1‑150 (MQ=255)
                                       gAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTc                                                                                                               >  1:2859451/1‑150 (MQ=255)
                                       gAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTc                                                                                                               <  2:2101937/150‑1 (MQ=255)
                                       gAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTc                                                                                                               <  2:428719/150‑1 (MQ=255)
                                       gAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGACAGGTGATGACGGCGTTGAATCACCACGAGTTTTc                                                                                                               >  1:2952633/1‑150 (MQ=255)
                                         aaTGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGc                                                                                                             >  1:4069429/1‑150 (MQ=255)
                                          aTGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGcc                                                                                                            >  2:468078/1‑150 (MQ=255)
                                          aTGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGcc                                                                                                            >  2:2752696/1‑150 (MQ=255)
                                           tGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCa                                                                                                           >  2:3335061/1‑150 (MQ=255)
                                            ggTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCAc                                                                                                          <  1:578272/150‑1 (MQ=255)
                                               gATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACTGCGTTGAATCACCACGAGTTTTCGCCACGCt                                                                                                       <  2:2584746/150‑1 (MQ=255)
                                                   tCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCggg                                                                                                   <  1:40039/150‑1 (MQ=255)
                                                    cACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGt                                                                                                  >  2:638651/1‑150 (MQ=255)
                                                    cACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGt                                                                                                  >  1:1489397/1‑150 (MQ=255)
                                                    cACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGt                                                                                                  >  1:2599365/1‑150 (MQ=255)
                                                     aCTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGtt                                                                                                 <  1:271909/150‑1 (MQ=255)
                                                      cTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGACTCACCACGAGTTTTCGCCACGCTCGGGTTg                                                                                                >  2:3870367/1‑150 (MQ=255)
                                                      cTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTg                                                                                                >  1:1683574/1‑150 (MQ=255)
                                                       tGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGTTGACGGCGTTGAATCACCACGCGTTTTCGCCACGCTCGGGTTGt                                                                                               <  2:3152877/150‑1 (MQ=255)
                                                       tGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGt                                                                                               >  1:3828114/1‑150 (MQ=255)
                                                       tGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGt                                                                                               >  2:3123884/1‑150 (MQ=255)
                                                       tGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGt                                                                                               >  1:221308/1‑150 (MQ=255)
                                                        gTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGta                                                                                              >  1:116151/1‑150 (MQ=255)
                                                          cAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGtaat                                                                                            <  1:745455/150‑1 (MQ=255)
                                                          cAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGtaat                                                                                            <  2:3462638/150‑1 (MQ=255)
                                                            aTGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGtaataa                                                                                          >  1:3567103/1‑150 (MQ=255)
                                                             tGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATaaa                                                                                         >  2:4539557/1‑150 (MQ=255)
                                                                gACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGat                                                                                      >  2:4345506/1‑150 (MQ=255)
                                                                 aCCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGata                                                                                     <  1:268428/150‑1 (MQ=255)
                                                                 aCCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGata                                                                                     >  2:1205361/1‑150 (MQ=255)
                                                                 aCCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTAGTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCGCGGGTTGTAATAAAGata                                                                                     >  1:3838255/1‑150 (MQ=255)
                                                                     aaTTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTg                                                                                 <  1:3028928/150‑1 (MQ=255)
                                                                     aaTTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTg                                                                                 <  2:1892860/150‑1 (MQ=255)
                                                                     aaTTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTg                                                                                 <  2:1037614/150‑1 (MQ=255)
                                                                     aaTTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTg                                                                                 <  1:2439418/150‑1 (MQ=255)
                                                                     aaTTTTTGACCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTg                                                                                 <  1:1449608/150‑1 (MQ=255)
                                                                      aTTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGa                                                                                >  1:1901845/1‑150 (MQ=255)
                                                                           tGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTc                                                                           >  1:3464949/1‑150 (MQ=255)
                                                                              ccccAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCgcc                                                                        >  1:1393218/1‑150 (MQ=255)
                                                                              ccccAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCgcc                                                                        >  1:3660079/1‑150 (MQ=255)
                                                                               cccAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCgcca                                                                       <  1:907456/150‑1 (MQ=255)
                                                                                ccAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCgccag                                                                      >  2:3557952/1‑150 (MQ=255)
                                                                                ccAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCgccag                                                                      >  2:3394388/1‑150 (MQ=255)
                                                                                ccAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCgccag                                                                      <  2:369930/150‑1 (MQ=255)
                                                                                  aGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCgccagcc                                                                    >  2:2315805/1‑150 (MQ=255)
                                                                                         gCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGgcgc                                                             >  1:1665718/1‑150 (MQ=255)
                                                                                          cACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGcgcg                                                            <  1:241551/150‑1 (MQ=255)
                                                                                            ctctGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGtt                                                          >  1:85659/1‑150 (MQ=255)
                                                                                                  gCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACg                                                    >  2:1175430/1‑150 (MQ=255)
                                                                                                      cGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATc                                                >  1:4398758/1‑150 (MQ=255)
                                                                                                      cGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATc                                                >  2:1586527/1‑150 (MQ=255)
                                                                                                       gTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCg                                               >  2:1337221/1‑150 (MQ=255)
                                                                                                        tctttGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGc                                              <  2:3509940/148‑1 (MQ=255)
                                                                                                         aTTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCt                                             >  2:296695/1‑150 (MQ=255)
                                                                                                            tGGGCGTCCGGCCAGCCACGGAAGTTAATTCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGc                                          <  2:886157/150‑1 (MQ=255)
                                                                                                            tGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGc                                          >  1:1845889/1‑150 (MQ=255)
                                                                                                                cGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTgg                                      <  1:3148871/150‑1 (MQ=255)
                                                                                                                 gTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGc                                     <  2:94662/150‑1 (MQ=255)
                                                                                                                 gTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGc                                     >  2:769245/1‑150 (MQ=255)
                                                                                                                     ggCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGTTATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTg                                 <  2:3201651/150‑1 (MQ=255)
                                                                                                                     ggCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTg                                 >  2:1734564/1‑150 (MQ=255)
                                                                                                                     ggCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTg                                 >  1:813622/1‑150 (MQ=255)
                                                                                                                      gccagccaCGGAAGTTAATGCGCTCCGTTGTAGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGa                                <  1:4019004/150‑1 (MQ=255)
                                                                                                                          gccaCGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTc                            >  1:2090948/1‑150 (MQ=255)
                                                                                                                            caCGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCAc                          >  2:427423/1‑150 (MQ=255)
                                                                                                                            caCGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCAc                          >  1:2969225/1‑150 (MQ=255)
                                                                                                                            caCGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCAc                          <  1:1642071/150‑1 (MQ=255)
                                                                                                                                 aaGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCt                     <  1:493546/150‑1 (MQ=255)
                                                                                                                                        tGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCtt              <  1:4105299/150‑1 (MQ=255)
                                                                                                                                          cgcTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCa            <  1:3064918/150‑1 (MQ=255)
                                                                                                                                          cgcTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCa            <  1:3064610/150‑1 (MQ=255)
                                                                                                                                          cgcTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCa            <  2:3784515/150‑1 (MQ=255)
                                                                                                                                             tCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCAAcc         <  2:4052610/150‑1 (MQ=255)
                                                                                                                                                gTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCAACCAgg      <  1:1564369/150‑1 (MQ=255)
                                                                                                                                                  tGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCAACCAGGGt    >  1:2369151/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTATCAATCAGTCCATTGGTCAGCGTTAACGTATCGCCAGAAATGGTGATATCACTGTCAATGCGACCGAATTTTTGCCCCCAGAACCAGCACTCTGTGCATCGTATTTGGGCGTCCGGCCAGCCACGGAAGTTAATGCGCTCCGTTGTCGGGAAAGGTGATGACGGCGTTGAATCACCACGAGTTTTCGCCACGCTCGGGTTGTAATAAAGATATTTGATATTCGCCAGCCACGGCGCGTTATTACGCATCGCTAGCGTGGCGTTGATTTCACGCCCTTGCGCTTCAACCAGGGTGC  >  CP009273/3386523‑3386820

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: