Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,219,790 A→T 93.0% L94Q (CTG→CAG)  pepE ← (alpha)‑aspartyl dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,219,7900AT93.0% 900.2 / 41.5 286L94Q (CTG→CAG) pepE(alpha)‑aspartyl dipeptidase
Reads supporting (aligned to +/- strand):  ref base A (13/7);  new base T (116/150);  total (129/157)
Fisher's exact test for biased strand distribution p-value = 1.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 286 total aligned reads displayed.

CTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAG  >  CP009273/4219644‑4219936
                                                                                                                                                  |                                                                                                                                                  
  tGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAAc                                                                                                                                               >  2:1281746/1‑150 (MQ=255)
    cGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTg                                                                                                                                             <  2:4374028/150‑1 (MQ=255)
     ggATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTgg                                                                                                                                            <  1:811417/150‑1 (MQ=255)
     ggATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTgg                                                                                                                                            >  1:4650575/1‑150 (MQ=255)
      gATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGa                                                                                                                                           >  1:3340128/1‑150 (MQ=255)
       aTCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGaa                                                                                                                                          <  1:4165075/150‑1 (MQ=255)
        tCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGaaa                                                                                                                                         <  2:4710945/150‑1 (MQ=255)
          gACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATAGACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGt                                                                                                                                       <  2:2740805/150‑1 (MQ=255)
           aCAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTa                                                                                                                                      <  1:769450/150‑1 (MQ=255)
            cAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTAt                                                                                                                                     >  1:1567443/1‑150 (MQ=255)
             aaTCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTAtt                                                                                                                                    <  1:4500509/150‑1 (MQ=255)
                cGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTccc                                                                                                                                 >  2:4311847/1‑150 (MQ=255)
                         gTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGAc                                                                                                                        >  2:4985951/1‑150 (MQ=255)
                         gTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGAc                                                                                                                        >  1:3204800/1‑150 (MQ=255)
                           tggtggTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACaa                                                                                                                      >  2:2054603/1‑150 (MQ=255)
                             gtggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATc                                                                                                                    <  1:1304149/150‑1 (MQ=255)
                             gtggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATc                                                                                                                    <  1:2278717/150‑1 (MQ=255)
                               ggtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCAc                                                                                                                  <  1:1385447/150‑1 (MQ=255)
                                gtACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACg                                                                                                                 <  1:2601798/150‑1 (MQ=255)
                                 tACGAATTGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGa                                                                                                                <  1:26523/150‑1 (MQ=255)
                                 tACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGa                                                                                                                >  1:3933720/1‑150 (MQ=255)
                                  aCGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAt                                                                                                               >  1:2087602/1‑150 (MQ=255)
                                    gAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtct                                                                                                             >  1:2340000/1‑150 (MQ=255)
                                     aaTAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGAtctc                                                                                                            >  1:902290/1‑150 (MQ=255)
                                      aTAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCa                                                                                                           <  1:194723/150‑1 (MQ=255)
                                      aTAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCa                                                                                                           <  2:2422219/150‑1 (MQ=255)
                                           tGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtt                                                                                                      <  2:1410879/150‑1 (MQ=255)
                                            gggCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAttt                                                                                                     <  1:1822838/150‑1 (MQ=255)
                                             ggCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCAtttt                                                                                                    <  1:2457007/150‑1 (MQ=255)
                                              gcaagcaaGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTc                                                                                                   <  1:1087108/150‑1 (MQ=255)
                                               caagcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCa                                                                                                  >  2:1689593/1‑150 (MQ=255)
                                                aagcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCaa                                                                                                 >  1:4683631/1‑150 (MQ=255)
                                                aagcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCaa                                                                                                 >  2:117198/1‑150 (MQ=255)
                                                 agcaagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAAt                                                                                                >  1:1685114/1‑150 (MQ=255)
                                                    aagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACACTCACGATCTCAGCATTTTCAATcgc                                                                                             >  2:2251803/1‑150 (MQ=255)
                                                    aagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgc                                                                                             <  1:348785/150‑1 (MQ=255)
                                                    aagGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgc                                                                                             >  2:5152484/1‑150 (MQ=255)
                                                     agGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGTCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgcg                                                                                            >  2:818228/1‑150 (MQ=255)
                                                     agGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATcgcg                                                                                            >  1:632716/1‑150 (MQ=255)
                                                       gTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGc                                                                                          >  2:3232515/1‑150 (MQ=255)
                                                         tAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCaa                                                                                        >  2:2073167/1‑150 (MQ=255)
                                                         tAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCaa                                                                                        >  2:1894302/1‑150 (MQ=255)
                                                         tAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCaa                                                                                        <  1:859393/150‑1 (MQ=255)
                                                          aGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAAc                                                                                       >  2:2658696/1‑150 (MQ=255)
                                                          aGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAAc                                                                                       >  2:398673/1‑150 (MQ=255)
                                                             gcCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACggg                                                                                    >  2:1578026/1‑150 (MQ=255)
                                                             gcCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACggg                                                                                    <  1:3588952/150‑1 (MQ=255)
                                                               cTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGAt                                                                                  <  2:3575708/150‑1 (MQ=255)
                                                               cTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGAt                                                                                  >  1:3584393/1‑150 (MQ=255)
                                                                tGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATc                                                                                 <  1:5047704/150‑1 (MQ=255)
                                                                  cgcTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATcta                                                                               >  1:1151266/1‑148 (MQ=255)
                                                                   gcTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTAc                                                                              <  1:4633287/150‑1 (MQ=255)
                                                                    cTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGAATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACa                                                                             >  2:4067805/1‑150 (MQ=255)
                                                                     tCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACaa                                                                            >  1:688042/1‑150 (MQ=255)
                                                                      ccAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAAc                                                                           <  1:107221/150‑1 (MQ=255)
                                                                         gCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGca                                                                        <  2:1759867/150‑2 (MQ=255)
                                                                            aaTATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATg                                                                     <  2:4972376/150‑1 (MQ=255)
                                                                             atatACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGa                                                                    <  2:3059790/150‑1 (MQ=255)
                                                                              tataCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGaa                                                                   >  1:3549774/1‑150 (MQ=255)
                                                                              tataCAGAGCGCCACGTTTAACCAAGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAAGCGCGGCACCGGGATCTACAACGCTATGaa                                                                   >  2:4945672/1‑150 (MQ=255)
                                                                                taCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATa                                                                 >  2:3748090/1‑150 (MQ=255)
                                                                                taCAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATa                                                                 <  2:2447421/150‑1 (MQ=255)
                                                                                   agagCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCg                                                              <  1:2703451/150‑1 (MQ=255)
                                                                                    gagCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCgg                                                             <  2:2820669/150‑1 (MQ=255)
                                                                                        gcCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGAc                                                         >  1:4090637/1‑150 (MQ=255)
                                                                                         ccACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACa                                                        >  2:4467557/1‑150 (MQ=255)
                                                                                         ccACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACa                                                        <  1:3760549/150‑1 (MQ=255)
                                                                                         ccACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACa                                                        >  1:2521251/1‑150 (MQ=255)
                                                                                              tttAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCACCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAac                                                   >  2:104601/1‑150 (MQ=255)
                                                                                                tAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAacac                                                 >  1:3002106/1‑150 (MQ=255)
                                                                                                  aCCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACAccc                                               <  1:4248703/150‑1 (MQ=255)
                                                                                                    cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAg                                             <  2:2162242/150‑1 (MQ=255)
                                                                                                    cACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAg                                             >  1:1441198/1‑150 (MQ=255)
                                                                                                       gTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCgg                                          <  1:4829413/150‑1 (MQ=255)
                                                                                                        tCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGa                                         <  1:5131411/150‑1 (MQ=255)
                                                                                                         cAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAg                                        <  2:4649772/150‑1 (MQ=255)
                                                                                                           gTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCg                                      <  2:2284020/150‑1 (MQ=255)
                                                                                                           gTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCg                                      <  2:4145301/150‑1 (MQ=255)
                                                                                                             aaTTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCgag                                    >  2:292296/1‑150 (MQ=255)
                                                                                                                 tgTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAAcc                                >  1:1924085/2‑150 (MQ=255)
                                                                                                                    gCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCa                             <  1:4961339/150‑1 (MQ=255)
                                                                                                                    gCCAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCa                             <  2:3275465/150‑1 (MQ=255)
                                                                                                                     ccAGCAGCCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAg                            <  1:3064343/150‑1 (MQ=255)
                                                                                                                       agcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                          <  1:2054603/150‑1 (MQ=255)
                                                                                                                       agcagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                          <  1:4311847/150‑1 (MQ=255)
                                                                                                                       agcagcCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGcc                          <  2:930978/150‑1 (MQ=255)
                                                                                                                         cagcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGt                        >  1:4736553/1‑150 (MQ=255)
                                                                                                                          agcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                       <  2:1488047/150‑1 (MQ=255)
                                                                                                                          agcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                       >  1:3662503/1‑150 (MQ=255)
                                                                                                                          agcCCGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtt                       >  2:1059345/1‑150 (MQ=255)
                                                                                                                            cccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtttt                     <  2:5084079/150‑1 (MQ=255)
                                                                                                                            cccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGtttt                     <  1:1245047/150‑1 (MQ=255)
                                                                                                                            cccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCCCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGACCCGCAGCCGtttt                     >  2:2979450/1‑150 (MQ=255)
                                                                                                                             ccGCGCTCGCGGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTc                    <  1:585176/150‑1 (MQ=255)
                                                                                                                                     gcgGCACTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTa            >  2:1832575/1‑150 (MQ=255)
                                                                                                                                       ggCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAat          >  1:4287476/1‑150 (MQ=255)
                                                                                                                                          aCTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACGGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcat       <  2:2865053/150‑1 (MQ=255)
                                                                                                                                           cTGTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAatcatc      >  2:738683/1‑150 (MQ=255)
                                                                                                                                             gTTTCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATccc    >  1:2796798/1‑150 (MQ=255)
                                                                                                                                               ttCTGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCTACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAg  <  1:2898570/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTTGCGGATCGACAATCGGCATATCGTTGGTGGTACGAATAGTTGGGCAAGCAAGGTTAGCGCCTGCGCTCCAGCCAATATACAGAGCGCCACGTTTAACCACGTCAGTAATTGGTGCCAGCAGCCCGCGCTCGCGGCACTGTTTCAGCAACTGGAAAGTATTCCCGCCGCCGACAATCACGATCTCAGCATTTTCAATCGCGGCAACGGGATCGACAACGCTATGAATACCGGTGACAGAAACACCCAGCGGAGCGAGAACCGCAGCCGTTTTCGCTGTGTAATCATCCCAG  >  CP009273/4219644‑4219936

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: