Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 170,868 C→A 12.2% 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,8680CA12.2% 1075.3 / 93.1 319V168F (GTT→TTT) hemLglutamate‑1‑semialdehyde aminotransferase (aminomutase)
Reads supporting (aligned to +/- strand):  ref base C (144/136);  new base A (21/18);  total (165/154)
Fisher's exact test for biased strand distribution p-value = 8.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 323 total aligned reads displayed.

TGGCAGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTCG  >  CP009273/170719‑171017
                                                                                                                                                     |                                                                                                                                                     
   cAGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGcc                                                                                                                                                    <  1:5378128/150‑1 (MQ=255)
    aGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGccc                                                                                                                                                   >  1:2938795/1‑150 (MQ=255)
    aGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGccc                                                                                                                                                   >  2:272449/1‑150 (MQ=255)
        gcggAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCg                                                                                                                                               <  2:871901/150‑1 (MQ=255)
         cggAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGa                                                                                                                                              <  2:3545615/150‑1 (MQ=255)
               acaGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTtgg                                                                                                                                        >  1:227003/1‑150 (MQ=255)
                 aGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTtggct                                                                                                                                      >  1:5610712/1‑150 (MQ=255)
                  gTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTtggctg                                                                                                                                     <  2:4680496/150‑1 (MQ=255)
                     cATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGcc                                                                                                                                  >  2:4925003/1‑150 (MQ=255)
                        atTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTaa                                                                                                                               >  1:1908162/1‑150 (MQ=255)
                         ttGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAAc                                                                                                                              <  2:24696/150‑1 (MQ=255)
                           gcctgccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGt                                                                                                                            <  2:3291181/150‑1 (MQ=255)
                            cctgccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTg                                                                                                                           <  1:1204653/150‑1 (MQ=255)
                             ctgccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAACTCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACgaga                                                                                                                          >  2:3402696/1‑150 (MQ=255)
                             ctgccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGa                                                                                                                          <  1:4571141/150‑1 (MQ=255)
                              tgccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCTGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAg                                                                                                                         <  2:738296/150‑1 (MQ=255)
                                ccACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTg                                                                                                                       >  1:3485852/1‑150 (MQ=255)
                                  aCCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTgcg                                                                                                                     <  1:2617654/150‑1 (MQ=255)
                                  aCCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTgcg                                                                                                                     >  1:5227324/1‑150 (MQ=255)
                                  aCCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTgcg                                                                                                                     >  1:2438000/1‑150 (MQ=255)
                                     ggCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCa                                                                                                                  <  2:4740883/150‑1 (MQ=255)
                                       cTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGa                                                                                                                >  1:3999895/1‑150 (MQ=255)
                                         cgacgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAAc                                                                                                              >  1:1343595/1‑150 (MQ=255)
                                          gacgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAAcc                                                                                                             >  2:4096285/1‑150 (MQ=255)
                                          gacgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAAcc                                                                                                             >  2:5200277/1‑150 (MQ=255)
                                           acgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCg                                                                                                            >  2:132964/1‑150 (MQ=255)
                                           acgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCg                                                                                                            <  1:3620474/150‑1 (MQ=255)
                                           acgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCg                                                                                                            >  1:3139001/1‑150 (MQ=255)
                                            cgaTAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCgg                                                                                                           >  2:2811352/1‑150 (MQ=255)
                                               taataCAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCtt                                                                                                        <  1:3814342/150‑1 (MQ=255)
                                                          aTCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAgcagg                                                                                             <  1:3345161/150‑1 (MQ=255)
                                                           tctcTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAgcaggc                                                                                            <  2:2686387/150‑1 (MQ=255)
                                                                tGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCa                                                                                       <  2:1711587/150‑1 (MQ=255)
                                                                 gCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGACAGAACCGGCTTTCACCAGCAGGCAGTCAg                                                                                      >  2:2142168/1‑150 (MQ=255)
                                                                    ggTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGt                                                                                   <  1:5130348/150‑1 (MQ=255)
                                                                     gTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTg                                                                                  <  2:974150/150‑1 (MQ=255)
                                                                      tATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGa                                                                                 <  2:4036275/150‑1 (MQ=255)
                                                                       aTTGATCAAATGCGGCGCGTACAGAAGACAGATCATTCTAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGCGTTTGGCTGGCATAACGTGAGTGCGCCAGACCCGGCTTTCACCAGCAGGCAGTCAGCGTGAc                                                                                >  2:3462008/1‑150 (MQ=255)
                                                                          gCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCAt                                                                             <  1:2101340/150‑1 (MQ=255)
                                                                           cTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATg                                                                            >  2:1123884/1‑150 (MQ=255)
                                                                            tCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATgg                                                                           <  1:4994384/150‑1 (MQ=255)
                                                                            tCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATgg                                                                           <  1:1690845/150‑1 (MQ=255)
                                                                              aaaTGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTa                                                                         >  2:4386250/1‑150 (MQ=255)
                                                                              aaaTGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTa                                                                         >  1:1823468/1‑150 (MQ=255)
                                                                               aaTGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTaa                                                                        >  2:5050554/1‑150 (MQ=255)
                                                                                  gcggcgCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAacac                                                                     <  1:4291622/150‑1 (MQ=255)
                                                                                  gcggcgCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAacac                                                                     <  2:2290695/150‑1 (MQ=255)
                                                                                   cggcgCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACAcc                                                                    <  2:3507868/150‑1 (MQ=255)
                                                                                    ggcgCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACAccc                                                                   <  2:4181066/150‑1 (MQ=255)
                                                                                     gcgcGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCt                                                                  <  2:6082/150‑1 (MQ=255)
                                                                                      cgcgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCtt                                                                 >  1:4397159/1‑150 (MQ=255)
                                                                                      cgcgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCtt                                                                 <  2:3176877/150‑1 (MQ=255)
                                                                                      cgcgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCtt                                                                 <  2:5244839/150‑1 (MQ=255)
                                                                                       gcgTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTc                                                                >  2:4817540/1‑150 (MQ=255)
                                                                                          tACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCaaa                                                             <  1:4416562/150‑1 (MQ=255)
                                                                                           aCAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAt                                                            >  2:128986/1‑150 (MQ=255)
                                                                                            cAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                           <  2:3577898/150‑1 (MQ=255)
                                                                                            cAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                           >  2:630617/1‑150 (MQ=255)
                                                                                            cAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                           <  1:1122784/150‑1 (MQ=255)
                                                                                            cAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAAtt                                                           <  1:4245943/150‑1 (MQ=255)
                                                                                                aGCCAGATCATTATACGTACAGGTTAAGGTATATTTGGCGAAATAAGCCGGAACGCCCGGCGAGTTTGGCTGGCCAAACGTGAGTGCTCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATCGTAAAAACCTTCAAAATtaat                                                       >  1:2158811/1‑150 (MQ=255)
                                                                                                aGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaat                                                       <  2:2216412/150‑1 (MQ=255)
                                                                                                aGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaat                                                       <  2:3072422/150‑1 (MQ=255)
                                                                                                 gCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTtaata                                                      <  1:4816/150‑1 (MQ=255)
                                                                                                   cAGATCATTATAAGTACAGGTTAAGTTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTaataat                                                    <  2:3029701/150‑1 (MQ=255)
                                                                                                      aTCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAAtttt                                                 <  2:2148430/150‑1 (MQ=255)
                                                                                                      aTCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAAtttt                                                 <  2:368266/150‑1 (MQ=255)
                                                                                                            ataAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTcgcg                                           >  1:4726431/1‑150 (MQ=255)
                                                                                                             taAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGa                                          <  2:1602382/150‑1 (MQ=255)
                                                                                                              aaGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGAc                                         >  1:1170668/1‑150 (MQ=255)
                                                                                                              aaGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGAc                                         <  2:1142690/150‑1 (MQ=255)
                                                                                                                 tACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCgg                                      >  2:5434175/1‑150 (MQ=255)
                                                                                                                  aCAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGt                                     <  2:7068/150‑1 (MQ=255)
                                                                                                                    aGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaa                                   <  1:4049688/150‑1 (MQ=255)
                                                                                                                     ggTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaaa                                  >  2:1161173/1‑150 (MQ=255)
                                                                                                                      gTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaaaa                                 >  2:2942795/1‑150 (MQ=255)
                                                                                                                      gTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTaaaa                                 <  1:976662/150‑1 (MQ=255)
                                                                                                                       ttAAGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAAc                                >  1:4332417/1‑150 (MQ=255)
                                                                                                                        tAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAAcc                               <  2:3168928/150‑1 (MQ=255)
                                                                                                                          aGGTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCAc                             >  1:638385/1‑150 (MQ=255)
                                                                                                                           ggTATATTTGGCGAAATCTGCCGGAAAGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACg                            >  1:1463423/1‑150 (MQ=255)
                                                                                                                            gTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACgg                           >  2:3903258/1‑150 (MQ=255)
                                                                                                                             tataTTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACggg                          <  1:5029892/150‑1 (MQ=255)
                                                                                                                              atatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGc                         >  2:4446358/1‑150 (MQ=255)
                                                                                                                               tatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTCAAACCACGGGcc                        >  2:4191983/1‑150 (MQ=255)
                                                                                                                               tatTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGcc                        <  2:2869073/150‑1 (MQ=255)
                                                                                                                                atTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCa                       <  1:375089/150‑1 (MQ=255)
                                                                                                                                   tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGc                    >  1:32031/1‑150 (MQ=255)
                                                                                                                                   tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGc                    >  1:31738/1‑150 (MQ=255)
                                                                                                                                   tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGc                    >  1:5269548/1‑150 (MQ=255)
                                                                                                                                   tGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGc                    <  2:3299888/150‑1 (MQ=255)
                                                                                                                                    ggCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCg                   <  2:1779471/150‑1 (MQ=255)
                                                                                                                                        aaaTCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATg               <  1:5526370/150‑1 (MQ=255)
                                                                                                                                          aTCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGgc             <  1:4101812/150‑1 (MQ=255)
                                                                                                                                          aTCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGgc             >  1:832352/1‑150 (MQ=255)
                                                                                                                                             tGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCt          >  1:3310253/1‑150 (MQ=255)
                                                                                                                                                  gAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATgg     >  1:3867119/1‑150 (MQ=255)
                                                                                                                                                  gAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATgg     >  1:783251/1‑150 (MQ=255)
                                                                                                                                                    aCGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTc   <  2:576568/150‑1 (MQ=255)
                                                                                                                                                    aCGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTc   <  1:200295/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCAGCGGCGGAACACAGTTCATATTGCCTGCCACCGGCTCGACGATAATACAGGCAATCTCTTGCGGGTATTGCTCAAATGCGGCGCGTACAGAAGCCAGATCATTATAAGTACAGGTTAAGGTATATTTGGCGAAATCTGCCGGAACGCCCGGCGAGTTTGGCTGGCCTAACGTGAGTGCGCCAGAACCGGCTTTCACCAGCAGGCAGTCAGCGTGACCATGGTAACACCCTTCAAATTTAATAATTTTGTCGCGACCGGTAAAACCACGGGCCAGGCGGATGGCGCTCATGGTCG  >  CP009273/170719‑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: