Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,848,002 T→A 92.6% C102* (TGT→TGA)  fhlA → formate hydrogenlyase transcriptional activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,848,0020TA92.6% 585.4 / 29.1 190C102* (TGT→TGA) fhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (8/6);  new base A (92/84);  total (100/90)
Fisher's exact test for biased strand distribution p-value = 7.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 194 total aligned reads displayed.

CGCAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTT  >  CP009273/2847857‑2848145
                                                                                                                                                 |                                                                                                                                               
cgcAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTgaga                                                                                                                                               <  2:1475815/148‑1 (MQ=255)
cgcAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTgaga                                                                                                                                               <  2:2835285/148‑1 (MQ=255)
  cAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGaaa                                                                                                                                             >  1:1159541/1‑148 (MQ=255)
    gACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                           >  2:974459/1‑148 (MQ=255)
    gACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                           <  2:2763779/148‑1 (MQ=255)
     aCTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATGTTGAGAAACCt                                                                                                                                          >  2:2085721/1‑148 (MQ=255)
     aCTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTACATTGCAGTTACGAAGAATTTTGAGAAACCt                                                                                                                                          >  2:667319/1‑148 (MQ=255)
       tCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTgg                                                                                                                                        >  2:966461/1‑148 (MQ=255)
       tCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTgg                                                                                                                                        >  2:1737344/1‑148 (MQ=255)
        cAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGc                                                                                                                                       >  2:3235584/1‑148 (MQ=255)
          aCGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCg                                                                                                                                     <  2:150638/148‑1 (MQ=255)
           cGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGc                                                                                                                                    >  2:1167062/1‑148 (MQ=255)
             tGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAg                                                                                                                                  <  2:426309/148‑1 (MQ=255)
              gCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGc                                                                                                                                 <  2:1217771/148‑1 (MQ=255)
              gCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGc                                                                                                                                 >  1:1993513/1‑148 (MQ=255)
              gCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGc                                                                                                                                 >  2:837737/1‑148 (MQ=255)
               cGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCt                                                                                                                                <  2:1422554/148‑1 (MQ=255)
                 tCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTgg                                                                                                                              >  2:3112865/1‑148 (MQ=255)
                      ttaCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCa                                                                                                                         <  1:1752949/148‑1 (MQ=255)
                       taCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAg                                                                                                                        >  1:892371/1‑148 (MQ=255)
                        aCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAgg                                                                                                                       >  2:1784074/1‑148 (MQ=255)
                        aCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAgg                                                                                                                       >  2:324444/1‑148 (MQ=255)
                          gcgtcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTg                                                                                                                     <  2:1040246/148‑1 (MQ=255)
                          gcgtcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTg                                                                                                                     <  2:1040282/148‑1 (MQ=255)
                          gcgtcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTg                                                                                                                     <  1:608596/148‑1 (MQ=255)
                             tcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGc                                                                                                                  >  2:255591/1‑148 (MQ=255)
                                cgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCtat                                                                                                               >  1:1611374/1‑148 (MQ=255)
                                  tGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATAc                                                                                                             >  1:389702/1‑148 (MQ=255)
                                       aaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaaa                                                                                                        >  1:2790571/1‑148 (MQ=255)
                                         gACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAt                                                                                                      <  1:2368618/148‑1 (MQ=255)
                                         gACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAt                                                                                                      <  1:1883021/148‑1 (MQ=255)
                                         gACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAt                                                                                                      >  2:2779772/1‑148 (MQ=255)
                                         gACACCCCCATTAAAAATGAAGACGAAACTCTTATGTCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATAATCTGCATTGCAGTTACGAAGAAATTTGAGAAACCTGGCCGCAGCTGTACGCAGGTGGGCTATACCCAAAAt                                                                                                      >  2:2505251/1‑148 (MQ=255)
                                          acacCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                     <  1:1976840/148‑1 (MQ=255)
                                             cccccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTgg                                                                                                  >  1:673735/1‑148 (MQ=255)
                                             cccccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTgg                                                                                                  <  2:1714163/148‑1 (MQ=255)
                                              ccccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGt                                                                                                 >  2:2340005/1‑148 (MQ=255)
                                               cccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTc                                                                                                <  1:45644/148‑1 (MQ=255)
                                                ccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCa                                                                                               >  1:436453/1‑148 (MQ=255)
                                                  aTTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACt                                                                                             <  1:337339/148‑1 (MQ=255)
                                                    tAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAt                                                                                           >  1:2256675/1‑148 (MQ=255)
                                                      aaTATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTg                                                                                         >  2:1795881/1‑148 (MQ=255)
                                                        tatGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGcc                                                                                       >  1:249380/1‑148 (MQ=255)
                                                         atGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCt                                                                                      >  1:10769/1‑148 (MQ=255)
                                                          tGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTg                                                                                     >  1:2232061/1‑148 (MQ=255)
                                                          tGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTg                                                                                     >  1:2232040/1‑148 (MQ=255)
                                                           gAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGa                                                                                    <  2:2280386/148‑1 (MQ=255)
                                                             aGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATg                                                                                  >  1:2737256/1‑148 (MQ=255)
                                                             aGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATg                                                                                  <  1:849190/148‑1 (MQ=255)
                                                               aCGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGcc                                                                                >  2:524764/1‑148 (MQ=255)
                                                               aCGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGcc                                                                                >  2:294903/1‑148 (MQ=255)
                                                                cGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCa                                                                               >  1:100743/1‑148 (MQ=255)
                                                                  aaaCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACt                                                                             >  1:1863080/1‑148 (MQ=255)
                                                                   aaCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTg                                                                            <  2:380068/148‑1 (MQ=255)
                                                                    aCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTgg                                                                           >  1:514443/1‑148 (MQ=255)
                                                                    aCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTgg                                                                           >  1:18172/1‑148 (MQ=255)
                                                                           tGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGa                                                                    <  1:2425604/148‑1 (MQ=255)
                                                                            ggCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGaa                                                                   <  2:1298293/148‑1 (MQ=255)
                                                                             gCACACGGTCCGGTACGCAGCAGTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTTGCGGCGGacg                                                                  >  2:2222804/1‑146 (MQ=255)
                                                                                cacGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGc                                                               <  2:2794323/148‑1 (MQ=255)
                                                                                      ccGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttt                                                         <  2:1136508/148‑1 (MQ=255)
                                                                                      ccGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttt                                                         <  2:3162878/148‑1 (MQ=255)
                                                                                      ccGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttt                                                         <  1:785582/148‑1 (MQ=255)
                                                                                          tACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtgg                                                     <  1:2087010/148‑1 (MQ=255)
                                                                                              cagcaTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGc                                                 >  1:2580933/1‑148 (MQ=255)
                                                                                              cagcaTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGc                                                 <  2:2693063/148‑1 (MQ=255)
                                                                                               agcaTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCt                                                <  2:2583721/148‑1 (MQ=255)
                                                                                                gcaTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtg                                               >  1:213978/1‑148 (MQ=255)
                                                                                                   ttttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGa                                            >  1:103769/1‑148 (MQ=255)
                                                                                                    tttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGaa                                           >  2:2028840/1‑148 (MQ=255)
                                                                                                    tttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGaa                                           <  1:2170120/148‑1 (MQ=255)
                                                                                                             cccTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCAGATTTTTGGTGGCTGTGAATTTATTCGt                                  >  2:2387553/1‑148 (MQ=255)
                                                                                                              ccTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGtt                                 <  1:939291/148‑1 (MQ=255)
                                                                                                              ccTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGtt                                 <  1:2016815/148‑1 (MQ=255)
                                                                                                               cTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTa                                <  2:1611374/148‑1 (MQ=255)
                                                                                                                 gATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATg                              <  2:2377874/148‑1 (MQ=255)
                                                                                                                  aTACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGa                             <  2:1998262/148‑1 (MQ=255)
                                                                                                                  aTACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGa                             >  1:1251622/1‑148 (MQ=255)
                                                                                                                   tACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGAc                            >  1:2806822/1‑148 (MQ=255)
                                                                                                                     cGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGa                          <  2:840238/148‑1 (MQ=255)
                                                                                                                      gCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGAt                         <  1:2977656/148‑1 (MQ=255)
                                                                                                                      gCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGAt                         <  1:1331077/148‑1 (MQ=255)
                                                                                                                       cTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATc                        >  1:171452/1‑148 (MQ=255)
                                                                                                                          cATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcc                     <  1:204468/148‑1 (MQ=255)
                                                                                                                           aTTTCAGTTACTAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCTATATTGCCTGATGCCACTGGCGGCGGAAGGGCTTATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGccc                    <  1:1987632/148‑1 (MQ=255)
                                                                                                                           aTTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGccc                    <  1:1313785/148‑1 (MQ=255)
                                                                                                                            ttGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcccc                   <  1:586934/148‑1 (MQ=255)
                                                                                                                                 gTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAg              <  2:325852/148‑1 (MQ=255)
                                                                                                                                 gTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAg              <  2:1993512/148‑1 (MQ=255)
                                                                                                                                 gTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAg              <  2:1993513/148‑1 (MQ=255)
                                                                                                                                   tACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCAGATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCg            >  2:400576/1‑148 (MQ=255)
                                                                                                                                      gaagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaa         <  1:1061980/148‑1 (MQ=255)
                                                                                                                                      gaagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaa         <  2:931818/148‑1 (MQ=255)
                                                                                                                                        agaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaaaa       >  1:1578313/1‑148 (MQ=255)
                                                                                                                                        agaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaaaa       >  1:749296/1‑148 (MQ=255)
                                                                                                                                         gaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAag      >  2:1359305/1‑148 (MQ=255)
                                                                                                                                         gaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAag      >  2:1472183/1‑148 (MQ=255)
                                                                                                                                            ttttGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGt   <  1:3224093/148‑1 (MQ=255)
                                                                                                                                            ttttGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGt   <  1:2189097/148‑1 (MQ=255)
                                                                                                                                             tttGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGGGtt  <  1:1491263/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CGCAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTT  >  CP009273/2847857‑2848145

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: ‑