Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,422,841 C→A 87.4% D293Y (GAT→TAT)  purF ← amidophosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,422,8410CA87.4% 716.2 / 101.9 254D293Y (GAT→TAT) purFamidophosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (16/16);  new base A (109/113);  total (125/129)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

Only 100 of 254 total aligned reads displayed.

CTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCAC  >  CP009273/2422695‑2422987
                                                                                                                                                  |                                                                                                                                                  
ctgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATAtt                                                                                                                                                  <  1:1507754/149‑1 (MQ=255)
ctgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATAtt                                                                                                                                                  <  2:4265630/149‑1 (MQ=255)
 tgctgGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTc                                                                                                                                                 <  2:3403774/149‑1 (MQ=255)
     ggCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCAt                                                                                                                                             <  2:1042079/149‑1 (MQ=255)
          ggCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTcgcg                                                                                                                                        <  1:3426273/149‑1 (MQ=255)
             aTGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGc                                                                                                                                     >  2:2118404/1‑149 (MQ=255)
               gATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCaa                                                                                                                                   >  1:2710844/1‑149 (MQ=255)
                    aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtc                                                                                                                              >  1:3245306/1‑149 (MQ=255)
                    aGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTtc                                                                                                                              <  2:2952290/149‑1 (MQ=255)
                     ggTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTtct                                                                                                                             <  1:178980/149‑1 (MQ=255)
                           gCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCa                                                                                                                       <  2:4111624/149‑1 (MQ=255)
                             cAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGt                                                                                                                     <  2:495269/149‑1 (MQ=255)
                             cAACATAGAGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGt                                                                                                                     <  2:4596935/149‑1 (MQ=255)
                              aaCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGtt                                                                                                                    >  2:25692/1‑149 (MQ=255)
                               aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGttt                                                                                                                   >  2:1139395/1‑149 (MQ=255)
                               aCATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGttt                                                                                                                   <  1:1335845/149‑1 (MQ=255)
                                 aTAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCg                                                                                                                 >  2:1869646/1‑149 (MQ=255)
                                  tAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGt                                                                                                                <  1:1553686/149‑1 (MQ=255)
                                  tAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGt                                                                                                                >  1:4529333/1‑149 (MQ=255)
                                  tAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGt                                                                                                                <  1:3889357/149‑1 (MQ=255)
                                    gCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGc                                                                                                              >  2:2501796/1‑149 (MQ=255)
                                         ttttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCata                                                                                                         <  1:1682998/149‑1 (MQ=255)
                                           ttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATAtt                                                                                                       <  2:1162233/149‑1 (MQ=255)
                                           ttAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATAtt                                                                                                       <  2:1296684/149‑1 (MQ=255)
                                            tAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTc                                                                                                      >  1:80152/1‑149 (MQ=255)
                                             aaCGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTca                                                                                                     <  1:1405107/149‑1 (MQ=255)
                                              aCGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTcac                                                                                                    <  1:4010355/149‑1 (MQ=255)
                                              aCGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTcac                                                                                                    <  2:4462065/149‑1 (MQ=255)
                                               cGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTcaca                                                                                                   >  2:839750/1‑149 (MQ=255)
                                                gAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCacac                                                                                                  >  1:4112636/1‑149 (MQ=255)
                                                 aaGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcg                                                                                                 >  1:183142/1‑149 (MQ=255)
                                                   gCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACAcgcg                                                                                               <  2:1001409/149‑1 (MQ=255)
                                                      cTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACgcgcgc                                                                                            >  2:1379700/1‑149 (MQ=255)
                                                      cTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACgcgcgc                                                                                            <  1:3370741/149‑1 (MQ=255)
                                                           ggTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTaa                                                                                       >  1:2896367/1‑149 (MQ=255)
                                                            gTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTaaa                                                                                      <  1:1597428/149‑1 (MQ=255)
                                                            gTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTaaa                                                                                      <  1:1597472/149‑1 (MQ=255)
                                                             tACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAAc                                                                                     >  1:420774/1‑149 (MQ=255)
                                                                ggTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGa                                                                                  >  1:3971932/1‑149 (MQ=255)
                                                                 gTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaa                                                                                 >  2:650876/1‑149 (MQ=255)
                                                                 gTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaa                                                                                 >  2:3041502/1‑149 (MQ=255)
                                                                  tttGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGaaa                                                                                >  2:1510975/1‑149 (MQ=255)
                                                                   ttGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAAt                                                                               <  1:3569693/149‑1 (MQ=255)
                                                                     gCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAAttt                                                                             <  1:3453855/149‑1 (MQ=255)
                                                                      cccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAAtttt                                                                            <  1:2117361/149‑1 (MQ=255)
                                                                       ccAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTg                                                                           >  1:3729393/1‑149 (MQ=255)
                                                                        cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                          <  2:2438735/149‑1 (MQ=255)
                                                                        cAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGt                                                                          <  1:1888523/149‑1 (MQ=255)
                                                                          gAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCg                                                                        >  1:3655094/1‑149 (MQ=255)
                                                                           aaTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGa                                                                       >  1:3791428/1‑149 (MQ=255)
                                                                            aTACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGAt                                                                      >  2:1775587/1‑149 (MQ=255)
                                                                             tACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATa                                                                     <  1:970126/149‑1 (MQ=255)
                                                                              aCGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATaa                                                                    <  2:3654338/149‑1 (MQ=255)
                                                                                gAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAAc                                                                  >  1:2579759/1‑149 (MQ=255)
                                                                                   cAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAg                                                               <  2:2262341/149‑1 (MQ=255)
                                                                                   cAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAg                                                               <  2:3054506/149‑1 (MQ=255)
                                                                                    aaTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGt                                                              >  2:4407286/1‑149 (MQ=255)
                                                                                     aTTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTc                                                             <  1:605838/149‑1 (MQ=255)
                                                                                        tCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgg                                                          <  2:4175493/149‑1 (MQ=255)
                                                                                         ccAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcggg                                                         <  2:2143648/149‑1 (MQ=255)
                                                                                          cAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggc                                                        >  1:1956834/1‑149 (MQ=255)
                                                                                            gcgcGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcgg                                                      >  1:635007/1‑149 (MQ=255)
                                                                                            gcgcGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCcgggcgg                                                      >  2:2197465/1‑149 (MQ=255)
                                                                                              gcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCgggcgggc                                                    <  2:3952587/149‑1 (MQ=255)
                                                                                              gcgATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCgggcgggc                                                    <  1:3062579/149‑1 (MQ=255)
                                                                                                 atatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCaaa                                                 >  1:4621665/1‑149 (MQ=255)
                                                                                                  tatCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAg                                                <  2:275149/149‑1 (MQ=255)
                                                                                                   atCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGt                                               <  2:3659367/149‑1 (MQ=255)
                                                                                                    tCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGta                                              <  1:706644/149‑1 (MQ=255)
                                                                                                    tCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGta                                              <  1:1120515/149‑1 (MQ=255)
                                                                                                    tCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGta                                              >  2:3932659/1‑149 (MQ=255)
                                                                                                      acacGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtata                                            >  1:1383757/1‑149 (MQ=255)
                                                                                                      acacGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGtata                                            >  1:2104815/1‑149 (MQ=255)
                                                                                                          gAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacata                                        <  1:4412564/149‑1 (MQ=255)
                                                                                                           aGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacatac                                       >  2:1533592/1‑149 (MQ=255)
                                                                                                           aGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacatac                                       >  1:2405317/1‑149 (MQ=255)
                                                                                                           aGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacatac                                       <  1:2672659/149‑1 (MQ=255)
                                                                                                           aGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTatacatac                                       <  2:3511692/149‑1 (MQ=255)
                                                                                                              tttCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCa                                    >  1:4039881/1‑149 (MQ=255)
                                                                                                               ttCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCaa                                   <  2:3730425/149‑1 (MQ=255)
                                                                                                                tCTGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCaaa                                  <  1:3507331/149‑1 (MQ=255)
                                                                                                                  tGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACa                                >  1:3391321/1‑149 (MQ=255)
                                                                                                                  tGGGATCGGGATCACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACa                                <  1:4388023/149‑1 (MQ=255)
                                                                                                                    ggatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAgg                              <  1:1429450/149‑1 (MQ=255)
                                                                                                                    ggatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAgg                              <  1:869821/149‑1 (MQ=255)
                                                                                                                     gatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGc                             >  2:1197631/1‑149 (MQ=255)
                                                                                                                     gatcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGc                             >  2:1784197/1‑149 (MQ=255)
                                                                                                                      atcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTCAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCa                            >  2:632578/1‑149 (MQ=255)
                                                                                                                      atcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCACACAGGCa                            >  2:3547577/1‑149 (MQ=255)
                                                                                                                      atcgggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCa                            <  2:2087076/149‑1 (MQ=255)
                                                                                                                       tcgggatcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCAc                           >  2:323928/1‑149 (MQ=255)
                                                                                                                          ggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGa                        >  1:448878/1‑149 (MQ=255)
                                                                                                                          ggatcACCACGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATAAATACTCAAAAAGGCACGGa                        >  1:2192228/1‑149 (MQ=255)
                                                                                                                             tcACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTg                     >  2:3218140/1‑149 (MQ=255)
                                                                                                                               accacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCt                   >  1:1007936/1‑149 (MQ=255)
                                                                                                                                 cacGTCGATATCCAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGa                 <  1:3200757/149‑1 (MQ=255)
                                                                                                                                            ccAGATATTCCCATTCGCGGGCAATTTTCTCGCCCCGTTTCGTGCCCATATTCACACGCGCGCTGGAAACGGAATTTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCa      <  1:2847782/149‑1 (MQ=255)
                                                                                                                                            ccAGATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCa      <  1:505292/149‑1 (MQ=255)
                                                                                                                                               gATATTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGca   >  2:1108195/1‑149 (MQ=255)
                                                                                                                                                aTCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGcac  >  1:1918507/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTGCTGGCCCGGCATGATAAAGGTGCGGCCAACATAGCGGTTTTTAACGAAGCCCTGGCGGTACGGTTTGCCCAGAATACGAGCAATTTCCAGCGCGATATCACACGAGGTTTCTGGGATCGGGATCACCACGTCGATATCCAGATCTTCCCATTCGCGGGCAATTTTCTCGCCCAGTTTCGTGCCCATATTCACACGCGCGCTGTAAACGGAAATTTTGTCGATAAACGAGTCCGGGCGGGCAAAGTATACATACTCAAACAGGCACGGATTGCTGACCGGATTGTCAGCAC  >  CP009273/2422695‑2422987

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: ‑