Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,575,198 G→A 6.2% intergenic (‑184/+26) pcaP ← / ← galP‑II porin/porin‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,575,1980GA 6.2% 290.8 / 5.8 96intergenic (‑184/+26)pcaP/galP‑IIporin/porin‑like protein
Reads supporting (aligned to +/- strand):  ref base G (45/45);  new base A (2/4);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

ACCTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCC  >  NC_002947/1575050‑1575341
                                                                                                                                                    |                                                                                                                                               
aCCTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCggg                                                                                                                                                 >  7:144460/1‑149 (MQ=255)
  cTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGtt                                                                                                                                               >  6:40907/1‑149 (MQ=255)
   tGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGtt                                                                                                                                               >  6:56805/1‑148 (MQ=255)
   tGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGtt                                                                                                                                               >  6:51214/1‑148 (MQ=255)
     aTAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTAt                                                                                                                                            >  8:162964/1‑149 (MQ=255)
     aTAGTGATAAGGATAGCCTTCATTCACCACAATCATCATCTCAGCGAATTGCGTCAGCTCGCACTTGGCAATTCATGGCAGGGTTTTCATCAAGCTGTTGCAAATGATGCCTGCAAAAGGGAGGGAGTCTTACCTATTCTCGGGTTTa                                                                                                                                             >  5:188257/1‑148 (MQ=255)
         tGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCcgc                                                                                                                                        <  1:249210/149‑1 (MQ=255)
           aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCcgcg                                                                                                                                       >  5:1438/1‑148 (MQ=255)
           aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCcgcg                                                                                                                                       >  2:322865/1‑148 (MQ=255)
           aTCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATAcgcg                                                                                                                                       >  7:308191/1‑148 (MQ=255)
             cAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCcgcg                                                                                                                                       >  2:331857/1‑146 (MQ=255)
              aGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAag                                                                                                                                    <  8:144460/148‑1 (MQ=255)
               ggATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAatag                                                                                                                                  <  4:349859/149‑1 (MQ=255)
               ggATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAAATCATCGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAagag                                                                                                                                  >  6:83749/1‑149 (MQ=255)
                 aTAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                <  2:10363/149‑1 (MQ=255)
                    gCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCaa                                                                                                                             <  5:40907/149‑1 (MQ=255)
                     ccTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAt                                                                                                                            >  6:9768/1‑149 (MQ=255)
                      cTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTATCATCACGCTGTTGCCAATGATGCCTTCAAAAGGCAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCaaat                                                                                                                           >  1:337321/1‑147 (MQ=255)
                       tGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATTCGTTCAAAAGGTAGGGCGTCTGGCCTATTGGCGGGt                                                                                                                                                <  5:374401/127‑1 (MQ=255)
                       tGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGt                                                                                                                                                >  6:374401/1‑127 (MQ=255)
                       tGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAAtt                                                                                                                           >  6:167807/1‑148 (MQ=255)
                         cAGTTACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAg                                                                                                                        <  4:331090/149‑1 (MQ=255)
                         cAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAg                                                                                                                        <  5:9768/149‑1 (MQ=255)
                            tCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         >  5:152016/1‑145 (MQ=255)
                            tCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         <  6:152016/145‑1 (MQ=255)
                            tCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTca                                                                                                                     >  2:343517/1‑149 (MQ=255)
                            tCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTca                                                                                                                     >  4:164650/1‑149 (MQ=255)
                            tCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTca                                                                                                                     <  7:162964/149‑1 (MQ=255)
                              accaccATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         <  2:194673/143‑1 (MQ=255)
                              accaccATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         >  1:194673/1‑143 (MQ=255)
                                caccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACag                                                                                                                  >  2:332439/1‑148 (MQ=255)
                                  ccaTCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGATGAGGCCAATTGAGTCACagaag                                                                                                               <  7:239042/149‑1 (MQ=255)
                                   catcatCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaaga                                                                                                              <  2:109270/149‑1 (MQ=255)
                                   catcatCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACagaaga                                                                                                              <  7:271901/149‑1 (MQ=255)
                                   aatcatcAGCTCAGCGAATGGCGGAAGCTCGCACTGGGCAATTCATGGAAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGTGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGTCCAATTGAGTCACagaaga                                                                                                              <  7:223647/148‑1 (MQ=255)
                                       atcaGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGaata                                                                                                                                   <  6:350427/124‑3 (MQ=255)
                                       atcaGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCAAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGaata                                                                                                                                   >  5:350427/1‑122 (MQ=255)
                                         cagctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                >  8:74959/1‑125 (MQ=255)
                                         cagctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGc                                                                                                                                <  7:74959/125‑1 (MQ=255)
                                         cagctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTcac                                                                                                                    >  6:215194/1‑137 (MQ=255)
                                         cagctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTcac                                                                                                                    <  5:215194/137‑1 (MQ=255)
                                            ctcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGa                                                                                                     >  4:85444/1‑149 (MQ=255)
                                             tcagcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGa                                                                                                     >  2:363025/1‑148 (MQ=255)
                                                gcGAATGGCGGCAGCTCGCACTTGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAAc                                                                                                 <  5:51214/149‑1 (MQ=255)
                                                gcGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAAc                                                                                                 <  5:56805/149‑1 (MQ=255)
                                                   aaTGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTca                                                                                              <  6:1438/149‑1 (MQ=255)
                                                       gcggcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         >  6:78195/1‑118 (MQ=255)
                                                       gcggcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         <  5:78195/118‑1 (MQ=255)
                                                        cggcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATc                                                                                         <  3:169399/149‑1 (MQ=255)
                                                          gcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTCTTCGCGGATAAACCCGCGAATAgg                                                                                                                                 <  8:306635/107‑1 (MQ=38)
                                                          gcAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCaa                                                                                       <  1:343517/149‑1 (MQ=255)
                                                               tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGt                                                                                  >  6:347030/1‑149 (MQ=255)
                                                               tCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGt                                                                                  >  8:174903/1‑149 (MQ=255)
                                                                 gCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTtc                                                                                >  5:327690/1‑149 (MQ=255)
                                                                 gCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTtc                                                                                >  6:91407/1‑149 (MQ=255)
                                                                  cACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTtct                                                                               >  8:168832/1‑149 (MQ=255)
                                                                    cTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCa                                                                             <  1:363025/149‑1 (MQ=255)
                                                                     tGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCAt                                                                            <  3:85444/149‑1 (MQ=255)
                                                                       ggCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCAt                                                                            <  7:174903/147‑1 (MQ=255)
                                                                       ggCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATc                                                                           <  1:331857/148‑1 (MQ=255)
                                                                               aTGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGtt                                                                  >  1:50399/1‑149 (MQ=255)
                                                                               aTGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGtc                                                                  <  3:164650/149‑2 (MQ=255)
                                                                                 ggCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGtt                                                                                 >  6:389597/1‑132 (MQ=255)
                                                                                 ggCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGtt                                                                                 <  5:389597/132‑1 (MQ=255)
                                                                                 ggCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGtt                                                                                 >  6:124897/1‑132 (MQ=255)
                                                                                 ggCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGtt                                                                                 <  5:124897/132‑1 (MQ=255)
                                                                                  gCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGt                                                               >  5:154345/1‑149 (MQ=255)
                                                                                   cAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         >  2:113585/1‑90 (MQ=255)
                                                                                   cAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGa                                                                                                                         <  1:113585/90‑1 (MQ=255)
                                                                                        ttttCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                         >  3:367226/1‑149 (MQ=255)
                                                                                        ttttCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                         <  2:50399/149‑1 (MQ=255)
                                                                                        ttttCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGaa                                                         <  7:168832/149‑1 (MQ=255)
                                                                                          ttCATACCGCTTTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAg                                                        <  8:308191/148‑1 (MQ=255)
                                                                                            catcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGt                                                                                                                       >  5:380041/1‑83 (MQ=255)
                                                                                            catcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGt                                                                                                                       <  6:380041/83‑1 (MQ=255)
                                                                                             atcaCGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTcgc                                                    >  7:94409/1‑149 (MQ=255)
                                                                                                  gCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACg                                                >  7:349949/1‑148 (MQ=255)
                                                                                                      ttGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCAt                                             <  8:349949/147‑1 (MQ=255)
                                                                                                               atgCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                  <  5:347030/149‑1 (MQ=255)
                                                                                                                tgCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgc                                  <  5:91407/148‑1 (MQ=255)
                                                                                                                 gCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgcgc                                >  8:191448/1‑149 (MQ=255)
                                                                                                                 gCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTgcgc                                >  2:313140/1‑149 (MQ=255)
                                                                                                                        aaaGGGAGGGCGTCTGGCCTATTCGGGGGTTTATCCGCGAAGAGGCCAATTGAGGCACTGTATATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACTGGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCa                         <  2:337321/149‑1 (MQ=255)
                                                                                                                        aaaGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCa                         <  1:322865/149‑1 (MQ=255)
                                                                                                                        aaaGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCa                         <  5:167807/149‑1 (MQ=255)
                                                                                                                         aaGGGAGGGCGGCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAg                        <  5:83749/149‑1 (MQ=255)
                                                                                                                          agggagggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGc                       <  6:154345/149‑1 (MQ=255)
                                                                                                                              agggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                             <  1:334062/139‑1 (MQ=255)
                                                                                                                              agggCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCAt                             >  2:334062/1‑139 (MQ=255)
                                                                                                                                          ctcttCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGc                       <  1:140774/130‑1 (MQ=255)
                                                                                                                                          ctcttCGCGGATAAACCCGCGAATAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGc                       >  2:140774/4‑133 (MQ=255)
                                                                                                                                          cTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGc       >  3:360606/1‑149 (MQ=255)
                                                                                                                                           tATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTg           <  5:8028/144‑1 (MQ=255)
                                                                                                                                           tATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTg           >  6:8028/1‑144 (MQ=255)
                                                                                                                                              tCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGGAGTCGCGACGCATGGTCGAGTTGCGCCATTTc                          >  3:126767/1‑126 (MQ=255)
                                                                                                                                              tCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGGAGTCGCGACGCATGGTCGAGTTGCGCCATTTc                          <  4:126767/126‑1 (MQ=255)
                                                                                                                                               cgcgGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGccc  >  5:39925/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACCTGATAGTGATCAGGATAGCCTGCAGTCACCACCATCATCAGCTCAGCGAATGGCGGCAGCTCGCACTGGGCAATTCATGGCAGGGTTTTCATCACGCTGTTGCCAATGATGCCTTCAAAAGGGAGGGCGTCTGGCCTATTCGCGGGTTTATCCGCGAAGAGGCCAATTGAGTCACAGAAGATTCAGCGGATAACTCACAATCAAGCGGTTCTCATCGAACGCGTTGGTGCTGAAGTCGCGACGCATGGTCGAGTTGCGCCATTTCAGCGACAGGTCCTTGAAGCTGCCC  >  NC_002947/1575050‑1575341

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑