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% 732.4 / 6.5 226C102* (TGT→TGA) fhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (1/5);  new base A (109/111);  total (110/116)
Fisher's exact test for biased strand distribution p-value = 2.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 229 total aligned reads displayed.

CGCAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTT  >  CP009273/2847857‑2848145
                                                                                                                                                 |                                                                                                                                               
cgcAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTgaga                                                                                                                                               >  1:2953093/1‑148 (MQ=255)
 gcAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGaa                                                                                                                                              <  1:486892/148‑1 (MQ=255)
  cAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGaaa                                                                                                                                             <  1:1048708/148‑1 (MQ=255)
  cAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGaaa                                                                                                                                             <  1:28259/148‑1 (MQ=255)
   aGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAc                                                                                                                                            <  1:267590/148‑1 (MQ=255)
    gACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                           <  1:2740793/148‑1 (MQ=255)
     aCTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCt                                                                                                                                          <  2:2727440/148‑1 (MQ=255)
     aCTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCt                                                                                                                                          <  1:2782813/148‑1 (MQ=255)
     aCTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCt                                                                                                                                          >  1:2361595/1‑148 (MQ=255)
      cTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTg                                                                                                                                         <  1:2717002/148‑1 (MQ=255)
        cAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGc                                                                                                                                       <  2:1420481/148‑1 (MQ=255)
             tGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAg                                                                                                                                  >  2:2892647/1‑148 (MQ=255)
             tGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAg                                                                                                                                  >  2:3416524/1‑148 (MQ=255)
              gCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGc                                                                                                                                 >  1:1994154/1‑148 (MQ=255)
                gTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTg                                                                                                                               >  1:1803680/1‑148 (MQ=255)
                gTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTg                                                                                                                               <  1:1627487/148‑1 (MQ=255)
                 tCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTgg                                                                                                                              <  1:2091999/148‑1 (MQ=255)
                  cTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGa                                                                                                                             <  1:733413/148‑1 (MQ=255)
                    tattaCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACg                                                                                                                           <  1:1855625/148‑1 (MQ=255)
                    tattaCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACg                                                                                                                           <  1:55294/148‑1 (MQ=255)
                       taCGCGTCGCGTGAAAAAGACACCCCCAGTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAg                                                                                                                        <  2:2657881/148‑1 (MQ=255)
                        aCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAgg                                                                                                                       >  1:588876/1‑148 (MQ=255)
                            gtcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTggg                                                                                                                   <  2:3308226/148‑1 (MQ=255)
                             tcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGc                                                                                                                  <  1:813699/148‑1 (MQ=255)
                             tcgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGACGCAGCTGGACGCAGGTGGGc                                                                                                                  >  1:1602581/1‑148 (MQ=255)
                              cgcgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCt                                                                                                                 >  2:2327440/1‑148 (MQ=255)
                                    aaaaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATAccc                                                                                                           >  2:2083301/1‑148 (MQ=255)
                                    aaaaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATAccc                                                                                                           >  2:3299897/1‑148 (MQ=255)
                                     aaaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACTGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACTCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCa                                                                                                          >  1:866678/1‑148 (MQ=255)
                                     aaaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCa                                                                                                          >  2:309135/1‑148 (MQ=255)
                                      aaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaa                                                                                                         >  1:2445536/1‑148 (MQ=255)
                                       aaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaaa                                                                                                        >  1:400052/1‑148 (MQ=255)
                                        aGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaaaa                                                                                                       <  1:67129/148‑1 (MQ=255)
                                        aGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaaaa                                                                                                       <  1:1596674/148‑1 (MQ=255)
                                         gACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAt                                                                                                      >  2:452465/1‑148 (MQ=255)
                                         gACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAt                                                                                                      >  1:3459625/1‑148 (MQ=255)
                                          acacCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                     >  2:1050257/1‑148 (MQ=255)
                                          acacCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                     <  2:2463578/148‑1 (MQ=255)
                                              ccccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGt                                                                                                 <  1:976911/148‑1 (MQ=255)
                                                ccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCa                                                                                               >  2:1682904/1‑148 (MQ=255)
                                                 cATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCAc                                                                                              >  1:3520622/1‑148 (MQ=255)
                                                 cATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCAc                                                                                              <  2:2624113/148‑1 (MQ=255)
                                                    tAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAt                                                                                           <  2:3451650/148‑1 (MQ=255)
                                                     aaaTATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAtt                                                                                          <  2:1189004/148‑1 (MQ=255)
                                                      aaTATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTg                                                                                         <  1:3082585/148‑1 (MQ=255)
                                                        tatGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTTATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGcc                                                                                       <  2:536422/148‑1 (MQ=255)
                                                        tatGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGcc                                                                                       >  1:2224199/1‑148 (MQ=255)
                                                         atGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCt                                                                                      >  1:2740147/1‑148 (MQ=255)
                                                          tGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTg                                                                                     >  2:340609/1‑148 (MQ=255)
                                                                cGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCa                                                                               <  1:299136/148‑1 (MQ=255)
                                                                  aaaCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACt                                                                             <  2:3521356/148‑1 (MQ=255)
                                                                   aaCTGTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAATAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTg                                                                            >  1:2634391/1‑148 (MQ=255)
                                                                       gTTCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggcgg                                                                        <  1:698674/148‑1 (MQ=255)
                                                                        ttCTGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggcggc                                                                       <  2:3044127/148‑1 (MQ=255)
                                                                           tGGCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGa                                                                    <  2:3201735/148‑1 (MQ=255)
                                                                            ggCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGaa                                                                   >  1:2538000/1‑148 (MQ=255)
                                                                             gCACACGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAg                                                                  >  1:627629/1‑148 (MQ=255)
                                                                                cacGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGc                                                               >  2:2373924/1‑148 (MQ=255)
                                                                                 acGGTCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCa                                                              <  1:1523454/148‑1 (MQ=255)
                                                                                     tCCGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAtt                                                          <  1:133041/148‑1 (MQ=255)
                                                                                      ccGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttt                                                         <  1:1196564/148‑1 (MQ=255)
                                                                                       cGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAtttt                                                        <  2:1271497/148‑1 (MQ=255)
                                                                                       cGGTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAtttt                                                        <  2:2241613/148‑1 (MQ=255)
                                                                                         gTACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtg                                                      >  1:938397/1‑148 (MQ=255)
                                                                                          tACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtgg                                                     <  2:134374/148‑1 (MQ=255)
                                                                                          tACGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtgg                                                     <  2:2281426/148‑1 (MQ=255)
                                                                                           aCGCAGCATTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTAGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtggt                                                    >  2:3260399/1‑148 (MQ=255)
                                                                                               agcaTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCt                                                >  1:2579525/1‑148 (MQ=255)
                                                                                                 caTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtgt                                              >  1:3523153/1‑148 (MQ=255)
                                                                                                   ttttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGa                                            <  2:630364/148‑1 (MQ=255)
                                                                                                   ttttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGa                                            <  1:110679/148‑1 (MQ=255)
                                                                                                          gcgcCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACACAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTAtt                                     >  2:3508042/1‑148 (MQ=255)
                                                                                                              ccTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGtt                                 >  1:1464549/1‑148 (MQ=255)
                                                                                                               cTGATCCGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTa                                <  2:1677466/148‑1 (MQ=255)
                                                                                                               cTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTa                                <  2:1192171/148‑1 (MQ=255)
                                                                                                                tGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTAt                               >  2:371928/1‑148 (MQ=255)
                                                                                                                 gATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATg                              <  2:1928634/148‑1 (MQ=255)
                                                                                                                  aTACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGa                             <  2:2359223/148‑1 (MQ=255)
                                                                                                                    aCGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACg                           <  2:3175127/148‑1 (MQ=255)
                                                                                                                     cGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGa                          >  2:1363629/1‑148 (MQ=255)
                                                                                                                        tGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCg                       >  1:734559/1‑148 (MQ=255)
                                                                                                                           aTTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGccc                    <  1:1447328/148‑1 (MQ=255)
                                                                                                                                aGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGCGCTATACCCAAACGATGGTCACTATTGCCTGAAGCCACTGGCGGCGGAAGTCCATATTTTTGGTTGCTGTGAAATTATTCGTTATGACGATCGCCCCTGGa               <  1:1733666/148‑1 (MQ=255)
                                                                                                                                 gTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAg              <  1:2327440/148‑1 (MQ=255)
                                                                                                                                  ttACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGc             <  1:1127147/148‑1 (MQ=255)
                                                                                                                                  ttACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGc             >  1:604560/1‑148 (MQ=255)
                                                                                                                                    aCGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGa           <  2:2445536/148‑1 (MQ=255)
                                                                                                                                     cGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaa          <  1:3416524/148‑1 (MQ=255)
                                                                                                                                       aagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaaa        >  2:2853160/1‑148 (MQ=255)
                                                                                                                                       aagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaaaa        >  1:2401523/1‑148 (MQ=255)
                                                                                                                                         gaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAag      <  2:866678/148‑1 (MQ=255)
                                                                                                                                         gaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAag      <  2:1493532/148‑1 (MQ=255)
                                                                                                                                         gaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAag      >  2:1951854/1‑148 (MQ=255)
                                                                                                                                          aaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAaga     >  1:3009467/1‑148 (MQ=255)
                                                                                                                                          aaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAaga     >  1:2164242/1‑148 (MQ=255)
                                                                                                                                          aaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAaga     <  2:1639853/148‑1 (MQ=255)
                                                                                                                                           aTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAagag    <  1:1980478/148‑1 (MQ=255)
                                                                                                                                           aTTTTGAGAAACCTAGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAagag    <  1:3260399/148‑1 (MQ=255)
                                                                                                                                            ttttGAGAAACCTGGCCGGAGGTGGACGCAGGTGGGATATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGTCGGCGGAAGGGCTTATTTTTGGTTGCTGTTAATTTATTCGTTATGACGATGGCCGCTGGAGCGATAAAGAGt   <  2:2634391/148‑1 (MQ=255)
                                                                                                                                             tttGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGtt  >  2:2855613/1‑148 (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: ‑