Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,902,059 C→T 12.6% V323I (GTC→ATC)  pyrG ← CTP synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,902,0590CT12.6% 825.9 / 82.3 270V323I (GTC→ATC) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (105/131);  new base T (9/25);  total (114/156)
Fisher's exact test for biased strand distribution p-value = 6.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 274 total aligned reads displayed.

ATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACA  >  CP009273/2901910‑2902208
                                                                                                                                                     |                                                                                                                                                     
atTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGAc                                                                                                                                                       >  1:1992473/1‑150 (MQ=255)
   gTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCt                                                                                                                                                    <  1:441159/150‑1 (MQ=255)
      ctcACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGac                                                                                                                                                 <  2:1538212/150‑1 (MQ=255)
       tcACGCGCAAAACTCGCGGTCGTAATCAAGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGACCAGAACTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGaca                                                                                                                                                >  1:3808174/1‑150 (MQ=255)
        cACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGacac                                                                                                                                               <  2:1965355/150‑1 (MQ=255)
         aCGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACg                                                                                                                                              >  1:635508/1‑150 (MQ=255)
         aCGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACg                                                                                                                                              <  1:265961/150‑1 (MQ=255)
            cgcAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAtt                                                                                                                                           >  1:671346/1‑150 (MQ=255)
                 aaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCa                                                                                                                                      <  2:657954/150‑1 (MQ=255)
                  aCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAg                                                                                                                                     <  1:3271323/150‑1 (MQ=255)
                    gcgcgGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGcc                                                                                                                                   <  2:75008/150‑1 (MQ=255)
                           cGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCgtg                                                                                                                            >  1:1301404/1‑150 (MQ=255)
                            gTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCgtgt                                                                                                                           <  2:6077/150‑1 (MQ=255)
                             tAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGtt                                                                                                                          <  1:3135576/150‑1 (MQ=255)
                              aaTCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAGTCTTCAGCCCACCGTGttt                                                                                                                         >  2:3347680/1‑150 (MQ=255)
                                  aTGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGt                                                                                                                     >  1:3175715/1‑150 (MQ=255)
                                      cTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCtt                                                                                                                 <  1:4563346/150‑1 (MQ=255)
                                       ttCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTc                                                                                                                <  2:1100172/150‑1 (MQ=255)
                                       ttCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTc                                                                                                                <  1:1818271/150‑1 (MQ=255)
                                        tCTACGCCCCGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTTAGGATTTCAACGCCGCGCGTTTCACCATCTTGTGAATCGATCAGTTTGATGTTGTTGCTGACACGATTGTTCAGCCCACCGTGTTTCAGTGCTTCg                                                                                                               <  1:871838/150‑1 (MQ=255)
                                        tCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCg                                                                                                               <  2:4043943/150‑1 (MQ=255)
                                        tCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCATTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCg                                                                                                               <  2:3963749/150‑1 (MQ=255)
                                          tACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGAt                                                                                                             <  1:4129100/150‑1 (MQ=255)
                                           aCGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATc                                                                                                            <  2:3840810/150‑1 (MQ=255)
                                           aCGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATc                                                                                                            >  1:1013950/1‑150 (MQ=255)
                                              ccACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACt                                                                                                         <  2:2297052/150‑1 (MQ=255)
                                                 cGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGAt                                                                                                      <  2:2812041/150‑1 (MQ=255)
                                                 cGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGAt                                                                                                      <  2:574041/150‑1 (MQ=255)
                                                  gATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGAtt                                                                                                     >  1:3309792/1‑150 (MQ=255)
                                                  gATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGAtt                                                                                                     >  1:38379/1‑150 (MQ=255)
                                                     aGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTtat                                                                                                  <  2:379214/150‑1 (MQ=255)
                                                         gAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGc                                                                                              >  1:176001/1‑150 (MQ=255)
                                                         gAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGc                                                                                              >  1:1174780/1‑150 (MQ=255)
                                                         gAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGGTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGc                                                                                              >  1:3422406/1‑150 (MQ=255)
                                                            aCCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATc                                                                                           <  2:1067147/150‑1 (MQ=255)
                                                            aCCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATc                                                                                           <  2:210931/150‑1 (MQ=255)
                                                              cGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCg                                                                                         <  1:1374852/150‑1 (MQ=255)
                                                                ccAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGc                                                                                       <  1:3071214/150‑1 (MQ=255)
                                                                 cAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCa                                                                                      <  1:1038322/150‑1 (MQ=255)
                                                                 cAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCa                                                                                      >  2:233567/1‑150 (MQ=255)
                                                                  aGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAg                                                                                     >  2:4347065/1‑150 (MQ=255)
                                                                  aGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAg                                                                                     <  2:3417615/150‑1 (MQ=255)
                                                                   ggTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGt                                                                                    <  1:827194/150‑1 (MQ=255)
                                                                    gTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGtt                                                                                   <  2:1733805/150‑1 (MQ=255)
                                                                     tACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTNGATCACTGATTTATAAGCATCCGGCAGTTc                                                                                  <  1:4277063/150‑1 (MQ=255)
                                                                          ggATTGCGTCCAGACCTTTAAGGATTTCAACGCNGCGCGTTTCNACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGt                                                                             <  1:1421985/150‑1 (MQ=255)
                                                                              tGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACtt                                                                         <  1:94148/150‑1 (MQ=255)
                                                                                 gTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGcc                                                                      >  2:2019307/1‑150 (MQ=255)
                                                                                 gTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGcc                                                                      >  1:335395/1‑150 (MQ=255)
                                                                                  tCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCg                                                                     >  2:3431595/1‑150 (MQ=255)
                                                                                  tCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCg                                                                     <  1:4488842/150‑1 (MQ=255)
                                                                                  tCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCg                                                                     >  1:223795/1‑150 (MQ=255)
                                                                                    cAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGAc                                                                   <  1:1449693/150‑1 (MQ=255)
                                                                                    cAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGAc                                                                   <  1:4635370/150‑1 (MQ=255)
                                                                                     aGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGAcc                                                                  <  1:2437792/150‑1 (MQ=255)
                                                                                      gACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCa                                                                 >  1:1076229/1‑150 (MQ=255)
                                                                                         cTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAc                                                              <  1:2448744/150‑1 (MQ=255)
                                                                                         cTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAc                                                              >  1:1777545/1‑150 (MQ=255)
                                                                                          tttAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAcc                                                             <  1:2326080/150‑1 (MQ=255)
                                                                                                  tttCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGa                                                     >  1:4618593/1‑150 (MQ=255)
                                                                                                  tttCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGa                                                     <  2:755803/150‑1 (MQ=255)
                                                                                                   ttCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGAc                                                    >  1:1647059/1‑150 (MQ=255)
                                                                                                    tCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACt                                                   <  1:1212638/150‑1 (MQ=255)
                                                                                                    tCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACt                                                   <  2:3912256/150‑1 (MQ=255)
                                                                                                     cAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACtt                                                  <  1:2106025/150‑1 (MQ=255)
                                                                                                     cAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACtt                                                  >  1:3697927/1‑150 (MQ=255)
                                                                                                         gccgcGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGTCGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACt                                              <  2:4579878/150‑1 (MQ=255)
                                                                                                           cgcgcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTa                                            <  1:2710547/150‑1 (MQ=255)
                                                                                                           cgcgcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTa                                            <  2:1267737/150‑1 (MQ=255)
                                                                                                           cgcgcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTa                                            <  2:3222671/150‑1 (MQ=255)
                                                                                                            gcgcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTAc                                           <  1:2826463/150‑1 (MQ=255)
                                                                                                             cgcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCTCCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTAcc                                          <  1:2912392/150‑1 (MQ=255)
                                                                                                             cgcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTAcc                                          <  1:598331/150‑1 (MQ=255)
                                                                                                               cgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCgg                                        >  1:2581785/1‑150 (MQ=255)
                                                                                                                gTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCggg                                       <  2:2090253/150‑1 (MQ=255)
                                                                                                                gTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCggg                                       >  2:1467558/1‑150 (MQ=255)
                                                                                                                   tCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTc                                    <  2:2401432/150‑1 (MQ=255)
                                                                                                                   tCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTc                                    <  2:4570604/150‑1 (MQ=255)
                                                                                                                   tCAACATCGTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTc                                    <  2:2639115/150‑1 (MQ=255)
                                                                                                                   acAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTc                                    <  2:3591397/149‑1 (MQ=255)
                                                                                                                    cAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCg                                   >  2:3036830/1‑150 (MQ=255)
                                                                                                                       cATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGctt                                <  1:4201106/150‑1 (MQ=255)
                                                                                                                        aTCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttc                               <  2:2628558/150‑1 (MQ=255)
                                                                                                                          cTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttctt                             <  2:1653302/150‑1 (MQ=255)
                                                                                                                          cTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttctt                             <  1:740257/150‑1 (MQ=255)
                                                                                                                          cTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttctt                             >  1:1537128/1‑150 (MQ=255)
                                                                                                                               gAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAg                        >  1:815407/1‑150 (MQ=255)
                                                                                                                               gAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAg                        >  1:3736136/1‑150 (MQ=255)
                                                                                                                                 aTCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGAt                      >  2:4634320/1‑150 (MQ=255)
                                                                                                                                  tCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATa                     <  2:2121738/150‑1 (MQ=255)
                                                                                                                                    gATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAAc                   >  2:2552713/1‑150 (MQ=255)
                                                                                                                                     aTCAGTTTGATGTTGACGCTGACACGTTTCTTCAGCACACCGTGTTTCAGTGCGTCGATCACTGATTTATAAGAATCCGGCAGTTCAAGGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAAcc                  <  2:1923434/150‑1 (MQ=255)
                                                                                                                                       cAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTg                <  1:4094101/150‑1 (MQ=255)
                                                                                                                                         gTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTGCTTCGAAGATAACCTGtt              >  2:414481/1‑150 (MQ=255)
                                                                                                                                             gATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCa          >  1:556087/1‑150 (MQ=255)
                                                                                                                                              aTGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCAt         <  1:4387439/150‑1 (MQ=255)
                                                                                                                                              aTGTTGACGCTGACAAGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGTCTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCAt         <  1:2047511/150‑1 (MQ=255)
                                                                                                                                                  tGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCgg     >  2:1365728/1‑150 (MQ=255)
                                                                                                                                                     cGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACa  <  1:1945432/150‑1 (MQ=255)
                                                                                                                                                     cGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACGTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACa  <  2:1779708/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACA  >  CP009273/2901910‑2902208

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