Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,811,276 A→T 11.7% D78V (GAC→GTC)  dsbG → protein disulfide isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,811,2760AT11.7% 239.0 / 17.8 84D78V (GAC→GTC) dsbGprotein disulfide isomerase
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base T (5/5);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 7.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

CAGCTGCAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTACCTGTTCAGCGAC  >  NC_002947/4811131‑4811424
                                                                                                                                                 |                                                                                                                                                    
cagcTGCAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTTGACGaa                                                                                                                                                  <  2:1198738/150‑1 (MQ=255)
   cTGCAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAg                                                                                                                                               >  1:1285257/1‑150 (MQ=255)
    tGCAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGTCGAACAgg                                                                                                                                              <  1:762654/150‑1 (MQ=255)
    tGCAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAgg                                                                                                                                              >  2:2061092/1‑150 (MQ=255)
     gCAAGCCAAGGGCTCCGCGATCAAGGGCAGCTTCGATGCGCCCAACGGCGTGGACGGTTAACCCGCCGAGTACGAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGTCAGCCTGTTCGTCGAACAggc                                                                                                                                             <  1:281214/150‑2 (MQ=255)
     gCAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAggg                                                                                                                                             <  2:2311930/150‑1 (MQ=255)
      cAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGTCGAACAggct                                                                                                                                            <  2:1524777/150‑3 (MQ=255)
      cAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGc                                                                                                                                            <  2:878341/150‑1 (MQ=255)
         gCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAg                                                                                                                                         <  2:1305246/150‑1 (MQ=255)
          ccAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGCCGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAgg                                                                                                                                        >  2:353131/1‑150 (MQ=255)
          ccAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAgg                                                                                                                                        <  2:1651169/150‑1 (MQ=255)
           cAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGa                                                                                                                                       >  2:1126966/1‑150 (MQ=255)
            aaGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGAc                                                                                                                                      <  1:764507/150‑1 (MQ=255)
                gcgcCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGa                                                                                                                                  >  1:1914689/1‑150 (MQ=255)
                            aGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGc                                                                                                                      >  2:1675404/1‑150 (MQ=255)
                            aGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGc                                                                                                                      <  2:682533/150‑1 (MQ=255)
                                     tCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGc                                                                                                             <  2:923396/150‑1 (MQ=255)
                                       gATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGTCCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTg                                                                                                           <  2:2423052/150‑1 (MQ=255)
                                       gATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGCCCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGAAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTg                                                                                                           <  1:1250722/150‑1 (MQ=255)
                                        aTGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTgg                                                                                                          <  2:1651236/150‑1 (MQ=255)
                                         tGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGt                                                                                                         >  2:1489120/1‑150 (MQ=255)
                                               cAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTcgac                                                                                                                                                     >  1:2216084/1‑100 (MQ=255)
                                               cAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTcgac                                                                                                                                                     <  2:2216084/100‑1 (MQ=255)
                                                 aCGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACgcgc                                                                                                 <  1:1771752/150‑1 (MQ=255)
                                                 aCGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACgcgc                                                                                                 <  1:86825/150‑1 (MQ=255)
                                                    gCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACGGGCGGGGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGccag                                                                                              >  1:1451970/1‑148 (MQ=255)
                                                       tGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGa                                                                                           <  2:126463/150‑1 (MQ=255)
                                                       tGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGa                                                                                           <  2:1712563/150‑1 (MQ=255)
                                                       tGCGCGGTTATGCCGCCGAGTAACAGAACAACGCACTGCCCCAATACATGATGACCGACGGAAAGCGGATGCTAGCCTGCAGCCGGTTTTTCGGACGGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGa                                                                                           >  2:866988/1‑150 (MQ=255)
                                                           cgattAGGCCGCCGATGACCACAACCACGCCCTGGCCCTGTACCTGACCCCCGAGGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTCCAAAAGCTGGTATACGCGCCGATGAGCaa                                                                                       <  2:817415/147‑1 (MQ=255)
                                                                 tGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTcga                                                                                                                                                      <  2:1699491/81‑1 (MQ=255)
                                                                 tGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTcga                                                                                                                                                      >  1:1699491/1‑81 (MQ=255)
                                                                 tGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAAt                                                                                 >  1:73750/1‑150 (MQ=255)
                                                                   ccgccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGTCGAACAggc                                                                                                                                             >  2:824923/1‑87 (MQ=255)
                                                                   ccgccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGTCGAACAggc                                                                                                                                             <  1:824923/88‑2 (MQ=255)
                                                                    cgcngAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTg                                                                              >  2:1501152/1‑150 (MQ=255)
                                                                    cgccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTg                                                                              >  2:2382905/1‑150 (MQ=255)
                                                                    cgccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTg                                                                              >  1:94134/1‑150 (MQ=255)
                                                                     gccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCaa                                                                                                                                          <  2:861116/89‑1 (MQ=255)
                                                                     gccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCaa                                                                                                                                          >  1:861116/1‑89 (MQ=255)
                                                                     gccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTgg                                                                             <  1:1715874/150‑1 (MQ=255)
                                                                     gccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTgg                                                                             <  2:1416731/150‑1 (MQ=255)
                                                                     gccgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTgg                                                                             <  1:1126966/150‑1 (MQ=255)
                                                                       cgAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGc                                                                           >  1:1006062/1‑150 (MQ=255)
                                                                        gAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTTTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACgcgc                                                                                                 >  1:1994173/1‑127 (MQ=255)
                                                                        gAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACgcgc                                                                                                 <  2:1994173/127‑1 (MQ=255)
                                                                            aCCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGaaaa                                                                      <  1:827354/150‑1 (MQ=255)
                                                                             ccAGAACAACGCCCTTGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAAt                                                                     >  1:322169/1‑150 (MQ=255)
                                                                              cAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATg                                                                    >  2:1445941/1‑150 (MQ=255)
                                                                                    aacGCCCTGGCCATGTACGTGACCCCCTACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGaaaaa                                                              >  2:2040482/1‑150 (MQ=255)
                                                                                        ccctggccctgTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACcg                                                          >  1:565897/1‑150 (MQ=255)
                                                                                        ccctggccctgTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACcg                                                          <  2:1349812/150‑1 (MQ=255)
                                                                                           gggccctgTACCGGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGt                                                       <  2:203304/149‑1 (MQ=255)
                                                                                              ccctgTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACcgc                                                         <  1:753217/145‑1 (MQ=255)
                                                                                              ccctgTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACcgc                                                         >  2:753217/1‑145 (MQ=255)
                                                                                                ctgTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTgg                                                                                                          >  1:1065971/1‑94 (MQ=255)
                                                                                                ctgTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTgg                                                                                                          <  2:1065971/94‑1 (MQ=255)
                                                                                                   tACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGcc                                               <  1:2382905/150‑1 (MQ=255)
                                                                                                   tACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGcc                                               <  2:73750/150‑1 (MQ=255)
                                                                                                     ccTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGa                                             <  2:2001773/150‑1 (MQ=255)
                                                                                                       tGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACg                                           >  2:1396589/1‑150 (MQ=255)
                                                                                                       tGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACg                                           <  2:909020/150‑1 (MQ=255)
                                                                                                             ccGACGGAGAGAAGGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAgg                                     >  2:99237/1‑150 (MQ=255)
                                                                                                              cGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGc                                    >  1:837423/1‑150 (MQ=255)
                                                                                                              cGACGGCAACCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGc                                    <  1:1563443/150‑1 (MQ=255)
                                                                                                                 cGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGa                                 <  1:1396589/150‑1 (MQ=255)
                                                                                                                   gCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATg                               >  1:1661046/1‑150 (MQ=255)
                                                                                                                    cAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATgc                              <  1:1489120/150‑1 (MQ=255)
                                                                                                                     aaGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGAt                                                                                             >  2:1848439/1‑86 (MQ=255)
                                                                                                                     aaGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGAt                                                                                             <  1:1848439/86‑1 (MQ=255)
                                                                                                                         aCGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGcgccgc                         >  1:2089818/1‑150 (MQ=255)
                                                                                                                           gTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCgcgc                       >  1:312585/1‑150 (MQ=255)
                                                                                                                              cTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAg                    <  2:94134/150‑1 (MQ=255)
                                                                                                                                 gTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAAt                                                                     >  1:2206857/1‑98 (MQ=255)
                                                                                                                                 gTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAAt                                                                     <  2:2206857/98‑1 (MQ=255)
                                                                                                                                    ggCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTAc              <  1:2040482/150‑1 (MQ=255)
                                                                                                                                    ggCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTAc              <  2:393279/150‑1 (MQ=255)
                                                                                                                                    ggCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTAc              <  2:2089818/150‑1 (MQ=255)
                                                                                                                                     gCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTAcc             >  2:116662/1‑150 (MQ=255)
                                                                                                                                        cccTGTTCGTCGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTAc              <  1:2050250/145‑1 (MQ=255)
                                                                                                                                        cccTGTTCGTCGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTAc              >  2:2050250/2‑146 (MQ=255)
                                                                                                                                         ccTGTTCGTCGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTACCTGtt         >  2:758708/1‑150 (MQ=255)
                                                                                                                                         ccTGTTCGTCGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGAATGGGGTACCTGtt         >  2:1636487/1‑150 (MQ=255)
                                                                                                                                           tGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGgtg                 <  1:528536/140‑1 (MQ=255)
                                                                                                                                           tGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGgtg                 >  2:528536/1‑140 (MQ=255)
                                                                                                                                            gTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGa                                 <  1:1743877/123‑1 (MQ=255)
                                                                                                                                            gTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGa                                 >  2:1743877/1‑123 (MQ=255)
                                                                                                                                             ttCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTACCTGTTCAGc     >  1:1068013/1‑150 (MQ=255)
                                                                                                                                               cgacgaACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTACCTGTTCAGCGa   <  2:1562783/150‑1 (MQ=255)
                                                                                                                                                gacgaACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTACCTGTTCAGc     >  1:1374231/1‑147 (MQ=255)
                                                                                                                                                gacgaACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTACCTGTTCAGc     <  2:1374231/147‑1 (MQ=255)
                                                                                                                                                gacgaACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGGACCTGTTCAGCGAc  >  1:492168/1‑150 (MQ=255)
                                                                                                                                                gacgaACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCCGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCTCGCAAGGTGTACCTGTTCAGCGAc  <  1:65162/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CAGCTGCAAGCCAAGGGCGCCGTGATCAAGGGCAGCTTCGATGCGCCCAACGGCCTGCGCGGTTATGCCGCCGAGTACCAGAACAACGCCCTGGCCCTGTACCTGACCCCCGACGGCAAGCACGTGCTGGTCGGCAGCCTGTTCGACGAACAGGGCAAGGACCTGAGCGCCGAGCCGCTGCAAAAGCTGGTATACGCGCCGATGAGCAAGGAAATCTGGGCGAAAATGGAAAAAACCGCGTGGATCGCCGACGGCAAGGCCGATGCGCCGCGCAAGGTGTACCTGTTCAGCGAC  >  NC_002947/4811131‑4811424

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