Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,523,001 T→A 100% 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,0010TA95.3% 513.7 / 12.7 169intergenic (‑1/+55)yjhP/yjhQputative methyltransferase/putative acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base A (71/90);  total (75/94)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 172 total aligned reads displayed.

CGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCA  >  CP009273/4522856‑4523145
                                                                                                                                                 |                                                                                                                                                
cGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCataat                                                                                                                                                <  1:2032883/148‑1 (MQ=255)
 gCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCataatt                                                                                                                                               >  1:1140516/1‑148 (MQ=255)
   ccGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTa                                                                                                                                             <  2:2437758/148‑1 (MQ=255)
      aGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCAtaatttaatt                                                                                                                                          >  1:1076855/1‑148 (MQ=255)
        ccGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTcc                                                                                                                                        >  2:3438961/1‑148 (MQ=255)
          gCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTg                                                                                                                                      <  1:1076651/148‑1 (MQ=255)
           cTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGt                                                                                                                                     <  2:2461547/148‑1 (MQ=255)
           cTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGt                                                                                                                                     <  1:32640/148‑1 (MQ=255)
             gCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTa                                                                                                                                   >  1:2586585/1‑148 (MQ=255)
               cGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTAtt                                                                                                                                 <  2:2198607/148‑1 (MQ=255)
                  ggTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTggg                                                                                                                              <  2:2976115/148‑1 (MQ=255)
                      aaGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTaa                                                                                                                          >  1:3082978/1‑148 (MQ=255)
                      aaGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTaa                                                                                                                          <  1:2630116/148‑1 (MQ=255)
                         aaTGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGcc                                                                                                                       <  1:3192626/148‑1 (MQ=255)
                         aaTGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGcc                                                                                                                       >  2:1919619/1‑148 (MQ=255)
                             cGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaa                                                                                                                   >  1:891006/1‑148 (MQ=255)
                              gggTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaaa                                                                                                                  >  2:1892871/1‑148 (MQ=255)
                               ggTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaaaa                                                                                                                 >  2:3493400/1‑148 (MQ=255)
                               ggTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaaaa                                                                                                                 >  1:3024094/1‑148 (MQ=255)
                                gTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTGTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTaaaaa                                                                                                                <  1:2915346/148‑1 (MQ=255)
                                 tGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAg                                                                                                               <  1:3385114/148‑1 (MQ=255)
                                  gCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATACTTTAATTCCTGTTAGTGGGTTAAGCCTGAAAAAgg                                                                                                              >  2:1628218/1‑148 (MQ=255)
                                  gCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAgg                                                                                                              >  2:602173/1‑148 (MQ=255)
                                     cGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGaca                                                                                                           <  1:852155/148‑1 (MQ=255)
                                       gCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGAcaca                                                                                                         <  2:1905272/148‑1 (MQ=255)
                                        cTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCGGTTATTGGGTTAAGCCTTAAAAAGGACacaa                                                                                                        >  2:3031630/1‑148 (MQ=255)
                                         ttCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaac                                                                                                       <  2:1353347/148‑1 (MQ=255)
                                          tCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACacaaca                                                                                                      >  2:1117318/1‑148 (MQ=255)
                                            aTGCGCAGCACGCGACCCATCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATTTTAAATATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAg                                                                                                    <  1:2587989/148‑1 (MQ=255)
                                             tGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGa                                                                                                   <  1:1178971/148‑1 (MQ=255)
                                             tGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGa                                                                                                   <  1:73948/148‑1 (MQ=255)
                                                gcagcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTg                                                                                                <  1:2122574/148‑1 (MQ=255)
                                                 cagcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGt                                                                                               <  1:2953674/148‑1 (MQ=255)
                                                 cagcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGt                                                                                               >  2:3405214/1‑148 (MQ=255)
                                                 cagcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGt                                                                                               <  1:1884681/148‑1 (MQ=255)
                                                 cagcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGGGGGATATCCATAATTTAATTCATGTTATTGGGTTAAGCCTTAAAAAGGACACACCAAGATTGt                                                                                               >  1:3339238/1‑148 (MQ=255)
                                                  agcaCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGtt                                                                                              <  1:2224122/148‑1 (MQ=255)
                                                    caCGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTAt                                                                                            >  1:3483618/1‑148 (MQ=255)
                                                          aCCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGt                                                                                      >  1:909062/1‑148 (MQ=255)
                                                          aCCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGt                                                                                      <  1:2956476/148‑1 (MQ=255)
                                                           cccAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTc                                                                                     >  2:1145083/1‑148 (MQ=255)
                                                               gCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTg                                                                                 >  2:2348367/1‑148 (MQ=255)
                                                               gCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTg                                                                                 <  1:3043744/148‑1 (MQ=255)
                                                                cGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGAGTGTTATTCACGTCAGTGc                                                                                >  2:2057495/1‑148 (MQ=255)
                                                                  tGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTc                                                                              <  1:1449619/148‑1 (MQ=255)
                                                                     cGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCgg                                                                           <  2:2586585/148‑1 (MQ=255)
                                                                      gTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGt                                                                          <  1:667072/148‑1 (MQ=255)
                                                                       tACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGtt                                                                         <  2:584312/148‑1 (MQ=255)
                                                                        aCTTCTCTTCGGTGAACGGGTTATGGATGCTGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTTTTATAGTGTTAAGCCTTAAAAACGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGttt                                                                        >  2:3133736/1‑148 (MQ=255)
                                                                          ttCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTtca                                                                      <  2:486262/148‑1 (MQ=255)
                                                                              ctTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAg                                                                  <  1:2899921/148‑1 (MQ=255)
                                                                                   gTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCt                                                             <  1:2814093/148‑1 (MQ=255)
                                                                                   gTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATAAGGGTCt                                                             >  1:2269592/1‑148 (MQ=255)
                                                                                   gTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACAAAACAATATTGTTATTCAAGTAAGTGCTCCTGATTCATCAGCGTCt                                                             >  1:386732/1‑148 (MQ=255)
                                                                                       aCGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTgcgc                                                         >  1:232234/1‑148 (MQ=255)
                                                                                            ttATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCgg                                                    >  2:1831984/1‑148 (MQ=255)
                                                                                             tATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGa                                                   >  2:678963/1‑148 (MQ=255)
                                                                                              aTGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAt                                                  <  2:1076855/148‑1 (MQ=255)
                                                                                               tGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtg                                                 >  2:3324745/1‑148 (MQ=255)
                                                                                                ggATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgt                                                <  2:252878/148‑1 (MQ=255)
                                                                                                 gATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgtg                                               <  1:487062/148‑1 (MQ=255)
                                                                                                 gATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGAtgtg                                               <  1:3369218/148‑1 (MQ=255)
                                                                                                     cGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCg                                           >  2:1368422/1‑148 (MQ=255)
                                                                                                      ggTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCgg                                          <  1:381883/148‑1 (MQ=255)
                                                                                                               tttCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCa                                 >  2:2995975/1‑148 (MQ=255)
                                                                                                                ttCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCAt                                <  1:171638/148‑1 (MQ=255)
                                                                                                                ttCGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCAt                                <  2:143734/148‑1 (MQ=255)
                                                                                                                  cGCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATgg                              <  2:136377/148‑1 (MQ=255)
                                                                                                                   gCTGATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTgg                                              >  1:3148108/1‑131 (MQ=255)
                                                                                                                      gATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTg                          >  2:2841001/1‑148 (MQ=255)
                                                                                                                      gATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTg                          >  2:2798063/1‑148 (MQ=255)
                                                                                                                      gATGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTg                          >  2:2799018/1‑148 (MQ=255)
                                                                                                                       aTGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGa                         >  2:2879978/1‑148 (MQ=255)
                                                                                                                       aTGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGa                         <  2:2838044/148‑1 (MQ=255)
                                                                                                                        tGGTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAg                        <  1:2451284/148‑1 (MQ=255)
                                                                                                                         ggTAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGc                       >  2:2044012/1‑148 (MQ=255)
                                                                                                                           tAAAAATACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCgg                     <  2:891006/148‑1 (MQ=255)
                                                                                                                            aaaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgt                    >  1:922462/1‑148 (MQ=255)
                                                                                                                            aaaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgt                    <  1:467085/148‑1 (MQ=255)
                                                                                                                             aaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACCACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtg                   <  2:802988/148‑1 (MQ=255)
                                                                                                                             aaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtg                   <  2:1106090/148‑1 (MQ=255)
                                                                                                                             aaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtg                   <  2:802998/148‑1 (MQ=255)
                                                                                                                             aaaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtg                   <  2:802851/148‑1 (MQ=255)
                                                                                                                              aaaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtga                  >  1:707128/1‑148 (MQ=255)
                                                                                                                               aaTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgag                 <  1:1202140/148‑1 (MQ=255)
                                                                                                                                aTACGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgagt                <  1:1919619/148‑1 (MQ=255)
                                                                                                                                  aCGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgagtga              >  1:1352093/1‑148 (MQ=255)
                                                                                                                                  aCGTGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGgtgagtga              <  1:1892871/148‑1 (MQ=255)
                                                                                                                                     tGGGATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTg           <  1:593831/148‑1 (MQ=255)
                                                                                                                                      gggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGc          >  2:1510628/1‑148 (MQ=255)
                                                                                                                                      gggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGc          >  2:370441/1‑148 (MQ=255)
                                                                                                                                      gggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGc          >  2:370609/1‑148 (MQ=255)
                                                                                                                                      gggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGc          >  2:1510458/1‑148 (MQ=255)
                                                                                                                                      gggATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGc          >  1:2352765/1‑148 (MQ=255)
                                                                                                                                        gATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcat        <  1:1320512/148‑1 (MQ=255)
                                                                                                                                        gATATCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcat        <  1:240968/148‑1 (MQ=255)
                                                                                                                                         atatCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc       >  2:2802656/1‑148 (MQ=255)
                                                                                                                                         atatCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGcatc       >  2:3257463/1‑148 (MQ=255)
                                                                                                                                            tCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCatcatc    >  1:504889/1‑148 (MQ=255)
                                                                                                                                            tCCATAATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCatcatc    <  1:2261116/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCCCGAGCCGCTGCCGAGGTCAAGAATGCGGGTGCCCGGCTTCATGCGCAGCACGCGACCCAGCGTGGCGTACTTCTCTTCGGTGAACGGGTTATGGATGCGGTGTTCACTTTCGCTGATGGTAAAAATACGTGGGATATCCATTATTTAATTCCTGTTATTGGGTTAAGCCTTAAAAAGGACACAACAAGATTGTTATTCACGTCAGTGCTCCGGTTTCATCAGGGTCTGCGCACCGGATGTGGCCGGTCAGGGTCATGGGCTGAGCGGTGAGTGATTGCATCATCCA  >  CP009273/4522856‑4523145

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: