Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,586,871 G→T 86.2% intergenic (‑172/+51) hipB ← / ← yneO antitoxin of HipAB toxin‑antitoxin system/pseudogene, AidA homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,586,8710GT86.2% 542.2 / 90.9 236intergenic (‑172/+51)hipB/yneOantitoxin of HipAB toxin‑antitoxin system/pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base G (11/22);  new base T (107/96);  total (118/118)
Fisher's exact test for biased strand distribution p-value = 5.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 238 total aligned reads displayed.

CTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGT  >  CP009273/1586727‑1587016
                                                                                                                                                |                                                                                                                                                 
    aGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCgg                                                                                                                                            >  1:963714/1‑148 (MQ=255)
    aGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCAGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCgg                                                                                                                                            >  2:55139/1‑148 (MQ=255)
     ggggATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGc                                                                                                                                           >  1:73007/1‑148 (MQ=255)
      gggATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCa                                                                                                                                          <  2:879488/148‑1 (MQ=255)
         atatTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAata                                                                                                                                       <  2:2110599/148‑1 (MQ=255)
         atatTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAata                                                                                                                                       >  2:2694183/1‑148 (MQ=255)
           attatAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATAtt                                                                                                                                     >  1:1082734/1‑148 (MQ=255)
            ttatAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTa                                                                                                                                    <  1:1470838/148‑1 (MQ=255)
            ttatAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTa                                                                                                                                    >  2:413962/1‑148 (MQ=255)
            ttatAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTa                                                                                                                                    >  1:3103935/1‑148 (MQ=255)
             tataAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTaa                                                                                                                                   >  1:535020/1‑148 (MQ=255)
               taAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGc                                                                                                                                 <  2:2290573/148‑1 (MQ=255)
               taAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGc                                                                                                                                 <  2:384817/148‑1 (MQ=255)
               taAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGc                                                                                                                                 >  1:2197416/1‑148 (MQ=255)
                  gTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGAGAAAAAGGAGGATTTAAATAAGGCGATGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGc                                                                                                                              >  2:3155882/1‑148 (MQ=255)
                   ttttATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCa                                                                                                                             >  1:9324/1‑148 (MQ=255)
                   ttttATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCa                                                                                                                             >  1:3137474/1‑148 (MQ=255)
                     ttATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATg                                                                                                                           <  1:3304927/148‑1 (MQ=255)
                      tATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGc                                                                                                                          >  1:2218344/1‑148 (MQ=255)
                       aTCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGcc                                                                                                                         >  1:2159778/1‑148 (MQ=255)
                         ccTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTAGCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCAt                                                                                                                       <  1:3119828/148‑1 (MQ=255)
                          cTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCAtc                                                                                                                      <  2:2863296/148‑1 (MQ=255)
                          cTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCAtc                                                                                                                      <  2:963714/148‑1 (MQ=255)
                           tttAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCAtct                                                                                                                     >  1:238695/1‑148 (MQ=255)
                               gTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGAc                                                                                                                 >  1:1191224/1‑148 (MQ=255)
                                tGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACa                                                                                                                >  2:3032168/1‑148 (MQ=255)
                                tGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACa                                                                                                                >  2:113383/1‑148 (MQ=255)
                                tGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACa                                                                                                                >  1:3084162/1‑148 (MQ=255)
                                      tAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCTGATAAAAAGGAGGATTTAACTAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGc                                                                                                          <  2:170459/148‑1 (MQ=255)
                                       aaGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCt                                                                                                         >  1:2448037/1‑148 (MQ=255)
                                         gTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTta                                                                                                       <  1:616977/148‑1 (MQ=255)
                                         gTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTta                                                                                                       <  2:2605181/148‑1 (MQ=255)
                                         gTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTta                                                                                                       <  2:430859/148‑1 (MQ=255)
                                           cAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTtata                                                                                                     >  2:1261412/1‑148 (MQ=255)
                                               cTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGt                                                                                                 >  2:2062552/1‑148 (MQ=255)
                                                 ggAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGttt                                                                                               <  1:2959433/148‑1 (MQ=255)
                                                  gAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTa                                                                                              >  2:3038669/1‑148 (MQ=255)
                                                    aaCACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTAcc                                                                                            <  1:497986/148‑1 (MQ=255)
                                                     acacATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCa                                                                                           >  2:486806/1‑148 (MQ=255)
                                                     acacATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCa                                                                                           >  2:1593884/1‑148 (MQ=255)
                                                     acacATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCa                                                                                           <  1:503307/148‑1 (MQ=255)
                                                      cacaTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAg                                                                                          <  1:1654918/148‑1 (MQ=255)
                                                      cacaTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAg                                                                                          <  1:2789699/148‑1 (MQ=255)
                                                       acaTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGc                                                                                         <  1:525226/148‑1 (MQ=255)
                                                        caTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCt                                                                                        >  2:3168625/1‑148 (MQ=255)
                                                        caTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCt                                                                                        >  1:1442910/1‑148 (MQ=255)
                                                        caTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCt                                                                                        >  2:3204566/1‑148 (MQ=255)
                                                        caTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCt                                                                                        <  2:1817412/148‑1 (MQ=255)
                                                                    tctACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTcc                                                                            >  1:950555/1‑148 (MQ=255)
                                                                    tctACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTcc                                                                            <  1:176880/148‑1 (MQ=255)
                                                                     ctACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTccc                                                                           >  1:2379187/1‑148 (MQ=255)
                                                                       aCGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAg                                                                         >  2:1597084/1‑148 (MQ=255)
                                                                         gggATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCa                                                                       <  2:3313365/148‑1 (MQ=255)
                                                                         gggATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCa                                                                       <  1:2014002/148‑1 (MQ=255)
                                                                            aTAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCCGCATGc                                                                    >  2:2561817/1‑148 (MQ=255)
                                                                             tAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGcc                                                                   <  1:496631/148‑1 (MQ=255)
                                                                              aaaCGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGccc                                                                  >  1:2188644/1‑148 (MQ=255)
                                                                              aaaCGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGccc                                                                  <  2:1754357/148‑1 (MQ=255)
                                                                                aCGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTg                                                                <  2:131495/148‑1 (MQ=255)
                                                                                aCGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTAACTCCCAGCATGCCCTg                                                                >  1:84985/1‑148 (MQ=255)
                                                                                   cTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGt                                                             >  1:657034/1‑148 (MQ=255)
                                                                                   cTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGt                                                             <  2:736203/148‑1 (MQ=255)
                                                                                   cTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGAGTTGATGCGCCGTTTTAATTCGGCAAAAATAAGCCGCATGCCCTATACACATGAGGCGTATAACGTTTACCACCAAAATTTCAATCCCAGCAAGCCCAGAGt                                                             >  1:1498505/1‑148 (MQ=255)
                                                                                         tatCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCt                                                       <  1:485434/148‑1 (MQ=255)
                                                                                          atCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGCCATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCta                                                      >  1:1122155/1‑148 (MQ=255)
                                                                                            ccGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCtata                                                    >  2:226464/1‑148 (MQ=255)
                                                                                            ccGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCtata                                                    >  2:465149/1‑148 (MQ=255)
                                                                                            ccGCGATCGCGGAAAAAAAGGAGGATTTAAATAAGGAGTTGATGCGCCGTTTTAAGTAGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATAATTCACAACCAGCATGACCTGAGTAACGcaaaa                                                    >  2:1225028/1‑144 (MQ=255)
                                                                                               cgATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGcc                                                 <  2:203294/148‑1 (MQ=255)
                                                                                                gATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCt                                                <  2:2402897/148‑1 (MQ=255)
                                                                                                gATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCt                                                <  1:1563853/148‑1 (MQ=255)
                                                                                                 aTCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCtt                                               <  1:509440/148‑1 (MQ=255)
                                                                                                  tCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCttt                                              <  1:1357563/148‑1 (MQ=255)
                                                                                                  tCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCttt                                              <  2:2197416/148‑1 (MQ=255)
                                                                                                   cgcgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTa                                             <  1:934217/148‑1 (MQ=255)
                                                                                                   cgcgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTa                                             <  1:1946369/148‑1 (MQ=255)
                                                                                                     cgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATc                                           <  2:2998819/148‑1 (MQ=255)
                                                                                                          aaaaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTa                                      >  2:2799949/1‑148 (MQ=255)
                                                                                                           aaaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAg                                     >  2:3248999/1‑148 (MQ=255)
                                                                                                            aaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGt                                    <  1:2060926/148‑1 (MQ=255)
                                                                                                            aaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGt                                    <  1:396032/148‑1 (MQ=255)
                                                                                                             aaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGtt                                   <  2:1900100/148‑1 (MQ=255)
                                                                                                                gaggaTTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGca                                <  2:3103935/148‑1 (MQ=255)
                                                                                                                     tttAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCg                           >  1:74410/1‑148 (MQ=255)
                                                                                                                     tttAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCg                           <  1:3168625/148‑1 (MQ=255)
                                                                                                                     tttAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTAACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATAACCTAGTTGCACACCg                           >  1:2498037/1‑148 (MQ=255)
                                                                                                                      ttAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGa                          >  1:1066065/1‑148 (MQ=255)
                                                                                                                       tAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGAc                         <  2:1094662/148‑1 (MQ=255)
                                                                                                                         aaTAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACAt                       <  2:73007/148‑1 (MQ=255)
                                                                                                                           tAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTc                     >  2:2625183/1‑148 (MQ=255)
                                                                                                                             aGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTccc                   <  1:1624967/148‑1 (MQ=255)
                                                                                                                              ggCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTcccc                  >  2:2146496/1‑148 (MQ=255)
                                                                                                                              ggCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTcccc                  <  1:2702665/148‑1 (MQ=255)
                                                                                                                               gCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTccccc                 >  2:1487322/1‑148 (MQ=255)
                                                                                                                                cGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCca                >  2:1038891/1‑148 (MQ=255)
                                                                                                                                     aTGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGc           >  1:35460/1‑148 (MQ=255)
                                                                                                                                        cgcCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTaa        <  2:2188644/148‑1 (MQ=255)
                                                                                                                                        cgcCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTaa        <  1:1582323/148‑1 (MQ=255)
                                                                                                                                             ttttAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGgttg   <  1:792139/148‑1 (MQ=255)
                                                                                                                                              ttGTAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGTCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGgttgt  <  2:2209954/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGT  >  CP009273/1586727‑1587016

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: