Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,498,239 A→C 6.2% intergenic (+97/+55) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,498,2390AC 6.2% 730.4 / 9.5 290intergenic (+97/+55)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (107/156);  new base C (21/6);  total (128/162)
Fisher's exact test for biased strand distribution p-value = 3.52e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 305 total aligned reads displayed.

TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  CP009273/3498090‑3498388
                                                                                                                                                     |                                                                                                                                                     
tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCa                                                                                                                                                       <  2:3079258/150‑1 (MQ=255)
tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACTACCCa                                                                                                                                                       >  2:5064176/1‑150 (MQ=255)
tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAAACCTCTCCCCa                                                                                                                                                       >  1:1762749/1‑150 (MQ=255)
    ctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgg                                                                                                                                                   <  1:1857847/150‑1 (MQ=255)
     tGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                                                  <  2:1451114/150‑1 (MQ=255)
      gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAgggg                                                                                                                                                 <  1:4215484/150‑1 (MQ=255)
          cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                             <  1:3744755/150‑1 (MQ=255)
            ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCTAAACCACTCCCCAAAGGGGCGAggg                                                                                                                                           >  1:5158261/1‑150 (MQ=255)
            ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                           <  1:4491962/150‑1 (MQ=255)
            ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                           <  2:604420/150‑1 (MQ=255)
              aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGa                                                                                                                                         <  1:1242166/150‑1 (MQ=255)
               aTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGAc                                                                                                                                        <  1:4961478/150‑1 (MQ=255)
               aTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCACTCCCCAAAGGGGCGAGGGGAc                                                                                                                                        >  1:3134020/1‑150 (MQ=255)
                tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCAAAGGGGCGAGGGGAcc                                                                                                                                       >  1:2932383/1‑150 (MQ=255)
                 cGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGTATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCTTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCg                                                                                                                                      <  1:1962221/150‑1 (MQ=255)
                     tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                  >  2:4610178/1‑150 (MQ=255)
                     tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCACTCCCCAAACGGGCGAGGGGGCCGGTCa                                                                                                                                  >  1:2946340/1‑150 (MQ=39)
                      cACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCAAAACCCAAAGGGGCGAGGGGACCGTTCAg                                                                                                                                 >  1:2872940/1‑150 (MQ=38)
                      cACCTTAACTATTGATAATAAAAGACATTAACCGTAAGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACACCCCAAAGGGGCCAGGGGACCGGTCAg                                                                                                                                 >  2:205930/1‑150 (MQ=39)
                        ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAAACCACACCCCAAATGGGCGCGGGGACCGTTGAGTa                                                                                                                               >  1:263312/1‑150 (MQ=38)
                               tATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAtt                                                                                                                        <  1:5038665/150‑1 (MQ=255)
                                aTTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAACCCACACCCCCAAGGGGGGCGGGGACCCTTCAGTAAgggaatt                                                                                                                       >  1:1033640/1‑145 (MQ=37)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                                                      >  1:3862396/1‑150 (MQ=255)
                                 ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                                                      <  1:1020305/150‑1 (MQ=255)
                                  tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                                     <  1:238342/150‑1 (MQ=255)
                                    ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCt                                                                                                                   <  1:3724018/150‑1 (MQ=255)
                                    ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTAGAATTGTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCt                                                                                                                   <  2:1723618/150‑1 (MQ=255)
                                     taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATGTAGCtt                                                                                                                  >  2:1973292/1‑150 (MQ=255)
                                      aataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTg                                                                                                                 <  1:2844974/150‑1 (MQ=255)
                                      aataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTACCCAAAGGGGCGAGGGGACCGGGCAGTACGGAATTTAGCTTg                                                                                                                 >  2:2532344/1‑150 (MQ=255)
                                       ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCGAACCCGCACCCCCAAGGGGGGAGGGGGCCCGTCAGTACCGAAATTATCTTgt                                                                                                                >  1:94877/1‑149 (MQ=37)
                                       ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGa                                                                                                                <  2:249967/150‑1 (MQ=255)
                                       ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCACACCCCAAAGGGGCGAGGGGACCGTTCAGTAAGGAATTTAGcatta                                                                                                                >  1:4536292/1‑146 (MQ=38)
                                       ataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAATTAGcattg                                                                                                                >  1:2033093/1‑146 (MQ=39)
                                        taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                               >  2:4671244/1‑150 (MQ=255)
                                         aaaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCTCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACg                                                                                                              <  2:869447/150‑1 (MQ=255)
                                          aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACTCCCCAAAGGGGCGAGGGGCCCGTTCAGTACGGAATTTAGCTTGACgg                                                                                                             >  2:1625451/1‑150 (MQ=255)
                                           aaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGt                                                                                                            <  2:2201601/150‑1 (MQ=255)
                                            aGACATTCACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           <  1:1803817/150‑1 (MQ=255)
                                            aGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGAGTCTGGCGCGTTGAATTTGACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTg                                                                                                           <  2:1026978/150‑1 (MQ=255)
                                              aCATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCt                                                                                                         <  1:1597795/150‑1 (MQ=255)
                                              aCATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCt                                                                                                         >  2:3007950/1‑150 (MQ=255)
                                              aCATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACAACCAAACCCTCTACCCAAAGGGGCGAGGGGACCGTTAAGTACGGAATTTAGCTTTACGTTGCt                                                                                                         >  2:1791797/1‑150 (MQ=38)
                                                aTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAt                                                                                                       <  2:3315638/150‑1 (MQ=255)
                                                aTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTAt                                                                                                       >  2:4492115/1‑150 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAACCCTCTACCCCAAGGGGCGCGGGGACCGTTCAGTACGGAAATTAGCTTGACGGTGCTATc                                                                                                      >  2:5011263/1‑150 (MQ=39)
                                                   aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCACTACCCAAAGGGGCGAGGGGACCGTTCAGTTCGGAATTTAGATTTACGGTGCAATACa                                                                                                    >  2:317121/1‑150 (MQ=38)
                                                   aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  1:2285652/150‑1 (MQ=255)
                                                     ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCACCCCAAATGGGCGATGGGACCGTTCAGTACGGAATTTAGCTTGAAGGTGCTATCCAGc                                                                                                  >  1:856918/1‑150 (MQ=255)
                                                      cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                 <  2:153270/150‑1 (MQ=255)
                                                       gTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTAAACCAAACCCTCTCCCCAAACGGGCGAGGGGACCGTTCCGTACGGAAATTAGCTTGACGGTGCTATCCAGCAt                                                                                                >  2:2376039/1‑150 (MQ=39)
                                                        tAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAtt                                                                                               <  2:2296526/150‑1 (MQ=255)
                                                        tAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAtt                                                                                               <  1:4708620/150‑1 (MQ=255)
                                                        tAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAtt                                                                                               <  1:3660817/150‑1 (MQ=255)
                                                             cGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCCGCATTAAATa                                                                                          >  2:3486419/1‑150 (MQ=255)
                                                              ggATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATACCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGCTTGAATGTGCGATCCAGCATTAAATAc                                                                                         >  2:4333021/1‑150 (MQ=39)
                                                               gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                        <  2:4208183/150‑1 (MQ=255)
                                                                     aTGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAAt                                                                                  <  1:4791611/150‑1 (MQ=255)
                                                                     aTGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAAACCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTTGCTTGACGGTGCTATCCAGCCTTAAATACTGGCAAt                                                                                  >  2:2861192/1‑150 (MQ=255)
                                                                     aTGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACCCCCGAACCCACTCCCCCAAGGGGCGCGGGGACCGGTCCGTACCGAAGTTAGCGTGGCGGTGCGATCCCGCCATAAATACGGGCAAt                                                                                  >  1:4637959/1‑150 (MQ=255)
                                                                          ccAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActct                                                                             <  1:569696/150‑1 (MQ=255)
                                                                             gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  2:3634178/150‑1 (MQ=255)
                                                                                tcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGTGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCa                                                                       <  1:781662/150‑1 (MQ=255)
                                                                                   catcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACCCCCTAACCCACAAACCAAAGGGGCGACGGGAACCATAACTACCGAAAATAGCTTTAACGTTCAAAACACCATAAAAAAAAAGCAATAATCACACGCaaaa                                                                    >  1:3269045/1‑150 (MQ=255)
                                                                                      ccGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTACCCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATcc                                                                 <  2:757103/150‑1 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAAAACTGGCCATACTCTCACGCAAAATCCgc                                                               >  2:4866262/1‑150 (MQ=255)
                                                                                          cGATGCTGGCGCGTTGAATTTTACATCCCGTACTTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcag                                                             <  2:4993089/150‑1 (MQ=255)
                                                                                          cGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCACTCCCCCAAGGGGCGAGGGGACCGTTCAGTAAGGAATTTAGCTTGACGGTGCTATACAGCATTAAATACAGGCAATACTCTCACGCAAAATCCgcag                                                             >  2:1510945/1‑150 (MQ=39)
                                                                                             tGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACCGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAc                                                          <  2:3405425/150‑1 (MQ=255)
                                                                                             tGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAc                                                          <  2:1308042/150‑1 (MQ=255)
                                                                                             tGCTGGCGCGATGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAc                                                          <  2:4804535/150‑1 (MQ=255)
                                                                                              gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAcc                                                         <  1:1400060/150‑1 (MQ=255)
                                                                                                tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAAACCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGTCTTTAGCTTTAAGGTGCTATCCAGCATTAAATACTGGGAAAACTCTCACGCAAAATCCGCAGCCCCGt                                                       >  2:5135299/1‑150 (MQ=38)
                                                                                                 ggCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAACCCTATCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATAATGGCAATACTCTCACGCAAAATCCGCAGCACCGtc                                                      >  2:2075728/1‑150 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  2:1355879/150‑1 (MQ=255)
                                                                                                          aaTTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCg                                             <  2:3815608/150‑1 (MQ=255)
                                                                                                             tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                          <  2:1162634/150‑1 (MQ=255)
                                                                                                                aCATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCATTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGa                                       <  2:4176857/150‑1 (MQ=255)
                                                                                                                 cATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGAt                                      <  2:4928571/150‑1 (MQ=255)
                                                                                                                 cATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGAt                                      <  2:2235324/150‑1 (MQ=255)
                                                                                                                      cGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCggg                                 <  2:1094072/150‑1 (MQ=255)
                                                                                                                       gTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                                >  1:2100544/1‑150 (MQ=255)
                                                                                                                        tACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGtt                               >  1:4941654/1‑150 (MQ=255)
                                                                                                                          cGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                             <  1:291426/150‑1 (MQ=255)
                                                                                                                          cGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                             >  1:2613100/1‑150 (MQ=255)
                                                                                                                           gTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAt                            >  1:670385/1‑150 (MQ=255)
                                                                                                                            ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                           <  1:2510964/150‑1 (MQ=255)
                                                                                                                              ccccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCg                         <  1:3377277/150‑1 (MQ=255)
                                                                                                                               cccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGa                        >  1:1374689/1‑150 (MQ=255)
                                                                                                                                  tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACTGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc                     <  2:191615/150‑1 (MQ=255)
                                                                                                                                  tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc                     <  2:4430164/150‑1 (MQ=255)
                                                                                                                                  tCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc                     >  1:2755440/1‑150 (MQ=255)
                                                                                                                                        taccccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACaa               <  2:554930/147‑1 (MQ=255)
                                                                                                                                           cccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  1:3534783/1‑150 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:1489671/1‑150 (MQ=255)
                                                                                                                                              tctcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTGAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCCCGCAAACTCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         <  1:4863440/150‑1 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCGCAATGTTCTGAtt     >  1:4039782/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:521924/1‑150 (MQ=255)
                                                                                                                                                     aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAAGGTTCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg  <  2:4575758/150‑1 (MQ=255)
                                                                                                                                                     aaaGGGGCGAGGGGACCGATCAGTACGCAATTTAGCTTGACGGAGCTAACCAGCATTAAAAACTGGCAATACTCACACGCAAAATCCGCAGCACCCTCTCTTCACCGCTGCGAAACGCGTAAAAGCGAATCTCACACTGTTCAGATtgtc  >  2:2734681/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  CP009273/3498090‑3498388

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: