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,0020TA99.7% 1014.0 / ‑5.1 305C102* (TGT→TGA) fhlAformate hydrogenlyase transcriptional activator
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (153/151);  total (153/152)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 306 total aligned reads displayed.

AGCGCAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTTC  >  CP009273/2847855‑2848146
                                                                                                                                                    |                                                                                                                                                
   gcAGACTCAACGTGCGTCTTATTCCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAc                                                                                                                                             <  1:4674660/150‑1 (MQ=255)
   gcAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAc                                                                                                                                             >  2:2950245/1‑150 (MQ=255)
    cAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                            <  1:4511759/150‑1 (MQ=255)
    cAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAAcc                                                                                                                                            >  2:3722585/1‑150 (MQ=255)
      gACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTg                                                                                                                                          <  2:630855/150‑1 (MQ=255)
       aCTCAACGTGCTTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTgg                                                                                                                                         <  1:3798039/150‑1 (MQ=255)
       aCTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTgg                                                                                                                                         <  1:4423411/150‑1 (MQ=255)
        cTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGACACCTGGc                                                                                                                                        >  1:998038/1‑150 (MQ=255)
        cTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGc                                                                                                                                        >  1:2796235/1‑150 (MQ=255)
        cTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGc                                                                                                                                        >  1:4950041/1‑150 (MQ=255)
         tCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGcc                                                                                                                                       >  1:4269105/1‑150 (MQ=255)
           aaCGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGc                                                                                                                                     >  2:55404/1‑150 (MQ=255)
               tGCGTCTTTTTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGCAGACACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCATTTACGAAGAATTTTGAGAAACCTGGCCTCAGCt                                                                                                                                 <  1:1666405/150‑1 (MQ=255)
               tGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCt                                                                                                                                 <  2:840426/150‑1 (MQ=255)
                  gTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGa                                                                                                                              >  1:4750138/1‑150 (MQ=255)
                  gTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGa                                                                                                                              >  2:1386420/1‑150 (MQ=255)
                    cTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACg                                                                                                                            >  2:3280936/1‑150 (MQ=255)
                      tattaCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCa                                                                                                                          <  1:4626412/150‑1 (MQ=255)
                       attaCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAg                                                                                                                         >  2:397806/1‑150 (MQ=255)
                        ttaCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGACACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAgg                                                                                                                        <  2:2838088/150‑1 (MQ=255)
                        ttaCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAgg                                                                                                                        >  1:344336/1‑150 (MQ=255)
                         taCGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGt                                                                                                                       >  2:74274/1‑150 (MQ=255)
                           cgcgtcgagtGAAAAAGACAACCCCATTAAATATGAAGACGAAACTTTTCTGTCACACGGTCCGGTACGCAGCA‑TTTTGTCTCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTgg                                                                                                                     <  1:79989/150‑1 (MQ=255)
                                  cgtGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATAc                                                                                                              <  2:4047364/150‑1 (MQ=255)
                                    tGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATAccc                                                                                                            >  2:1701691/1‑150 (MQ=255)
                                     gAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCa                                                                                                           >  2:1069508/1‑150 (MQ=255)
                                       aaaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaaa                                                                                                         >  2:5045983/1‑150 (MQ=255)
                                        aaaGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCaaaa                                                                                                        <  1:3102799/150‑1 (MQ=255)
                                          aGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                      >  2:1466656/1‑150 (MQ=255)
                                          aGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                      <  2:4968344/150‑1 (MQ=255)
                                          aGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAAtt                                                                                                      >  1:2005265/1‑150 (MQ=255)
                                                ccccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCa                                                                                                <  1:704247/150‑1 (MQ=255)
                                                ccccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCa                                                                                                <  2:5293376/150‑1 (MQ=255)
                                                 cccATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCAc                                                                                               <  2:1485516/150‑1 (MQ=255)
                                                   cATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTa                                                                                             >  2:2721899/1‑150 (MQ=255)
                                                    aTTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAt                                                                                            >  2:5082022/1‑150 (MQ=255)
                                                    aTTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAt                                                                                            <  2:5213167/150‑1 (MQ=255)
                                                     ttAAATATGCAGACGACACTGTACTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAtt                                                                                           <  1:2731182/150‑1 (MQ=255)
                                                     ttAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTAtt                                                                                           >  2:2067286/1‑150 (MQ=255)
                                                      tAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTg                                                                                          <  1:1085451/150‑1 (MQ=255)
                                                       aaaTATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGc                                                                                         >  1:1745267/1‑150 (MQ=255)
                                                        aaTATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGcc                                                                                        >  1:4563539/1‑150 (MQ=255)
                                                        aaTATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGcc                                                                                        >  2:729336/1‑150 (MQ=255)
                                                         atatGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCt                                                                                       >  2:1856910/1‑150 (MQ=255)
                                                          tatGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTg                                                                                      >  2:4135165/1‑150 (MQ=255)
                                                          tatGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTg                                                                                      >  1:1258313/1‑150 (MQ=255)
                                                           atGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGa                                                                                     >  2:2442726/1‑150 (MQ=255)
                                                            tGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGAt                                                                                    >  2:2888454/1‑150 (MQ=255)
                                                            tGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGAt                                                                                    <  2:935023/150‑1 (MQ=255)
                                                             gAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATg                                                                                   <  1:1832437/150‑1 (MQ=255)
                                                              aaGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCATTTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATg                                                                                   <  1:3198187/150‑1 (MQ=255)
                                                              aaGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGc                                                                                  >  1:3407080/1‑150 (MQ=255)
                                                              aaGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGc                                                                                  >  2:2402329/1‑150 (MQ=255)
                                                               aGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGTCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGcc                                                                                 >  2:4859890/1‑150 (MQ=255)
                                                                gACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCa                                                                                <  2:4112702/150‑1 (MQ=255)
                                                                  cGAAACTGTTCTGGCACACGGTCCGGTACGCTGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACt                                                                              <  2:437390/150‑1 (MQ=255)
                                                                   gAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTg                                                                             >  1:1210778/1‑150 (MQ=255)
                                                                    aaaCTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTgg                                                                            >  2:233152/1‑150 (MQ=255)
                                                                     aaCTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggc                                                                           >  2:3089835/1‑150 (MQ=255)
                                                                     aaCTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggc                                                                           <  1:2141058/150‑1 (MQ=255)
                                                                     aaCTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTggc                                                                           <  1:3042883/150‑1 (MQ=255)
                                                                           tCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGa                                                                     <  1:4078829/150‑1 (MQ=255)
                                                                                cacaCGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGc                                                                <  1:2165400/150‑1 (MQ=255)
                                                                                 acacGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCa                                                               >  1:643247/1‑150 (MQ=255)
                                                                                 acacGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCa                                                               >  1:5183807/1‑150 (MQ=255)
                                                                                  cacGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCat                                                              <  1:911274/150‑1 (MQ=255)
                                                                                        ccGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttttt                                                        >  2:2142573/1‑150 (MQ=255)
                                                                                        ccGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttttt                                                        <  1:3505839/150‑1 (MQ=255)
                                                                                        acGGTACGCAGCA‑TTGTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTTGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATAttttt                                                        <  1:1842704/149‑1 (MQ=255)
                                                                                          ggTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtgg                                                      >  1:1071064/1‑150 (MQ=255)
                                                                                           gTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtggt                                                     <  1:4159118/150‑1 (MQ=255)
                                                                                            tACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtggtg                                                    >  2:4797143/1‑150 (MQ=255)
                                                                                             aCGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTtggtgg                                                   <  1:4377906/150‑1 (MQ=255)
                                                                                                cagca‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCATAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtg                                                <  2:1315493/150‑1 (MQ=255)
                                                                                                cagca‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtg                                                <  2:4640635/150‑1 (MQ=255)
                                                                                                  gca‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCtgtg                                              <  1:4984997/150‑1 (MQ=255)
                                                                                                      ttttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAAt                                           <  1:3823307/150‑1 (MQ=255)
                                                                                                       tttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAAtt                                          >  1:5079669/1‑150 (MQ=255)
                                                                                                       tttGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAAtt                                          >  1:2435336/1‑150 (MQ=255)
                                                                                                          gTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTAt                                       >  2:4697480/1‑150 (MQ=255)
                                                                                                          gTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTAt                                       >  1:4582997/1‑150 (MQ=255)
                                                                                                           tCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTAtt                                      <  1:3348700/150‑1 (MQ=255)
                                                                                                              cgcCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGt                                   >  1:357710/1‑150 (MQ=255)
                                                                                                              cgcCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGt                                   <  1:2048843/150‑1 (MQ=255)
                                                                                                                 ccTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTAt                                <  1:2884320/150‑1 (MQ=255)
                                                                                                                 ccTGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTAt                                <  1:1346148/150‑1 (MQ=255)
                                                                                                                   tGATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGTTGGCTGTGAATTTATTCGTTATGa                              >  2:3349574/1‑150 (MQ=255)
                                                                                                                    gATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGAc                             <  1:3280936/150‑1 (MQ=255)
                                                                                                                    gATACGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGAc                             <  1:4905822/150‑1 (MQ=255)
                                                                                                                        cGCTGCATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATc                         <  2:49086/150‑1 (MQ=255)
                                                                                                                             cATTGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGcccc                    <  1:1726050/150‑1 (MQ=255)
                                                                                                                                tGCAGTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTgg                 >  2:115922/1‑150 (MQ=255)
                                                                                                                                    gTTACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCg             >  2:2568490/1‑150 (MQ=255)
                                                                                                                                      tACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaa           <  2:4502526/150‑1 (MQ=255)
                                                                                                                                      tACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaa           >  1:2293685/1‑150 (MQ=255)
                                                                                                                                      tACGAAGAATTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGaa           >  2:3494897/1‑150 (MQ=255)
                                                                                                                                          aagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAag       >  2:1641659/1‑150 (MQ=255)
                                                                                                                                          aagaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAag       <  2:2712086/150‑1 (MQ=255)
                                                                                                                                           agaaTTTTGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAaga      >  2:2271136/1‑150 (MQ=255)
                                                                                                                                               ttttGAGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTTc  <  1:2624196/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
AGCGCAGACTCAACGTGCGTCTTATTACGCGTCGCGTGAAAAAGACACCCCCATTAAATATGAAGACGAAACTGTTCTGGCACACGGTCCGGTACGCAGCA‑TTTTGTCGCGCCCTGATACGCTGCATTGCAGTTACGAAGAATTTTGTGAAACCTGGCCGCAGCTGGACGCAGGTGGGCTATACCCAAAATTTGGTCACTATTGCCTGATGCCACTGGCGGCGGAAGGGCATATTTTTGGTGGCTGTGAATTTATTCGTTATGACGATCGCCCCTGGAGCGAAAAAGAGTTC  >  CP009273/2847855‑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: ‑