Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,984,342 G→T 16.6% intergenic (‑216/‑171) hemC ← / → cyaA hydroxymethylbilane synthase/adenylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,984,3420GT16.6% 948.6 / 110.6 298intergenic (‑216/‑171)hemC/cyaAhydroxymethylbilane synthase/adenylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (117/131);  new base T (15/35);  total (132/166)
Fisher's exact test for biased strand distribution p-value = 2.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 304 total aligned reads displayed.

GGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGG  >  CP009273/3984193‑3984491
                                                                                                                                                     |                                                                                                                                                     
      aCGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTa                                                                                                                                                 <  2:2563089/150‑1 (MQ=255)
       cGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAACCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTaa                                                                                                                                                <  1:3685626/150‑1 (MQ=255)
       cGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTaa                                                                                                                                                <  1:1005992/150‑1 (MQ=255)
        gTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTaaa                                                                                                                                               >  1:668309/1‑150 (MQ=255)
         tAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAAt                                                                                                                                              >  1:1813115/1‑150 (MQ=255)
           gACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTg                                                                                                                                            <  2:4023100/150‑1 (MQ=255)
            aCGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGa                                                                                                                                           >  1:1629185/1‑150 (MQ=255)
            aCGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGa                                                                                                                                           <  2:2565571/150‑1 (MQ=255)
              gCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATc                                                                                                                                         <  2:2847785/150‑1 (MQ=255)
               ccTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCa                                                                                                                                        <  2:2229244/150‑1 (MQ=255)
               ccTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCa                                                                                                                                        <  2:1482249/150‑1 (MQ=255)
               ccTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCa                                                                                                                                        <  2:3211656/150‑1 (MQ=255)
                   cgcAAACCGTACAATGAGGTCTGGCAGTGGATCCTGACCGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGtt                                                                                                                                    <  2:167533/150‑1 (MQ=255)
                     cAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGGAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGtttt                                                                                                                                  <  2:3770086/150‑1 (MQ=255)
                      aaaCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTa                                                                                                                                 >  2:2649629/1‑150 (MQ=255)
                        aCCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGa                                                                                                                               <  2:2986054/150‑1 (MQ=255)
                           gTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCa                                                                                                                            <  2:1394678/150‑1 (MQ=255)
                            tAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAt                                                                                                                           >  1:492131/1‑150 (MQ=255)
                             aaaaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCAtt                                                                                                                          <  2:2411457/150‑1 (MQ=255)
                              aaaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAttt                                                                                                                         <  1:2721846/150‑1 (MQ=255)
                              aaaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAttt                                                                                                                         <  1:1473735/150‑1 (MQ=255)
                              aaaTGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAttt                                                                                                                         <  1:1903580/150‑1 (MQ=255)
                                 tGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCAtttttt                                                                                                                      >  2:140610/1‑150 (MQ=255)
                                  gAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtc                                                                                                                     <  2:3082415/150‑1 (MQ=255)
                                  gAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTtc                                                                                                                     >  1:373724/1‑150 (MQ=255)
                                       cTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTcgtcgt                                                                                                                >  1:165499/1‑150 (MQ=255)
                                        tGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTg                                                                                                               >  2:3267627/1‑150 (MQ=255)
                                        tGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTg                                                                                                               >  1:2108101/1‑150 (MQ=255)
                                           cAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGaaa                                                                                                            >  2:523702/1‑150 (MQ=255)
                                               ggATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaa                                                                                                        <  1:2674367/150‑1 (MQ=255)
                                                    cTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTaaaaaaa                                                                                                   >  1:1298457/1‑150 (MQ=255)
                                                    cTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTAGTATGTAGCGCATCTTTCTTTACGGTCAAACAGCAATGTGTTAAATTGATCACGTTTTATACCATTTTTTCGTCGTGAAACTaaaaaaa                                                                                                   >  1:1172688/1‑150 (MQ=255)
                                                     tGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAc                                                                                                  >  1:1652035/1‑150 (MQ=255)
                                                     tGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAAc                                                                                                  <  1:2204325/150‑1 (MQ=255)
                                                          ggCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGgc                                                                                             <  1:3868925/150‑1 (MQ=255)
                                                             gTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGcgcg                                                                                          <  1:367813/150‑1 (MQ=255)
                                                               ttCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGaa                                                                                        <  2:840758/150‑1 (MQ=255)
                                                                tCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGaaa                                                                                       <  1:3471344/150‑1 (MQ=255)
                                                                  aCGCCGTTGTAATAAGGACTTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAg                                                                                     >  2:557470/1‑150 (MQ=255)
                                                                  aCGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAg                                                                                     <  1:3840355/150‑1 (MQ=255)
                                                                  aCGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAg                                                                                     >  2:516902/1‑150 (MQ=255)
                                                                     ccGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTgt                                                                                  >  1:4635175/1‑149 (MQ=255)
                                                                         tGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAAc                                                                              <  1:2298952/150‑1 (MQ=255)
                                                                          gTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACg                                                                             >  2:3245696/1‑150 (MQ=255)
                                                                             ataaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGtt                                                                          >  1:1724599/1‑150 (MQ=255)
                                                                             ataaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGtt                                                                          >  2:4016028/1‑150 (MQ=255)
                                                                              taaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTa                                                                         <  2:2493242/150‑1 (MQ=255)
                                                                              taaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTa                                                                         <  2:3920680/150‑1 (MQ=255)
                                                                              taaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTGTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGGAACGGTTa                                                                         >  2:3533420/1‑150 (MQ=255)
                                                                               aaGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTAc                                                                        <  2:1586979/150‑1 (MQ=255)
                                                                                acgAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTAcc                                                                       >  1:1149623/3‑150 (MQ=255)
                                                                                aGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTAcc                                                                       <  1:3083338/150‑1 (MQ=255)
                                                                                  gAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCtt                                                                     <  1:4235842/150‑1 (MQ=255)
                                                                                    aTTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTg                                                                   <  2:3804725/150‑1 (MQ=255)
                                                                                     tttACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGa                                                                  <  2:1750771/150‑1 (MQ=255)
                                                                                         cAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTACCGGTTACCTTTGACATa                                                              >  2:732843/1‑150 (MQ=255)
                                                                                          agagAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATAc                                                             >  1:4442667/1‑150 (MQ=255)
                                                                                              aaTAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGaaa                                                         >  1:3245695/1‑150 (MQ=255)
                                                                                                 aaaCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAatat                                                      <  2:2108101/150‑1 (MQ=255)
                                                                                                  aaCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATc                                                     <  1:3144745/150‑1 (MQ=255)
                                                                                                   aCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATcc                                                    <  2:1105690/150‑1 (MQ=255)
                                                                                                   aCGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATcc                                                    <  1:2551401/150‑1 (MQ=255)
                                                                                                     ggTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCg                                                  <  1:3214444/150‑1 (MQ=255)
                                                                                                     ggTGCTACACTTGTATGTAGCGCAACTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCg                                                  <  2:4263567/150‑1 (MQ=255)
                                                                                                        gCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAAt                                               <  2:3401066/150‑1 (MQ=255)
                                                                                                         cTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATg                                              >  1:3586653/1‑150 (MQ=255)
                                                                                                              cTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcgcg                                         >  2:3890189/1‑150 (MQ=255)
                                                                                                              cTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcgcg                                         >  2:12693/1‑150 (MQ=255)
                                                                                                              cTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcgcg                                         <  1:4166107/150‑1 (MQ=255)
                                                                                                              cTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCcgcg                                         <  2:1934017/150‑1 (MQ=255)
                                                                                                                  tatgtaGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGtt                                     >  1:2010921/1‑150 (MQ=255)
                                                                                                                   atgtaGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTa                                    >  1:4230335/1‑150 (MQ=255)
                                                                                                                         cgcATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTg                              >  1:902142/1‑150 (MQ=255)
                                                                                                                           cATCTTTCTTTACGTTCAATCAGCTAGGTGTTAAATTGATCACGTTTTTGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTATCGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGAt                            <  1:2849951/150‑1 (MQ=255)
                                                                                                                           cATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGAt                            <  1:4065486/150‑1 (MQ=255)
                                                                                                                           cATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGAt                            <  1:95737/150‑1 (MQ=255)
                                                                                                                           cATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGAt                            <  2:1813115/150‑1 (MQ=255)
                                                                                                                            aTCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATg                           >  1:297067/1‑150 (MQ=255)
                                                                                                                              ctttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGtt                         >  2:4330910/1‑150 (MQ=255)
                                                                                                                              ctttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGtt                         >  2:4474643/1‑150 (MQ=255)
                                                                                                                              ctttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGtt                         >  1:1880811/1‑150 (MQ=255)
                                                                                                                              ctttcttgACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGtt                         <  2:3150850/150‑1 (MQ=255)
                                                                                                                               tttctttACGGTCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTg                        <  1:924419/150‑1 (MQ=255)
                                                                                                                               tttctttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTg                        >  1:2390293/1‑150 (MQ=255)
                                                                                                                                   tttACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCgg                    >  1:543623/1‑150 (MQ=255)
                                                                                                                                      aCGGTCAATCAGCTATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAGTAGTGGTATCGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGTAAt                 <  1:2784429/150‑1 (MQ=255)
                                                                                                                                      aCGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAAt                 >  2:2825780/1‑150 (MQ=255)
                                                                                                                                        ggTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATca               >  2:3281776/1‑150 (MQ=255)
                                                                                                                                        ggTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATca               <  1:3303178/150‑1 (MQ=255)
                                                                                                                                        ggTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATca               >  2:4428005/1‑150 (MQ=255)
                                                                                                                                          tCAATCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATcaca             <  2:2557832/150‑1 (MQ=255)
                                                                                                                                             aTCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTc          <  2:2065677/150‑1 (MQ=255)
                                                                                                                                             aTCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTc          >  1:4242695/1‑150 (MQ=255)
                                                                                                                                              tCAGCAATGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa         <  2:4145762/150‑1 (MQ=255)
                                                                                                                                              tCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCa         >  2:584328/1‑150 (MQ=255)
                                                                                                                                               cAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCAt        >  1:4584412/1‑150 (MQ=255)
                                                                                                                                                  cAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGAc     <  1:795930/150‑1 (MQ=255)
                                                                                                                                                   aaGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACg    <  1:2215522/150‑1 (MQ=255)
                                                                                                                                                   aaGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACg    >  1:602101/1‑150 (MQ=255)
                                                                                                                                                     ggTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACggg  >  2:4560048/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGAAACGTAGACGCCTGCGCAAACCGTAAAATGAGGTCTGGCAGTGGATCCTGACAGGCGTTTCACGCCGTTGTAATAAGGAATTTACAGAGAATAAACGGTGCTACACTTGTATGTAGCGCATCTTTCTTTACGGTCAATCAGCAAGGTGTTAAATTGATCACGTTTTAGACCATTTTTTCGTCGTGAAACTAAAAAAACCAGGCGCGAAAAGTGGTAACGGTTACCTTTGACATACGAAATATCCCGAATGCCGCGTGTTACCGTTGATGTTGGCGGAATCACAGTCATGACGGG  >  CP009273/3984193‑3984491

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: