Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,629,424 T→A 100% intergenic (+171/+144) yhiM → / ← yhiN acid resistance protein, inner membrane/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,629,4240TA94.5% 543.2 / 18.4 184intergenic (+171/+144)yhiM/yhiNacid resistance protein, inner membrane/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (9/1);  new base A (76/95);  total (88/96)
Fisher's exact test for biased strand distribution p-value = 6.68e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 185 total aligned reads displayed.

GTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAATAGGT  >  CP009273/3629279‑3629568
                                                                                                                                                 |                                                                                                                                                
gTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGtatt                                                                                                                                                <  2:1332316/148‑1 (MQ=255)
  tGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCa                                                                                                                                              >  2:2327050/1‑148 (MQ=255)
  tGACCATTTGAGTAACGCGGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCACCCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCa                                                                                                                                              <  2:1706171/148‑1 (MQ=255)
    aCCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAAt                                                                                                                                            <  1:1177855/148‑1 (MQ=255)
      cATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCa                                                                                                                                          >  1:2313791/1‑148 (MQ=255)
         ttGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCt                                                                                                                                       <  1:579775/148‑1 (MQ=255)
         ttGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCt                                                                                                                                       <  2:1717887/148‑1 (MQ=255)
         ttGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCt                                                                                                                                       <  2:508870/148‑1 (MQ=255)
             gTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTtcc                                                                                                                                   <  1:1030645/148‑1 (MQ=255)
               aaCGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTtccat                                                                                                                                 >  2:3256280/1‑148 (MQ=255)
                 cGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAGAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTtccatcc                                                                                                                               >  2:681488/1‑148 (MQ=255)
                 cGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTtccatcc                                                                                                                               >  1:2393582/1‑148 (MQ=255)
                 cGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTtccatcc                                                                                                                               >  2:683919/1‑148 (MQ=255)
                    ttgTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCAcc                                                                                                                            <  2:3097363/146‑1 (MQ=255)
                    tGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCAcc                                                                                                                            >  2:1620661/1‑148 (MQ=255)
                    tGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCAcc                                                                                                                            <  1:2423802/148‑1 (MQ=255)
                     ggTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCt                                                                                                                           <  1:648973/148‑1 (MQ=255)
                       tttttCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCtga                                                                                                                         >  2:3373463/1‑146 (MQ=255)
                        ttttCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTtat                                                                                                                        <  2:766198/148‑1 (MQ=255)
                        ttttCTTAGGAATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTtat                                                                                                                        <  2:2334715/148‑1 (MQ=255)
                          ttCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTatat                                                                                                                      <  1:635395/148‑1 (MQ=255)
                            cTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTa                                                                                                                    <  1:2988099/148‑1 (MQ=255)
                              tAGGCATCATGAAATAACGCATCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAg                                                                                                                  <  2:2765360/148‑1 (MQ=255)
                              tAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAg                                                                                                                  >  2:2611654/1‑148 (MQ=255)
                              tAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAg                                                                                                                  <  2:3074574/148‑1 (MQ=255)
                               aGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGc                                                                                                                 <  1:2122938/148‑1 (MQ=255)
                               aGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGc                                                                                                                 <  1:3107154/148‑1 (MQ=255)
                                ggCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTTATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCa                                                                                                                >  1:234168/1‑148 (MQ=255)
                                ggCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCa                                                                                                                >  1:1258917/1‑148 (MQ=255)
                                   atcatGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATgg                                                                                                             <  1:2183953/148‑1 (MQ=255)
                                     catGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAg                                                                                                           >  2:2191639/1‑148 (MQ=255)
                                     catGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAg                                                                                                           >  2:1802662/1‑148 (MQ=255)
                                       tGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAggg                                                                                                         <  1:263689/148‑1 (MQ=255)
                                         aaaTAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGtt                                                                                                       >  1:1136067/1‑148 (MQ=255)
                                              aCGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGt                                                                                                  <  1:1676248/148‑1 (MQ=255)
                                              aCGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGAAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTAAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGt                                                                                                  <  1:3440146/148‑1 (MQ=255)
                                                    caTTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACgg                                                                                            >  1:2725023/1‑148 (MQ=255)
                                                    caTTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACgg                                                                                            >  1:2477444/1‑148 (MQ=255)
                                                    caTTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACgg                                                                                            >  1:209258/1‑148 (MQ=255)
                                                    caTTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGATACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAAACACCTTTCCATCCACCTAATATTAAGCATGGAGGGTTTCAGTTCACgg                                                                                            >  1:3381028/1‑148 (MQ=255)
                                                        aaTGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCCCCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTc                                                                                        <  2:996907/148‑1 (MQ=255)
                                                         aTGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCa                                                                                       >  1:1049214/1‑148 (MQ=255)
                                                                gTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGaaa                                                                                <  2:2056649/148‑1 (MQ=255)
                                                                      aaGTATAAAAAACAGCATAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATc                                                                          <  1:757636/148‑1 (MQ=255)
                                                                      aaGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTCTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATc                                                                          >  1:97982/1‑148 (MQ=255)
                                                                        gTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATcac                                                                        <  2:211886/148‑1 (MQ=255)
                                                                         tataAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATcaca                                                                       >  1:3120076/1‑148 (MQ=255)
                                                                           taAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCTGTTCACGGGCTCATTAGAAAATAATCACAAg                                                                     <  2:3236284/148‑1 (MQ=255)
                                                                           taAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAg                                                                     >  2:3137142/1‑148 (MQ=255)
                                                                            aaaaaaCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGa                                                                    <  2:3052810/148‑1 (MQ=255)
                                                                            aaaaaaCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGa                                                                    >  2:295958/1‑148 (MQ=255)
                                                                            aaaaaaCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGa                                                                    <  1:359822/148‑1 (MQ=255)
                                                                               aaaCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATaa                                                                 >  1:1010158/1‑148 (MQ=255)
                                                                               aaaCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCACCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATaa                                                                 <  2:2057561/148‑1 (MQ=255)
                                                                                aaCAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAAc                                                                >  1:2376383/1‑148 (MQ=255)
                                                                                  cAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCa                                                              >  2:3317418/1‑148 (MQ=255)
                                                                                  cAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCa                                                              >  2:901664/1‑148 (MQ=255)
                                                                                   aGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCAc                                                             >  2:1566015/1‑148 (MQ=255)
                                                                                     cAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTa                                                           <  1:730405/148‑1 (MQ=255)
                                                                                     cAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTa                                                           >  2:648347/1‑148 (MQ=255)
                                                                                       aaGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATc                                                         >  1:2660225/1‑148 (MQ=255)
                                                                                       aaGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATc                                                         >  1:2660140/1‑148 (MQ=255)
                                                                                          ttACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCaaa                                                      <  2:695413/148‑1 (MQ=255)
                                                                                          ttACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCaaa                                                      >  2:1995910/1‑148 (MQ=255)
                                                                                          ttACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCaaa                                                      <  1:1802662/148‑1 (MQ=255)
                                                                                             cTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAAtt                                                   >  1:2674517/1‑148 (MQ=255)
                                                                                              tGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTa                                                  >  1:3480144/1‑148 (MQ=255)
                                                                                              tGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTa                                                  >  2:2275904/1‑148 (MQ=255)
                                                                                              tGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTa                                                  >  2:1533160/1‑148 (MQ=255)
                                                                                                 ttttttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCt                                               >  1:38746/1‑148 (MQ=255)
                                                                                                 ttttttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCt                                               <  1:3137142/148‑1 (MQ=255)
                                                                                                 ttttttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCt                                               <  2:2361502/148‑1 (MQ=255)
                                                                                                 ttttttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCt                                               <  2:1049214/148‑1 (MQ=255)
                                                                                                    tttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCAt                                            <  1:1533160/148‑1 (MQ=255)
                                                                                                    tttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCAt                                            <  2:1136159/148‑1 (MQ=255)
                                                                                                     ttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCAtt                                           <  1:2522191/148‑1 (MQ=255)
                                                                                                     ttCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCAtt                                           <  2:1021112/148‑1 (MQ=255)
                                                                                                      tCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCAttt                                          <  1:3533457/148‑1 (MQ=255)
                                                                                                         aCCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTa                                       <  2:1692939/148‑1 (MQ=255)
                                                                                                            tGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATg                                    >  1:2797072/1‑148 (MQ=255)
                                                                                                                   atTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGActct                             >  1:2069888/1‑148 (MQ=255)
                                                                                                                   atTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGActct                             >  1:3557472/1‑148 (MQ=255)
                                                                                                                           aaGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGATCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAAttt                     >  2:3202385/1‑148 (MQ=255)
                                                                                                                            aGATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGTGGGTTTAAGTTCACGGGCTCATTAGAAAATAATCACAATATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAAtttt                    <  1:2281368/148‑1 (MQ=255)
                                                                                                                             gATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTc                   <  1:1447484/148‑1 (MQ=255)
                                                                                                                             gATCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTc                   >  2:189511/1‑148 (MQ=255)
                                                                                                                               tCTGGCCACTGAGTGTTTATTCAATCACCTTTCCATCCACCTTATAATAAGCATGGAGGGTTTCATTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCAt                 <  1:3001213/148‑1 (MQ=255)
                                                                                                                               tCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCAt                 >  1:1080451/1‑148 (MQ=255)
                                                                                                                               tCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCAt                 <  1:593892/148‑1 (MQ=255)
                                                                                                                               tCTGGCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCAt                 <  1:593381/148‑1 (MQ=255)
                                                                                                                                 tGGCCACTGAGTGATTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCt               <  2:2154760/148‑1 (MQ=255)
                                                                                                                                  ggCCACTGAGTGAGTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCtt              >  1:194811/1‑148 (MQ=255)
                                                                                                                                        tgagtgagTATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGaa        <  2:3495103/148‑1 (MQ=255)
                                                                                                                                         gagtgagtATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGaaa       >  2:1210564/1‑148 (MQ=255)
                                                                                                                                          agtgagtATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAAt      >  2:1782753/1‑148 (MQ=255)
                                                                                                                                             gagtATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAATAgg   <  2:3443159/148‑1 (MQ=255)
                                                                                                                                             gagtATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAATAgg   <  1:1995910/148‑1 (MQ=255)
                                                                                                                                             gagtATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAATAgg   <  2:2725023/148‑1 (MQ=255)
                                                                                                                                             gagtATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAATAgg   <  2:2477444/148‑1 (MQ=255)
                                                                                                                                              agtATTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAATAGGt  >  2:3552650/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GTTGACCATTTGCGTAACGCTGGTTTTTCTTAGGCATCATGAAATAACGCAGCATTAATGCATAGTGGTTAAGTATAAAAAACAGCAAAGTTACTGTTTTTTTCAACCTGTTCATATTTCATAAAGATCTGGCCACTGAGTGAGTTTTCAATCACCTTTCCATCCACCTTATATTAAGCATGGAGGGTTTCAGTTCACGGGCTCATTAGAAAATAATCACAAGATAACCACTATCAAAATTAGCTCATTTTTAATGCGACTCTAATAATTTTCATCTTTAGGAAATAGGT  >  CP009273/3629279‑3629568

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: