Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 510,357 A→T 9.2% I291I (ATT→ATA qmcA ← PHB domain membrane‑anchored putative protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273510,3570AT 9.2% 1236.6 / 88.9 358I291I (ATT→ATAqmcAPHB domain membrane‑anchored putative protease
Reads supporting (aligned to +/- strand):  ref base A (160/165);  new base T (18/15);  total (178/180)
Fisher's exact test for biased strand distribution p-value = 5.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 366 total aligned reads displayed.

CCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCGATAGC  >  CP009273/510209‑510506
                                                                                                                                                    |                                                                                                                                                     
 cgcTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAAt                                                                                                                                                     >  1:3781686/1‑150 (MQ=255)
   cTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATcc                                                                                                                                                   >  1:3288784/1‑150 (MQ=255)
   cTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATcc                                                                                                                                                   >  1:2975472/1‑150 (MQ=255)
    tCCACAACAAATAACCAGTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATccc                                                                                                                                                  <  2:4676691/150‑1 (MQ=255)
     ccACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCg                                                                                                                                                 <  2:3274810/150‑1 (MQ=255)
      cacaACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCgg                                                                                                                                                <  1:3160923/150‑1 (MQ=255)
           caaATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATc                                                                                                                                           >  2:854121/1‑150 (MQ=255)
            aaaTAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCg                                                                                                                                          <  2:1105206/150‑1 (MQ=255)
             aaTAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGa                                                                                                                                         >  1:949588/1‑150 (MQ=255)
               tAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGAcc                                                                                                                                       <  1:110577/150‑1 (MQ=255)
               tAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGAcc                                                                                                                                       >  2:3875862/1‑150 (MQ=255)
                 aCCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGAcccc                                                                                                                                     >  2:3803225/1‑150 (MQ=255)
                  ccATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCa                                                                                                                                    <  1:4503522/150‑1 (MQ=255)
                    aTTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATc                                                                                                                                  >  1:3502453/1‑150 (MQ=255)
                    aTTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATc                                                                                                                                  >  1:638471/1‑150 (MQ=255)
                      ttCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAg                                                                                                                                >  2:2234655/1‑150 (MQ=255)
                         cGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGgct                                                                                                                             <  2:3933349/150‑1 (MQ=255)
                         cGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGgct                                                                                                                             >  2:32317/1‑150 (MQ=255)
                         cGCCCAGAATCTCGGGTGCCTGATGCAACCCGCCGTGACTGAGCCATAAATTTTGTGGTTGAACGACCATTAACTCCATCTTGGCTGAGTCGGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGgct                                                                                                                             <  1:18800/150‑1 (MQ=255)
                            ccAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGgctgct                                                                                                                          <  1:156078/150‑1 (MQ=255)
                            ccAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGgctgct                                                                                                                          <  1:1652432/150‑1 (MQ=255)
                             cAGCATCTCGGCGGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAATTATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGctgctg                                                                                                                         <  1:3248532/150‑1 (MQ=255)
                              aGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTgg                                                                                                                        <  1:3061861/150‑1 (MQ=255)
                                cATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGcc                                                                                                                      >  2:4689148/1‑150 (MQ=255)
                                 aTCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCTATCCCGGCAATCGACCCCATCAGGCTGCTGGCct                                                                                                                     >  2:3634835/1‑150 (MQ=255)
                                 aTCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCct                                                                                                                     >  1:4065080/1‑150 (MQ=255)
                                 aTCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCct                                                                                                                     >  1:1678668/1‑150 (MQ=255)
                                 aTCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCct                                                                                                                     <  2:3288668/150‑1 (MQ=255)
                                  tctcGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCctc                                                                                                                    >  1:142118/1‑150 (MQ=255)
                                   ctcGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCctct                                                                                                                   >  2:3344624/1‑150 (MQ=255)
                                     cGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTaa                                                                                                                 >  2:1980479/1‑150 (MQ=255)
                                           ccAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCCGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGcat                                                                                                           <  1:816916/150‑1 (MQ=255)
                                             agcagcAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGcatca                                                                                                         >  1:1933386/1‑150 (MQ=255)
                                              gcagcaAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCTATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGcatcat                                                                                                        <  2:1205128/150‑1 (MQ=255)
                                                agcaAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCAtcatca                                                                                                      >  1:3996301/1‑150 (MQ=255)
                                                  caAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCact                                                                                                    <  2:141071/150‑1 (MQ=255)
                                                    aaCCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCactac                                                                                                  >  2:736051/1‑150 (MQ=255)
                                                         ccgAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGc                                                                                             >  2:1462631/1‑150 (MQ=255)
                                                          cgAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCt                                                                                            >  1:4457543/1‑150 (MQ=255)
                                                          cgAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCt                                                                                            >  1:4622853/1‑150 (MQ=255)
                                                           gagaCTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTg                                                                                           >  2:4314965/1‑150 (MQ=255)
                                                           gagaCTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTg                                                                                           >  2:254308/1‑150 (MQ=255)
                                                            agaCTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGt                                                                                          >  2:4409966/1‑150 (MQ=255)
                                                              aCTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTa                                                                                        >  1:2404047/1‑150 (MQ=255)
                                                                tGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCTATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACt                                                                                      <  2:3555644/150‑1 (MQ=255)
                                                                tGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACt                                                                                      <  2:4570203/150‑1 (MQ=255)
                                                                  aGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTgg                                                                                    >  2:2340031/1‑150 (MQ=255)
                                                                    ccAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTggag                                                                                  <  2:2193294/150‑1 (MQ=255)
                                                                    ccAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTggag                                                                                  <  2:3157629/150‑1 (MQ=255)
                                                                     cAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTggagg                                                                                 >  2:2591903/1‑150 (MQ=255)
                                                                      aGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTggagga                                                                                >  2:903416/1‑150 (MQ=255)
                                                                      aGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTggagga                                                                                >  1:272866/1‑150 (MQ=255)
                                                                         aaaTATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAAcc                                                                             <  1:884499/150‑1 (MQ=255)
                                                                          aaTATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCg                                                                            >  1:3330919/1‑150 (MQ=255)
                                                                             atGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATc                                                                         <  2:4059694/150‑1 (MQ=255)
                                                                             atGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATc                                                                         <  2:1669977/150‑1 (MQ=255)
                                                                               gtgGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATctg                                                                       <  2:4232980/150‑1 (MQ=255)
                                                                                 ggATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATctgct                                                                     <  2:640364/150‑1 (MQ=255)
                                                                                   aTGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGt                                                                   <  2:588831/150‑1 (MQ=255)
                                                                                     gAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTaa                                                                 <  1:2628348/150‑1 (MQ=255)
                                                                                          aCCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCtt                                                            <  1:4406610/150‑1 (MQ=255)
                                                                                           ccATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTc                                                           <  1:1594701/150‑1 (MQ=255)
                                                                                            cATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCg                                                          >  2:1527937/1‑150 (MQ=255)
                                                                                             aTTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCgg                                                         >  2:1696247/1‑150 (MQ=255)
                                                                                             aTTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCgg                                                         <  1:2065699/150‑1 (MQ=255)
                                                                                                      atcatGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCt                                                >  2:3476929/1‑150 (MQ=255)
                                                                                                       tcatGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTg                                               <  1:3891793/150‑1 (MQ=255)
                                                                                                        catGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTgc                                              <  1:2281581/150‑1 (MQ=255)
                                                                                                           ggCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCTATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCt                                           <  2:2102077/150‑1 (MQ=255)
                                                                                                           ggCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCt                                           <  2:300371/150‑1 (MQ=255)
                                                                                                               gAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGa                                       >  2:1848798/1‑150 (MQ=255)
                                                                                                               gAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGa                                       <  1:2155825/150‑1 (MQ=255)
                                                                                                                 gTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAag                                     >  2:3192618/1‑150 (MQ=255)
                                                                                                                   ccGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAagta                                   >  2:4279801/1‑150 (MQ=255)
                                                                                                                    cGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAagtag                                  >  1:1351780/1‑150 (MQ=255)
                                                                                                                     gTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAagtagt                                 <  1:4561957/150‑1 (MQ=255)
                                                                                                                       ttgttgGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTc                               <  2:791242/150‑1 (MQ=255)
                                                                                                                       ttgttgGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTc                               <  1:4651317/150‑1 (MQ=255)
                                                                                                                               gcTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTg                       <  1:579050/150‑1 (MQ=255)
                                                                                                                                 tGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGaa                     >  1:1776734/1‑150 (MQ=255)
                                                                                                                                  gTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAat                    <  1:52407/150‑1 (MQ=255)
                                                                                                                                   tCTTTCACCAGCTCGGCTATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAata                   >  1:131419/1‑150 (MQ=255)
                                                                                                                                    cTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAatat                  <  2:2946986/150‑1 (MQ=255)
                                                                                                                                     tttCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATc                 >  1:2585474/1‑150 (MQ=255)
                                                                                                                                     tttCACCAGCTCGGAAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATc                 >  1:2911414/1‑150 (MQ=255)
                                                                                                                                      ttCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCa                >  2:2743609/1‑150 (MQ=255)
                                                                                                                                      ttCACCAGCTCGGCAATCCCGGCAATAGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCa                <  1:1008942/150‑1 (MQ=255)
                                                                                                                                       tCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCAc               <  2:151116/150‑1 (MQ=255)
                                                                                                                                        cACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCAcc              >  2:3207548/1‑150 (MQ=255)
                                                                                                                                           cAGCTCGGCTATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGa           >  1:217158/1‑150 (MQ=255)
                                                                                                                                           cAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGa           >  2:4255892/1‑150 (MQ=255)
                                                                                                                                            aGCTCGGCTATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAg          <  1:3634835/150‑1 (MQ=255)
                                                                                                                                            aGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAg          <  1:1503325/150‑1 (MQ=255)
                                                                                                                                             gCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTGTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAgg         <  2:179578/150‑1 (MQ=255)
                                                                                                                                               tCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCg       >  2:100744/1‑150 (MQ=255)
                                                                                                                                                cGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCGa      >  1:1728164/1‑150 (MQ=255)
                                                                                                                                                cGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCGa      >  1:4116175/1‑150 (MQ=255)
                                                                                                                                                    aaTCCCGGCAATCGACCCCATCAGGCTGCTGGCGTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCGATAGc  >  2:2809266/1‑150 (MQ=255)
                                                                                                                                                    aaTCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCGATAGc  >  1:3403001/1‑150 (MQ=255)
                                                                                                                                                    aaTCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCGATAGc  >  2:849188/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCGCTCCACAACAAATAACCATTTCCGCCCAGCATCTCGGCTGCCAGCAGCAAACCGCCGAGACTGAGCCAGAAAATATGTGGATGAACGACCATTAACTCCATCATGGCTGAGTCCGTTTGTTGGCGCTGTCTTTCACCAGCTCGGCAATCCCGGCAATCGACCCCATCAGGCTGCTGGCCTCTAATGGCATCATCACTACTTTGCTGTTACTGGAGGAACCGATCTGCTGTAACGCTTCGGTGTATTTCTGCGCTACGAAGTAGTTCACCGCCTGAATATCACCGGAGGCGATAGC  >  CP009273/510209‑510506

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: