Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 170,868 C→A 7.3% V168F (GTT→TTT)  hemL ← glutamate‑1‑semialdehyde aminotransferase (aminomutase)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273170,8680CA 7.3% 1082.9 / 40.9 293V168F (GTT→TTT) hemLglutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base C (126/145);  new base A (9/13);  total (135/158)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 297 total aligned reads displayed.

GGCAGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTCG  >  CP009273/170720‑171017
                                                                                                                                                    |                                                                                                                                                     
   aGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGccc                                                                                                                                                   >  1:1302739/1‑150 (MQ=255)
     cggcggAACACAGTTCATATTGCCTGCCACCGGCTCGTCGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCgg                                                                                                                                                 <  2:3171022/150‑1 (MQ=255)
      ggcggAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGc                                                                                                                                                <  2:3035952/150‑1 (MQ=255)
          gAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGt                                                                                                                                            >  1:3229738/1‑150 (MQ=255)
           aaCACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGtt                                                                                                                                           >  2:1513622/1‑150 (MQ=255)
                     atatTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCt                                                                                                                                 <  2:4063552/150‑1 (MQ=255)
                        ttGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAAc                                                                                                                              <  2:856942/150‑1 (MQ=255)
                               ccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGTTCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTg                                                                                                                       >  2:1560043/1‑150 (MQ=255)
                               ccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTg                                                                                                                       >  2:3846771/1‑150 (MQ=255)
                                  ccGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTgcgc                                                                                                                    <  2:696459/150‑1 (MQ=255)
                                  ccGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTgcgc                                                                                                                    >  1:1573962/1‑150 (MQ=255)
                                   cGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGcc                                                                                                                   <  1:4623913/150‑1 (MQ=255)
                                    ggCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCa                                                                                                                  <  1:860093/150‑1 (MQ=255)
                                    ggCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCa                                                                                                                  <  2:1644224/150‑1 (MQ=255)
                                      cTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGa                                                                                                                >  2:1941017/1‑150 (MQ=255)
                                         gacgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAAcc                                                                                                             <  1:2126528/150‑1 (MQ=255)
                                            gaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGc                                                                                                          >  1:4699350/1‑150 (MQ=255)
                                             ataataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCt                                                                                                         <  1:5115335/150‑1 (MQ=255)
                                             ataataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCt                                                                                                         >  2:621014/1‑150 (MQ=255)
                                              taataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCtt                                                                                                        <  2:1195517/150‑1 (MQ=255)
                                               aataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCttt                                                                                                       <  2:4199128/150‑1 (MQ=255)
                                                ataCAGGCAATCTCTTGCGGTTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAATGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGACTAACGTGAGTGCGCCAGAACCGGCTTTc                                                                                                      <  1:4780831/150‑1 (MQ=255)
                                                ataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTc                                                                                                      >  2:2778851/1‑150 (MQ=255)
                                                 taCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCa                                                                                                     <  1:3674705/150‑1 (MQ=255)
                                                   cAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCAcc                                                                                                   <  2:4454256/150‑1 (MQ=255)
                                                    aGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACca                                                                                                  <  1:261821/150‑1 (MQ=255)
                                                       cAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCCGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGATGGCCTAACGTGAGTGCTCCAGAACCGGCTTTCACcagca                                                                                               >  1:21739/1‑150 (MQ=255)
                                                        aaTCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACcagcag                                                                                              <  1:4585899/150‑1 (MQ=255)
                                                           ctctTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGTGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAgcaggca                                                                                           <  1:3335165/150‑1 (MQ=255)
                                                           ctctTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAgcaggca                                                                                           <  1:3613290/150‑1 (MQ=255)
                                                            tctTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAgcaggcag                                                                                          <  2:240139/150‑1 (MQ=255)
                                                            tctTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAgcaggcag                                                                                          <  2:924643/150‑1 (MQ=255)
                                                             ctTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGt                                                                                         >  1:249584/1‑150 (MQ=255)
                                                               tGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCa                                                                                       <  2:1451395/150‑1 (MQ=255)
                                                               tGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCa                                                                                       <  2:1232504/150‑1 (MQ=255)
                                                                 cGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGc                                                                                     >  2:4295678/1‑150 (MQ=255)
                                                                    gTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTg                                                                                  >  1:1006022/1‑150 (MQ=255)
                                                                    gTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTg                                                                                  >  2:1035485/1‑150 (MQ=255)
                                                                      aTTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGAc                                                                                >  2:3397652/1‑150 (MQ=255)
                                                                        tgataaaaTGCTGCTCTTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGATGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCa                                                                              <  2:4127848/145‑1 (MQ=255)
                                                                         gCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCAt                                                                             <  1:1015653/150‑1 (MQ=255)
                                                                          cTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATg                                                                            <  1:422105/150‑1 (MQ=255)
                                                                             aaaTGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTa                                                                         <  2:4804785/150‑1 (MQ=255)
                                                                                tGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAaca                                                                      >  2:2120222/1‑150 (MQ=255)
                                                                                     cgcgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCTAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCtt                                                                 <  1:2548060/150‑1 (MQ=255)
                                                                                        gTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCaa                                                              >  1:912772/1‑150 (MQ=255)
                                                                                         tACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCaaa                                                             <  1:2989730/150‑1 (MQ=255)
                                                                                           cAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                           >  2:4397126/1‑150 (MQ=255)
                                                                                            aGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAttt                                                          >  2:274493/1‑150 (MQ=255)
                                                                                            aGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAttt                                                          <  2:2622524/150‑1 (MQ=255)
                                                                                             gAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCCGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTta                                                         >  2:914626/1‑150 (MQ=255)
                                                                                               aGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaat                                                       <  1:396591/150‑1 (MQ=255)
                                                                                               aGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaat                                                       >  2:4596955/1‑150 (MQ=255)
                                                                                                gCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaatc                                                      >  2:1050459/1‑149 (MQ=255)
                                                                                                  cAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTaataat                                                    >  1:2059865/1‑150 (MQ=255)
                                                                                                  cAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTaataat                                                    <  1:1156700/150‑1 (MQ=255)
                                                                                                    gATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAAttt                                                  >  1:4162544/1‑150 (MQ=255)
                                                                                                    gATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTAAATAAttt                                                  >  2:3644402/1‑150 (MQ=255)
                                                                                                      tCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTg                                                <  2:2520881/150‑1 (MQ=255)
                                                                                                      tCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTg                                                <  1:3461223/150‑1 (MQ=255)
                                                                                                       cATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGt                                               <  1:5186760/150‑1 (MQ=255)
                                                                                                        aTTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTc                                              >  2:3094973/1‑150 (MQ=255)
                                                                                                         ttATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTcg                                             >  1:4988024/1‑150 (MQ=255)
                                                                                                            taAGTACAGGTTAAGGTATAATTGGCGAAATCTGCCGGAACGCCCGGCGAGTGTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGa                                          <  2:3555270/150‑1 (MQ=255)
                                                                                                              aGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGAcc                                        <  1:831755/150‑1 (MQ=255)
                                                                                                               gTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCg                                       >  1:688148/1‑150 (MQ=255)
                                                                                                               gTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCg                                       >  1:862254/1‑150 (MQ=255)
                                                                                                                tACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCgg                                      >  2:294826/1‑150 (MQ=255)
                                                                                                                  cAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTa                                    >  2:3441504/1‑150 (MQ=255)
                                                                                                                  cAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTa                                    <  2:4011858/150‑1 (MQ=255)
                                                                                                                   aGGTTAAGGTATATTTTGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGAGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCTACCGGTaa                                   >  1:1945528/1‑150 (MQ=255)
                                                                                                                   aGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaa                                   <  1:4529420/150‑1 (MQ=255)
                                                                                                                   aGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaa                                   >  2:3818117/1‑150 (MQ=255)
                                                                                                                    ggTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaaa                                  <  2:4483945/150‑1 (MQ=255)
                                                                                                                        aaGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCa                              <  1:120510/150‑1 (MQ=255)
                                                                                                                        aaGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCa                              <  1:120652/150‑1 (MQ=255)
                                                                                                                         aGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCAc                             >  2:562105/1‑150 (MQ=255)
                                                                                                                         aGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCAc                             >  2:5027755/1‑150 (MQ=255)
                                                                                                                          ggTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACg                            <  2:328290/150‑1 (MQ=255)
                                                                                                                          ggTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCAACAGCAGGCAGTCAGCGTGACCATGGTAACTCCCTTCAAATTTATTTATTTTGTCGCGACCGGTAAACCCACg                            <  2:1386131/150‑1 (MQ=255)
                                                                                                                           gTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACgg                           >  1:4922048/1‑150 (MQ=255)
                                                                                                                             atatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGc                         >  1:598050/1‑150 (MQ=255)
                                                                                                                             atatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGc                         >  1:5064392/1‑150 (MQ=255)
                                                                                                                              tatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGcc                        <  2:4979655/150‑1 (MQ=255)
                                                                                                                              tatTTGGCGAAATCTGACGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACAGGATTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACAACGGGcc                        >  2:1744486/1‑150 (MQ=255)
                                                                                                                               ctttGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCa                       <  2:4844460/149‑1 (MQ=255)
                                                                                                                               atTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCa                       >  1:2352220/1‑150 (MQ=255)
                                                                                                                                tttGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAg                      <  1:4060726/150‑1 (MQ=255)
                                                                                                                                 ttGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAgg                     <  1:5171818/150‑1 (MQ=255)
                                                                                                                                  tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGc                    >  2:2183218/1‑150 (MQ=255)
                                                                                                                                  tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGc                    >  2:4164974/1‑150 (MQ=255)
                                                                                                                                  tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGc                    >  2:2287713/1‑150 (MQ=255)
                                                                                                                                     cGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGa                 <  2:3533806/150‑1 (MQ=255)
                                                                                                                                       aaaTCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCACATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATg               >  2:4686871/1‑150 (MQ=255)
                                                                                                                                       aaaTCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATg               >  1:2226693/1‑150 (MQ=255)
                                                                                                                                            tGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCAGGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCt          >  2:4793331/1‑150 (MQ=255)
                                                                                                                                            tGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGACTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCt          >  2:975529/1‑150 (MQ=255)
                                                                                                                                              ccGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCa        <  1:3914710/150‑1 (MQ=255)
                                                                                                                                                gtaaCGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGTCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATg      <  2:3109995/148‑1 (MQ=255)
                                                                                                                                                  aaCGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGt    <  1:1613822/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCAGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTCG  >  CP009273/170720‑171017

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: