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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,902,0590CT14.7% 908.9 / 116.2 312V323I (GTC→ATC) pyrGCTP synthetase
Reads supporting (aligned to +/- strand):  ref base C (130/136);  new base T (15/31);  total (145/167)
Fisher's exact test for biased strand distribution p-value = 5.40e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 319 total aligned reads displayed.

ATTGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACA  >  CP009273/2901910‑2902208
                                                                                                                                                     |                                                                                                                                                     
 ttGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACg                                                                                                                                                      >  2:1870294/1‑150 (MQ=255)
  tGTTCTCACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGc                                                                                                                                                     >  1:2387829/1‑150 (MQ=255)
      ctcACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGac                                                                                                                                                 <  2:3130998/150‑1 (MQ=255)
       tcACGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGaca                                                                                                                                                <  1:2842090/150‑1 (MQ=255)
         aCGCGCAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACg                                                                                                                                              >  1:1011623/1‑150 (MQ=255)
           gcgcAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGag                                                                                                                                            >  2:2849386/1‑149 (MQ=255)
            cgcAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAtt                                                                                                                                           <  1:3494436/150‑1 (MQ=255)
             gcAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttc                                                                                                                                          <  2:2734415/150‑1 (MQ=255)
              cAAAACGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttct                                                                                                                                         <  1:3302245/150‑1 (MQ=255)
                aaaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttcttc                                                                                                                                       >  1:1364505/1‑150 (MQ=255)
                aaaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGAttcttc                                                                                                                                       >  1:27067/1‑150 (MQ=255)
                 aaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCa                                                                                                                                      <  2:1603695/150‑1 (MQ=255)
                 aaCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCa                                                                                                                                      <  1:1798204/150‑1 (MQ=255)
                  aCGCGCGGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAg                                                                                                                                     <  1:3819971/150‑1 (MQ=255)
                   cgcgcgGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGc                                                                                                                                    <  1:2407738/150‑1 (MQ=255)
                     cgcgGTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGccc                                                                                                                                  <  2:1005209/150‑1 (MQ=255)
                         gTCGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCg                                                                                                                              >  1:4334505/1‑150 (MQ=255)
                           cGTAATCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCgtg                                                                                                                            >  2:1414685/1‑150 (MQ=255)
                           cGTAATCATGCCTTCTACGCCACGATAGCCGAAAACGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCTATCAGTTTGATTTTGACGCTGACACGATTCTTCAGCCCACCgtg                                                                                                                            >  2:1414476/1‑150 (MQ=255)
                              aaTCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGttt                                                                                                                         >  1:4211963/1‑150 (MQ=255)
                              aaTCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGttt                                                                                                                         <  1:4492623/150‑1 (MQ=255)
                               aTCATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTc                                                                                                                        >  1:2102385/1‑150 (MQ=255)
                                 cATGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAg                                                                                                                      <  1:3943086/150‑1 (MQ=255)
                                   tGCCTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTg                                                                                                                    <  1:123735/150‑1 (MQ=255)
                                     ccTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCt                                                                                                                  <  2:737152/150‑1 (MQ=255)
                                      cTTCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCtt                                                                                                                 <  1:746682/150‑1 (MQ=255)
                                       ttCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTc                                                                                                                <  1:539034/150‑1 (MQ=255)
                                        tCTACGCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCg                                                                                                               >  2:1758342/1‑150 (MQ=255)
                                             gCCACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCAc                                                                                                          >  1:1933488/1‑150 (MQ=255)
                                              ccACGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACt                                                                                                         <  1:1721124/150‑1 (MQ=255)
                                                 cGATAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGAt                                                                                                      >  1:4140425/1‑150 (MQ=255)
                                                   aTAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGAttt                                                                                                    <  2:1896621/150‑1 (MQ=255)
                                                    tAGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTta                                                                                                   <  1:940444/150‑1 (MQ=255)
                                                     aGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTCAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTtat                                                                                                  <  1:4352169/150‑1 (MQ=255)
                                                     aGCCGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTtat                                                                                                  >  2:4092276/1‑150 (MQ=255)
                                                       ccGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATaa                                                                                                >  2:2289549/1‑150 (MQ=255)
                                                        cGAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAg                                                                                               <  2:868405/150‑1 (MQ=255)
                                                         gAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGc                                                                                              >  1:3703605/1‑150 (MQ=255)
                                                         gAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGc                                                                                              >  1:309285/1‑150 (MQ=255)
                                                         gAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGc                                                                                              <  1:2504304/150‑1 (MQ=255)
                                                         gAAACCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGc                                                                                              >  2:3289609/1‑150 (MQ=255)
                                                          aaaCCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCa                                                                                             <  2:793140/150‑1 (MQ=255)
                                                          aaaCCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCa                                                                                             >  2:181693/1‑150 (MQ=255)
                                                           aaCCGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCAt                                                                                            >  2:694242/1‑150 (MQ=255)
                                                           aaCCGCCAGGTACGAGGATTGCGTCCAGACCATTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCAt                                                                                            <  1:1198256/150‑1 (MQ=255)
                                                              cGCCAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACTCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCg                                                                                         <  2:1734442/150‑1 (MQ=255)
                                                                 cAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCa                                                                                      <  1:169670/150‑1 (MQ=255)
                                                                 cAGGTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCa                                                                                      >  1:3869023/1‑150 (MQ=255)
                                                                    gTACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGtt                                                                                   <  1:2450970/150‑1 (MQ=255)
                                                                     tACGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTc                                                                                  <  1:2170406/150‑1 (MQ=255)
                                                                       cGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCaa                                                                                >  1:2069014/1‑150 (MQ=255)
                                                                       cGAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCaa                                                                                >  1:3410441/1‑150 (MQ=255)
                                                                        gAGGATTGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAAt                                                                               <  1:585151/150‑1 (MQ=255)
                                                                             ttGCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACt                                                                          <  1:3406800/150‑1 (MQ=255)
                                                                               gCGTCCAGACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTg                                                                        <  2:2131149/150‑1 (MQ=255)
                                                                                      gACCTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCa                                                                 >  1:4237681/1‑150 (MQ=255)
                                                                                         cTTTAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAc                                                              <  2:3133472/150‑1 (MQ=255)
                                                                                          tttAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAcc                                                             <  2:1335438/150‑1 (MQ=255)
                                                                                          tttAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATAcc                                                             <  2:196098/150‑1 (MQ=255)
                                                                                            tAAGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGa                                                           >  1:1529569/1‑150 (MQ=255)
                                                                                             aaGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGAt                                                          >  1:762868/1‑150 (MQ=255)
                                                                                             aaGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGAt                                                          <  2:2920616/150‑1 (MQ=255)
                                                                                              aGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATg                                                         <  1:146704/150‑1 (MQ=255)
                                                                                              aGGATTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATg                                                         >  2:4302325/1‑150 (MQ=255)
                                                                                                 aTTTCAACGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTg                                                      >  1:843391/1‑150 (MQ=255)
                                                                                                      aaCGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTc                                                 >  1:3999949/1‑150 (MQ=255)
                                                                                                       aCGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCa                                                >  2:4558287/1‑150 (MQ=255)
                                                                                                       aCGCCGCGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCa                                                >  2:3768071/1‑150 (MQ=255)
                                                                                                        cgccgcGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCAc                                               <  2:2187080/150‑1 (MQ=255)
                                                                                                         gccgcGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACt                                              <  1:639329/150‑1 (MQ=255)
                                                                                                         gccgcGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACt                                              >  1:514611/1‑150 (MQ=255)
                                                                                                         gccgcGCGTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACt                                              <  1:1088445/150‑1 (MQ=255)
                                                                                                          ccgcGCGTTTCAACATCTTGTGACTCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACtt                                             <  2:2376496/150‑1 (MQ=255)
                                                                                                           cgcgcgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTa                                            >  2:118950/1‑150 (MQ=255)
                                                                                                               cgTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCgg                                        <  1:2363912/150‑1 (MQ=255)
                                                                                                                gTTTCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCggg                                       >  2:2475514/1‑150 (MQ=255)
                                                                                                                   tCAACATCTTGTGAATCGATCAGTTTGATGTTGATGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTc                                    <  1:1870440/150‑1 (MQ=255)
                                                                                                                   tCAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTc                                    <  1:2614591/150‑1 (MQ=255)
                                                                                                                    cAACATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCg                                   >  2:601636/1‑150 (MQ=255)
                                                                                                                     aaCATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGc                                  >  2:2755035/1‑150 (MQ=255)
                                                                                                                     aaCATCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGc                                  >  2:3599120/1‑150 (MQ=255)
                                                                                                                        aTCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttc                               <  1:358/150‑1 (MQ=255)
                                                                                                                        aTCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttc                               <  1:2900441/150‑1 (MQ=255)
                                                                                                                         tCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttct                              >  1:4208356/1‑150 (MQ=255)
                                                                                                                         tCTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttct                              <  1:2654727/150‑1 (MQ=255)
                                                                                                                          cTTGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGcttctt                             <  2:1043428/150‑1 (MQ=255)
                                                                                                                           ttGTGAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCttcttc                            >  2:2872217/1‑150 (MQ=255)
                                                                                                                            tgtgAATCGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCg                           <  1:2274438/150‑1 (MQ=255)
                                                                                                                                   cGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATaa                    <  1:249952/150‑1 (MQ=255)
                                                                                                                                   cGATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATaa                    >  2:2593542/1‑150 (MQ=255)
                                                                                                                                    gATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAAc                   >  2:2414392/1‑150 (MQ=255)
                                                                                                                                    gATCAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAAc                   >  2:2662976/1‑150 (MQ=255)
                                                                                                                                       cAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTg                <  1:4471190/150‑1 (MQ=255)
                                                                                                                                       cAGTTTGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTg                >  2:2127472/1‑150 (MQ=255)
                                                                                                                                            tGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTccc           >  2:3948188/1‑150 (MQ=255)
                                                                                                                                            tGATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTccc           >  2:4398966/1‑150 (MQ=255)
                                                                                                                                             gATGTTGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCa          <  2:256254/150‑1 (MQ=255)
                                                                                                                                                 ttGACGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCg      <  2:4502677/150‑1 (MQ=255)
                                                                                                                                                     cGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACa  <  2:4230703/150‑1 (MQ=255)
                                                                                                                                                     cGCTGACACGATTCTTCAGCCCACCGTGTTTCAGTGCTTCGATCACTGATTTATAAGCATCCGGCAGTTCAATGTACTTGCCGACCATACCGATGGTGACTTCACTTACCGGGTTCGCTTCTTCGAAGATAACCTGTTCCCATTCGGACa  >  2:3461592/1‑150 (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: