Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,523,001 T→A 8.6% intergenic (‑1/+55) yjhP ← / ← yjhQ putative methyltransferase/putative acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,523,0010TA 8.6% 778.9 / 51.8 233intergenic (‑1/+55)yjhP/yjhQputative methyltransferase/putative acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (98/115);  new base A (9/11);  total (107/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

Only 100 of 243 total aligned reads displayed.

ATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAG  >  CP009273/4522852‑4523149
                                                                                                                                                     |                                                                                                                                                    
  ctcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCattat                                                                                                                                                    <  2:4266603/150‑1 (MQ=255)
  ctcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCattat                                                                                                                                                    >  1:1863263/1‑150 (MQ=255)
   tcGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCattatt                                                                                                                                                   >  1:1062987/1‑150 (MQ=255)
      cccGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTaa                                                                                                                                                <  2:3253191/150‑1 (MQ=255)
      cccGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGGTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTaa                                                                                                                                                <  2:3253005/150‑1 (MQ=255)
        cGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCAtaatttaatt                                                                                                                                              <  2:1097704/150‑1 (MQ=255)
        cGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAAtt                                                                                                                                              <  2:1451648/150‑1 (MQ=255)
        cGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAAtt                                                                                                                                              <  2:3736504/150‑1 (MQ=255)
         gAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTc                                                                                                                                             <  2:1671312/150‑1 (MQ=255)
         gAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCGTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTc                                                                                                                                             <  2:2869865/150‑1 (MQ=255)
           gCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCt                                                                                                                                           >  1:2092750/1‑150 (MQ=255)
            ccGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTg                                                                                                                                          <  2:2733445/150‑1 (MQ=255)
                tGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTAt                                                                                                                                      >  2:3263856/1‑150 (MQ=255)
                   cGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTgg                                                                                                                                   <  2:216863/150‑1 (MQ=255)
                   cGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTgg                                                                                                                                   >  1:427681/1‑150 (MQ=255)
                    gAGGTCTAGAATGCGGGTGCCCGGATTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTggg                                                                                                                                  <  1:2468863/150‑1 (MQ=255)
                    gAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTggg                                                                                                                                  <  1:2797952/150‑1 (MQ=255)
                     aGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTTAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGt                                                                                                                                 <  1:465780/150‑1 (MQ=255)
                     aGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGACCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGt                                                                                                                                 >  2:1393291/1‑150 (MQ=255)
                         cAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAg                                                                                                                             <  1:2038101/150‑1 (MQ=255)
                         cAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAg                                                                                                                             <  1:1506486/150‑1 (MQ=255)
                           aGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGcc                                                                                                                           >  2:3537368/1‑150 (MQ=255)
                            gAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCt                                                                                                                          <  2:3390450/150‑1 (MQ=255)
                            gAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCt                                                                                                                          <  2:2982429/150‑1 (MQ=255)
                                gCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTaaa                                                                                                                      <  2:3683588/150‑1 (MQ=255)
                                gCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTaaa                                                                                                                      <  1:3972910/150‑1 (MQ=255)
                                    gTGCCCGGCTTCATTCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGTGTTATGGATGCGGTGTTCACTTTCTCTGATGGTAAAAATACGTGGGATATCAATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAgg                                                                                                                  <  1:2566129/150‑1 (MQ=255)
                                    gTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAgg                                                                                                                  <  1:331495/150‑1 (MQ=255)
                                    gTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAgg                                                                                                                  <  1:3505634/150‑1 (MQ=255)
                                     tGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGa                                                                                                                 >  2:2792668/1‑150 (MQ=255)
                                     tGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGa                                                                                                                 <  1:1395860/150‑1 (MQ=255)
                                      gCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGac                                                                                                                <  1:2415765/150‑1 (MQ=255)
                                       cccGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGaca                                                                                                               >  1:1268090/1‑150 (MQ=255)
                                           gCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCGTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaac                                                                                                           <  2:4358881/150‑1 (MQ=255)
                                             ttCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACAcaacaa                                                                                                         >  1:989425/1‑150 (MQ=255)
                                              tCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAg                                                                                                        <  2:3568824/150‑1 (MQ=255)
                                              tCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAg                                                                                                        >  1:4485239/1‑150 (MQ=255)
                                                aTGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGAt                                                                                                      <  2:3134909/150‑1 (MQ=255)
                                                   cgcAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGt                                                                                                   >  2:1484253/1‑150 (MQ=255)
                                                     cagcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTa                                                                                                 <  2:4039422/150‑1 (MQ=255)
                                                       gcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTAtt                                                                                               <  1:2636637/150‑1 (MQ=255)
                                                         aCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCa                                                                                             <  2:3225598/150‑1 (MQ=255)
                                                         aCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCa                                                                                             <  1:698198/150‑1 (MQ=255)
                                                         aCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCa                                                                                             <  1:216619/150‑1 (MQ=255)
                                                            cgACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTACTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGt                                                                                          >  1:1887420/1‑150 (MQ=255)
                                                             gACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                         >  2:165175/1‑150 (MQ=255)
                                                             gACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                         <  1:3135753/150‑1 (MQ=255)
                                                              aCCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAATAAGGACACAACAAGATTGTTATTCACGTCa                                                                                        <  1:617964/150‑1 (MQ=255)
                                                              aCCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCa                                                                                        >  1:661537/1‑150 (MQ=255)
                                                              aCCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCa                                                                                        <  1:980433/150‑1 (MQ=255)
                                                              aCCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCa                                                                                        >  1:2874790/1‑150 (MQ=255)
                                                              aCCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCa                                                                                        <  1:2994085/150‑1 (MQ=255)
                                                               cccAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAg                                                                                       <  1:265184/150‑1 (MQ=255)
                                                                  aGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGc                                                                                    >  1:1741579/1‑150 (MQ=255)
                                                                   gCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCt                                                                                   <  2:2097536/150‑1 (MQ=255)
                                                                      tGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCg                                                                                >  1:4152363/1‑150 (MQ=255)
                                                                          gTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGttt                                                                            <  1:1001530/150‑1 (MQ=255)
                                                                               tctcTTCGGTGATCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtcatca                                                                       <  2:290719/150‑1 (MQ=255)
                                                                               tctcTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtcatca                                                                       <  2:291046/150‑1 (MQ=255)
                                                                                ctctTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAg                                                                      >  1:3538801/1‑150 (MQ=255)
                                                                                    tCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTc                                                                  <  1:3774164/150‑1 (MQ=255)
                                                                                     cGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCt                                                                 <  2:4032396/150‑1 (MQ=255)
                                                                                     cGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCt                                                                 >  1:2326991/1‑150 (MQ=255)
                                                                                             gggTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCg                                                         <  1:3293295/150‑1 (MQ=255)
                                                                                                 tATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtg                                                     <  1:432342/150‑1 (MQ=255)
                                                                                                  aTGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgt                                                    <  2:1778486/150‑1 (MQ=255)
                                                                                                     gATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGc                                                 <  2:2302335/150‑1 (MQ=255)
                                                                                                     gATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGc                                                 >  1:3182353/1‑150 (MQ=255)
                                                                                                      aTGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGcc                                                >  2:1338393/1‑150 (MQ=255)
                                                                                                      aTGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGcc                                                <  1:3871207/150‑1 (MQ=255)
                                                                                                       tGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCg                                               >  1:2805069/1‑150 (MQ=255)
                                                                                                       tGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCg                                               >  1:2942543/1‑150 (MQ=255)
                                                                                                        gCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCgg                                              >  2:2649174/1‑150 (MQ=255)
                                                                                                         cGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGt                                             >  2:2329427/1‑150 (MQ=255)
                                                                                                            tgtTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAg                                          >  1:79148/1‑150 (MQ=255)
                                                                                                             gtTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAgg                                         >  1:4041081/1‑150 (MQ=255)
                                                                                                             gtTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAgg                                         >  2:2443835/1‑150 (MQ=255)
                                                                                                                cACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTc                                      <  1:870040/150‑1 (MQ=255)
                                                                                                                cACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTc                                      >  2:2851422/1‑150 (MQ=255)
                                                                                                                 aCTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCa                                     >  1:1593015/1‑150 (MQ=255)
                                                                                                                 aCTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCa                                     >  1:1593074/1‑150 (MQ=255)
                                                                                                                 aCTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGAGCACCGGATGTGGCCGGTCAGGGTCa                                     >  2:4073144/1‑150 (MQ=255)
                                                                                                                  cTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCAt                                    <  2:460108/150‑1 (MQ=255)
                                                                                                                        cTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTg                              >  2:3653028/1‑150 (MQ=255)
                                                                                                                          gATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAg                            <  2:636913/150‑1 (MQ=255)
                                                                                                                          gATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAg                            <  2:2898069/150‑1 (MQ=255)
                                                                                                                               tAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtg                       >  2:3948281/1‑150 (MQ=255)
                                                                                                                                    aTACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgagtga                  >  2:4080221/1‑150 (MQ=255)
                                                                                                                                    aTACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgagtga                  >  2:383178/1‑150 (MQ=255)
                                                                                                                                    aTACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgagtga                  >  1:764308/1‑150 (MQ=255)
                                                                                                                                      aCGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGAtt                <  2:141823/150‑1 (MQ=255)
                                                                                                                                          gggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcat            >  2:4146206/1‑150 (MQ=255)
                                                                                                                                          gggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcat            <  1:2360022/150‑1 (MQ=255)
                                                                                                                                           ggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc           >  2:1262167/1‑150 (MQ=255)
                                                                                                                                           ggATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc           >  2:2677112/1‑150 (MQ=255)
                                                                                                                                               atCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATcc       <  2:4192770/150‑1 (MQ=255)
                                                                                                                                               atCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATcc       >  2:746632/1‑150 (MQ=255)
                                                                                                                                                tCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCca      <  1:3399139/150‑1 (MQ=255)
                                                                                                                                                 ccATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcag     <  2:1829373/150‑1 (MQ=255)
                                                                                                                                                  cATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCcagc    >  1:1762018/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCTCGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCAGCAG  >  CP009273/4522852‑4523149

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: