Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,586,871 G→T 26.7% 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,8710GT26.7% 705.6 / 212.0 306intergenic (‑172/+51)hipB/yneOantitoxin of HipAB toxin‑antitoxin system/pseudogene, AidA homolog
Reads supporting (aligned to +/- strand):  ref base G (118/106);  new base T (48/34);  total (166/140)
Fisher's exact test for biased strand distribution p-value = 3.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 310 total aligned reads displayed.

ATCCGCTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATT  >  CP009273/1586722‑1587020
                                                                                                                                                     |                                                                                                                                                     
aTCCGCTTAAGGGGATATTATAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGa                                                                                                                                                      <  1:3258619/150‑1 (MQ=255)
  ccGCTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTaa                                                                                                                                                     >  1:916301/1‑150 (MQ=255)
     cTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTc                                                                                                                                                  <  1:445392/150‑1 (MQ=255)
     cTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTc                                                                                                                                                  >  1:3350786/1‑150 (MQ=255)
     cTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTc                                                                                                                                                  <  1:1579963/150‑1 (MQ=255)
      ttAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCg                                                                                                                                                 <  1:3435755/150‑1 (MQ=255)
      ttAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCg                                                                                                                                                 <  1:2199417/150‑1 (MQ=255)
        aaGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGc                                                                                                                                               >  2:2853774/1‑150 (MQ=255)
        aaGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGc                                                                                                                                               >  1:1320367/1‑150 (MQ=255)
         aGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCAGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAATGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAATTCGGCa                                                                                                                                              >  2:3269547/1‑150 (MQ=255)
           gggATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAat                                                                                                                                            <  2:1085836/150‑1 (MQ=255)
           gggATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAat                                                                                                                                            <  1:3944254/150‑1 (MQ=255)
            ggATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAata                                                                                                                                           <  1:3094165/150‑1 (MQ=255)
             gATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAatat                                                                                                                                          <  1:2126810/150‑1 (MQ=255)
                attatAAGTTTTATCCTTTAGTGAGGATAAGTCACTCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCAGCGATATCCGCTATCGCGGATAAAAAGTAGGATTTAAATAAGGCGTTGATGCGCAGTTTTAAGTAGGAAATATTaa                                                                                                                                       >  2:1181187/1‑150 (MQ=255)
                       gTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCAt                                                                                                                                <  2:4544576/150‑1 (MQ=255)
                         tttATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGc                                                                                                                              >  2:1393667/1‑150 (MQ=255)
                          ttATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGcc                                                                                                                             >  1:3614092/1‑150 (MQ=255)
                           tATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGGCGGCAATATTAAGCCGCATGCCa                                                                                                                            >  2:1253160/1‑150 (MQ=255)
                            aTCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCAt                                                                                                                           >  2:489763/1‑150 (MQ=255)
                            aTCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCATATACGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCAt                                                                                                                           >  2:1389956/1‑150 (MQ=255)
                             tCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCTTTTGAAGTCGGCAATATTAAGCCGCATGCCAtc                                                                                                                          >  2:1675961/1‑150 (MQ=255)
                              ccTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCAtct                                                                                                                         >  2:4100063/1‑150 (MQ=255)
                               cTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCAtctc                                                                                                                        <  2:1836121/150‑1 (MQ=255)
                                tttAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCg                                                                                                                       <  1:691448/150‑1 (MQ=255)
                                  tAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGAc                                                                                                                     <  2:3684884/150‑1 (MQ=255)
                                   aGTGAGGATACGTCAATCTGGAAACACATCTTATCCGCTCTCCGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCAGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACa                                                                                                                    >  2:3478115/1‑150 (MQ=255)
                                      gAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGc                                                                                                                 <  1:2538080/150‑1 (MQ=255)
                                          aTAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCt                                                                                                             <  2:2902315/150‑1 (MQ=255)
                                          aTAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCt                                                                                                             <  1:1155846/150‑1 (MQ=255)
                                           tAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCtt                                                                                                            <  2:3438048/150‑1 (MQ=255)
                                           tAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCtt                                                                                                            <  1:2431797/150‑1 (MQ=255)
                                             aGTCTATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTtat                                                                                                          <  2:3962634/150‑1 (MQ=255)
                                             aGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTtat                                                                                                          >  1:568749/1‑150 (MQ=255)
                                               tCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATAc                                                                                                        >  2:4476570/1‑150 (MQ=255)
                                               tCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATAc                                                                                                        <  1:1773037/150‑1 (MQ=255)
                                                  aTCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGt                                                                                                     <  2:3614244/150‑1 (MQ=255)
                                                   tCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGtt                                                                                                    <  2:59336/150‑1 (MQ=255)
                                                     tGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTa                                                                                                  <  1:2130557/150‑1 (MQ=255)
                                                     tGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTa                                                                                                  <  2:4377281/150‑1 (MQ=255)
                                                      ggAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTAc                                                                                                 <  2:3086898/150‑1 (MQ=255)
                                                        aaaCACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCa                                                                                               >  2:1610293/1‑150 (MQ=255)
                                                          acacATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGc                                                                                             <  2:1022766/150‑1 (MQ=255)
                                                           cacaTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAACTAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAgat                                                                                            <  2:4322368/150‑3 (MQ=255)
                                                           cacaTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCt                                                                                            <  1:1858625/150‑1 (MQ=255)
                                                              aTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCtata                                                                                         >  1:4320785/1‑150 (MQ=255)
                                                              aTCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCtata                                                                                         <  1:2902190/150‑1 (MQ=255)
                                                                cTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATAtt                                                                                       <  2:3143742/150‑1 (MQ=255)
                                                                cTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATAtt                                                                                       >  1:504386/1‑150 (MQ=255)
                                                                  tATCCGCTCAACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAAAAAGGCGATGATGCGCCGTTTGAAGTCGGCAATATTAAGACGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTc                                                                                     >  2:4412235/1‑150 (MQ=255)
                                                                   aTCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCa                                                                                    >  2:2186486/1‑150 (MQ=255)
                                                                        ctctACGGGATACACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTccc                                                                               <  2:2663827/150‑1 (MQ=255)
                                                                         tctACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCa                                                                              >  2:4275473/1‑150 (MQ=255)
                                                                              gggATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCCGCATg                                                                         >  1:1269476/1‑150 (MQ=255)
                                                                                 aTAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGccc                                                                      >  2:631736/1‑150 (MQ=255)
                                                                                 aTAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGccc                                                                      >  1:354904/1‑150 (MQ=255)
                                                                                  tAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCt                                                                     <  2:1014806/150‑1 (MQ=255)
                                                                                   aaaCGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTg                                                                    >  1:2807995/1‑150 (MQ=255)
                                                                                   aaaCGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTg                                                                    <  1:2855807/150‑1 (MQ=255)
                                                                                    aaCGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCTGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGa                                                                   >  2:3300873/1‑150 (MQ=255)
                                                                                        cTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTAt                                                               <  1:4234602/150‑1 (MQ=255)
                                                                                        cTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTAt                                                               <  2:900336/150‑1 (MQ=255)
                                                                                          gCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCg                                                             <  1:1124693/150‑1 (MQ=255)
                                                                                              tatCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCtat                                                         <  2:787184/150‑1 (MQ=255)
                                                                                                 ccGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAATTCGGCAATATTAAGCCGCATGCCATCACGACAAGCGGCTTATACGGTTTACCAGCTATATTTCACACCCAGCATGACCTGAGTAACGCTATAGc                                                      >  1:4085007/1‑150 (MQ=255)
                                                                                                  cgcgATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCTACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGcc                                                     >  2:1298217/1‑150 (MQ=255)
                                                                                                  cgcgATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGcc                                                     <  1:3065486/150‑1 (MQ=255)
                                                                                                    cgATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCtt                                                   >  2:3004463/1‑150 (MQ=255)
                                                                                                          cgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCAc                                             >  2:1941821/1‑150 (MQ=255)
                                                                                                          cgGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCAc                                             >  1:2666853/1‑150 (MQ=255)
                                                                                                            gATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCt                                           >  2:2560002/1‑150 (MQ=255)
                                                                                                             aTAGAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTa                                          <  1:901359/150‑1 (MQ=255)
                                                                                                             aTAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTa                                          <  1:3322306/150‑1 (MQ=255)
                                                                                                             aTAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTa                                          <  1:437515/150‑1 (MQ=255)
                                                                                                              tAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAg                                         >  2:17911/1‑150 (MQ=255)
                                                                                                               aaaaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGt                                        >  1:1906759/1‑150 (MQ=255)
                                                                                                               aaaaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGt                                        >  2:1197874/1‑150 (MQ=255)
                                                                                                               aaaaaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGt                                        >  1:3152033/1‑150 (MQ=255)
                                                                                                                  aaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGc                                     >  1:1943890/1‑150 (MQ=255)
                                                                                                                  aaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGc                                     >  1:4270809/1‑150 (MQ=255)
                                                                                                                  aaGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGc                                     >  2:1007049/1‑150 (MQ=255)
                                                                                                                    ggaggaTTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGcac                                   >  2:2392275/1‑150 (MQ=255)
                                                                                                                      aggaTTTAAATAAGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCacac                                 <  2:2873484/150‑1 (MQ=255)
                                                                                                                         aTTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGa                              <  2:268420/150‑1 (MQ=255)
                                                                                                                            tAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACAt                           <  1:1003910/150‑1 (MQ=255)
                                                                                                                               aTAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCAAGCCATCTAGACATGCGGCTTATACGGTATACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGACTTTATCACCTAGTTGCACACCGACATTcc                        >  1:3127988/1‑150 (MQ=255)
                                                                                                                                 aaGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTcccc                      <  1:375749/150‑1 (MQ=255)
                                                                                                                                  aGGCGTTGATGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTccccc                     >  1:806607/1‑150 (MQ=255)
                                                                                                                                   ggCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCca                    >  1:3662923/1‑150 (MQ=255)
                                                                                                                                   ggCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCca                    >  1:2467843/1‑150 (MQ=255)
                                                                                                                                     cGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCcaca                  <  1:2715909/150‑1 (MQ=255)
                                                                                                                                     cGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCcaca                  <  1:4523072/150‑1 (MQ=255)
                                                                                                                                      gTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAg                 >  2:3752358/1‑150 (MQ=255)
                                                                                                                                         gATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCt              <  2:2092850/150‑1 (MQ=255)
                                                                                                                                          aTGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTAAATTTCACTCCCAGCAAGCCCTGAGTATCGCTATAGCCTTAATCACATAGTTGCACACCGACATTCCCCCACAGGCTa             >  2:480768/1‑150 (MQ=255)
                                                                                                                                           tGCGCCGTTTTAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTaa            >  2:1921530/1‑150 (MQ=255)
                                                                                                                                           tGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTaa            >  1:3294202/1‑150 (MQ=255)
                                                                                                                                           tGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTaa            <  1:3239092/150‑1 (MQ=255)
                                                                                                                                             cgcCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAgg          >  2:3408018/1‑150 (MQ=255)
                                                                                                                                                     gAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTAtt  >  1:2002976/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCCGCTTAAGGGGATATTATAAGTTTTATCCTTTAGTGAGGATAAGTCAATCTGGAAACACATCTTATCCGCTCTACGGGATAAACGCTGCGATATCCGCGATCGCGGATAAAAAGGAGGATTTAAATAAGGCGTTGATGCGCCGTTTGAAGTCGGCAATATTAAGCCGCATGCCATCTCGACATGCGGCTTATACGGTTTACCAGCTATATTTCACTCCCAGCATGCCCTGAGTATCGCTATAGCCTTTATCACCTAGTTGCACACCGACATTCCCCCACAGGCTAAGGTTGTTATT  >  CP009273/1586722‑1587020

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: