Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,848,002 T→A 100% 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,0020TA97.4% 634.6 / 4.6 194C102* (TGT→TGA) fhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base A (100/89);  total (102/92)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 195 total aligned reads displayed.

CGCAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTTC  >  CP009273/2847857‑2848146
                                                                                                                                                 |                                                                                                                                                
 gcAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGaa                                                                                                                                               <  2:1175058/148‑1 (MQ=255)
  cAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGaaa                                                                                                                                              >  1:400538/1‑148 (MQ=255)
  cAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGaaa                                                                                                                                              >  1:756485/1‑148 (MQ=255)
   aGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAc                                                                                                                                             >  2:1625672/1‑148 (MQ=255)
   aGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAc                                                                                                                                             <  2:290438/148‑1 (MQ=255)
    gACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                            >  2:2110171/1‑148 (MQ=255)
    gACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                            <  2:240785/148‑1 (MQ=255)
    gACTCAACGTGCGTCTTATAACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAACACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                            <  1:643647/148‑1 (MQ=255)
      cTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTg                                                                                                                                          <  2:2160190/148‑1 (MQ=255)
      cTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTg                                                                                                                                          <  1:2882644/148‑1 (MQ=255)
       tCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTgg                                                                                                                                         <  1:3503234/148‑1 (MQ=255)
        cAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGc                                                                                                                                        >  2:686146/1‑148 (MQ=255)
         aaCGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGcc                                                                                                                                       >  1:1210835/1‑148 (MQ=255)
             tGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAg                                                                                                                                   <  1:1017834/148‑1 (MQ=255)
              gCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGc                                                                                                                                  >  1:136803/1‑148 (MQ=255)
                  cTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGa                                                                                                                              >  2:544974/1‑148 (MQ=255)
                  cTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGa                                                                                                                              >  1:485721/1‑148 (MQ=255)
                       taCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAg                                                                                                                         <  1:1986639/148‑1 (MQ=255)
                            gtcgcgtGAAAAAGACACCCCCATTAACTATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTggg                                                                                                                    <  1:2731951/148‑1 (MQ=255)
                             tcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGc                                                                                                                   <  2:1846076/148‑1 (MQ=255)
                                cgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTAGCCGCAGCTGGACGCAGGTGGGCtat                                                                                                                <  2:3186011/148‑1 (MQ=255)
                                   gAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATAcc                                                                                                             <  2:926777/148‑1 (MQ=255)
                                     aaaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCa                                                                                                           >  1:2980518/1‑148 (MQ=255)
                                      aaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATCCGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaa                                                                                                          >  2:2729958/1‑148 (MQ=255)
                                         gACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAt                                                                                                       >  1:2435787/1‑148 (MQ=255)
                                         gACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAt                                                                                                       >  1:2435852/1‑148 (MQ=255)
                                          acacCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                      >  1:104700/1‑148 (MQ=255)
                                          acacCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                      >  1:686876/1‑148 (MQ=255)
                                          acacCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                      >  1:3276112/1‑148 (MQ=255)
                                           cacCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAttt                                                                                                     <  2:2841423/148‑1 (MQ=255)
                                               cccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTc                                                                                                 >  1:996211/1‑148 (MQ=255)
                                               cccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTc                                                                                                 >  1:3337689/1‑148 (MQ=255)
                                               cccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTc                                                                                                 >  1:3337662/1‑148 (MQ=255)
                                               cccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTc                                                                                                 <  2:81507/148‑1 (MQ=255)
                                                ccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCa                                                                                                >  2:2782404/1‑148 (MQ=255)
                                                    tAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAt                                                                                            >  1:2865471/1‑148 (MQ=255)
                                                    tAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAt                                                                                            <  2:466499/148‑1 (MQ=255)
                                                         atGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCt                                                                                       >  1:293660/1‑148 (MQ=255)
                                                         atGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCt                                                                                       >  1:1944463/1‑148 (MQ=255)
                                                           gAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGa                                                                                     <  1:898428/148‑1 (MQ=255)
                                                           gAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGa                                                                                     <  1:2686912/148‑1 (MQ=255)
                                                             aGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGAGACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATg                                                                                   <  2:3421177/148‑1 (MQ=255)
                                                              gACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGc                                                                                  >  1:1444250/1‑148 (MQ=255)
                                                                cGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCa                                                                                <  1:3162631/148‑1 (MQ=255)
                                                                   aaCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTg                                                                             <  2:2379936/148‑1 (MQ=255)
                                                                    aCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTgg                                                                            >  1:2546934/1‑148 (MQ=255)
                                                                    aCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTgg                                                                            <  2:2493082/148‑1 (MQ=255)
                                                                       gTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggcgg                                                                         >  1:69722/1‑148 (MQ=255)
                                                                       gTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggcgg                                                                         <  2:704978/148‑1 (MQ=255)
                                                                       gTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggcgg                                                                         >  1:1977446/1‑148 (MQ=255)
                                                                          cTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGcggcgg                                                                      >  2:177683/1‑148 (MQ=255)
                                                                                cacGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGc                                                                <  2:3433972/148‑1 (MQ=255)
                                                                                cacGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGc                                                                <  2:3318417/148‑1 (MQ=255)
                                                                                  cGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCat                                                              <  2:1050957/148‑1 (MQ=255)
                                                                                   ggTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCata                                                             <  2:1221901/148‑1 (MQ=255)
                                                                                     tCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAtt                                                           <  1:2583979/148‑1 (MQ=255)
                                                                                      ccGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttt                                                          <  1:2205732/148‑1 (MQ=255)
                                                                                        ggTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttttt                                                        <  1:1625500/148‑1 (MQ=255)
                                                                                         gTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtg                                                       <  1:3273369/148‑1 (MQ=255)
                                                                                               agcaTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCt                                                 <  2:1276334/148‑1 (MQ=255)
                                                                                                gcaTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtg                                                <  2:2435852/148‑1 (MQ=255)
                                                                                                 caTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtgt                                               >  1:2809273/1‑148 (MQ=255)
                                                                                                  aTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtgtg                                              <  1:2287103/148‑1 (MQ=255)
                                                                                                   ttttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCACAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCCTATTTTTGGTGGCTGTGa                                             >  2:2017102/1‑148 (MQ=255)
                                                                                                    tttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGaa                                            >  1:3133248/1‑148 (MQ=255)
                                                                                                     ttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAAt                                           <  1:1076004/148‑1 (MQ=255)
                                                                                                     ttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAAt                                           <  1:2428198/148‑1 (MQ=255)
                                                                                                      tGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAAtt                                          <  2:1771105/148‑1 (MQ=255)
                                                                                                         cgcgCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTAt                                       <  2:1137676/148‑1 (MQ=255)
                                                                                                          gcgcCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTAtt                                      <  2:254320/148‑1 (MQ=255)
                                                                                                           cgcCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTc                                     >  2:1025450/1‑148 (MQ=255)
                                                                                                             cccTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGt                                   <  2:1018666/148‑1 (MQ=255)
                                                                                                               cTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAGTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTa                                 <  2:3272189/148‑1 (MQ=255)
                                                                                                                tGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTAt                                >  2:3376164/1‑148 (MQ=255)
                                                                                                                tGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTAt                                >  1:2937207/1‑148 (MQ=255)
                                                                                                                 gATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATg                               <  2:130427/148‑1 (MQ=255)
                                                                                                                  aTACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGa                              <  1:1911868/148‑1 (MQ=255)
                                                                                                                    aCGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACg                            >  2:361351/1‑148 (MQ=255)
                                                                                                                    aCGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACg                            >  1:2237206/1‑148 (MQ=255)
                                                                                                                        tGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCg                        <  2:668948/148‑1 (MQ=255)
                                                                                                                         gCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGc                       >  1:3357150/1‑148 (MQ=255)
                                                                                                                          cATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcc                      >  2:810489/1‑148 (MQ=255)
                                                                                                                          cATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcc                      >  2:1552755/1‑148 (MQ=255)
                                                                                                                            ttGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcccc                    >  1:2903700/1‑148 (MQ=255)
                                                                                                                            ttGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcccc                    <  1:867412/148‑1 (MQ=255)
                                                                                                                            ttGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcccc                    <  1:2647533/148‑1 (MQ=255)
                                                                                                                            ttGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcccc                    >  1:1494859/1‑148 (MQ=255)
                                                                                                                               cAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTgg                 >  2:1332615/1‑148 (MQ=255)
                                                                                                                               cAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTgg                 <  2:1587478/148‑1 (MQ=255)
                                                                                                                                  ttACGAAGAATTTTTAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGc              >  1:2587995/1‑148 (MQ=255)
                                                                                                                                  ttACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGc              >  2:1581754/1‑148 (MQ=255)
                                                                                                                                     cGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaa           <  1:2471413/148‑1 (MQ=255)
                                                                                                                                      gaagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaa          <  2:1091841/148‑1 (MQ=255)
                                                                                                                                      gaagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGATATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaa          <  2:2679522/148‑1 (MQ=255)
                                                                                                                                        agaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaaaa        >  1:360178/1‑148 (MQ=255)
                                                                                                                                          aaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAaga      >  1:2388381/1‑148 (MQ=255)
                                                                                                                                          aaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAaga      >  1:1199720/1‑148 (MQ=255)
                                                                                                                                          aaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAaga      <  1:1898397/148‑1 (MQ=255)
                                                                                                                                            ttttGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGt    >  2:338864/1‑148 (MQ=255)
                                                                                                                                              ttGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTTc  <  1:1348294/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTTC  >  CP009273/2847857‑2848146

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