Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,198,622 A→G 9.1% D309G (GAC→GGC) ‡ gspE → protein secretion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,198,6220AG 9.1% 199.7 / 4.8 72D309G (GAC→GGC) ‡gspEprotein secretion protein
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base G (4/3);  total (34/38)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGAACAGCCTCAACGATGGCAGCCGCAACATCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCTGCGCGACATGGGCAT  >  NC_002947/1198477‑1198769
                                                                                                                                                 |                                                                                                                                                   
cTGAACAGCCTCAACGATGGCAGCCGCAACATCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCt                                                                                                                                                 <  1:292668/150‑1 (MQ=255)
        ccTCAACGATGGCAGCCGCAACATCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGAt                                                                                                                                         >  2:727777/1‑150 (MQ=255)
        ccTCAACGATGGCAGCCGCAACATCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGAt                                                                                                                                         >  2:49548/1‑150 (MQ=255)
          ttaaCGATGGCAGCCGCAACATCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCa                                                                                                                                       >  1:1292113/3‑150 (MQ=255)
              cGATGGCAGCCGCAACATCCTCACCGTCGAGGACCCGGTTGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCt                                                                                                                                   <  1:1340747/150‑1 (MQ=255)
                           aaCATCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCcgc                                                                                                                      >  1:382556/1‑150 (MQ=255)
                              aTCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGAt                                                                                                                   <  1:1245199/150‑1 (MQ=255)
                              aTCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGAt                                                                                                                   >  1:206249/1‑150 (MQ=255)
                               tcgtCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGGTCCTGACgca                                                                                                                                           >  1:1645703/4‑123 (MQ=255)
                               tCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGGTCCTGACgca                                                                                                                                           <  2:1645703/125‑3 (MQ=255)
                               tCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGAc                                                                                                                                              <  2:244751/122‑1 (MQ=255)
                               tCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGAc                                                                                                                                              >  1:244751/1‑122 (MQ=255)
                                        tCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCACCATCCTGCGTCAGGACCCTGACGTGATCATGc                                                                                                                                    <  2:450661/123‑1 (MQ=255)
                                        tCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCACCATCCTGCGTCAGGACCCTGACGTGATCATGc                                                                                                                                    >  1:450661/1‑123 (MQ=255)
                                            ggACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCa                                                                                                     >  1:1273335/1‑150 (MQ=255)
                                              aCCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGa                                                                                                   <  1:1166659/150‑1 (MQ=255)
                                               cccGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGt                                                                                                                                            <  2:1401562/108‑1 (MQ=255)
                                               cccGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGt                                                                                                                                            >  1:1401562/1‑108 (MQ=255)
                                                         tACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAg                                                                                        <  2:986540/150‑1 (MQ=255)
                                                              cATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACt                                                                                                          >  1:974772/1‑127 (MQ=255)
                                                              cATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACt                                                                                                          <  2:974772/127‑1 (MQ=255)
                                                                    cGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATc                                                                                                                                        <  2:452202/91‑1 (MQ=255)
                                                                    cGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATc                                                                                                                                        >  1:452202/1‑91 (MQ=255)
                                                                    cGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCAc                                                                             <  1:49117/150‑1 (MQ=255)
                                                                          tGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCa                                                                       <  1:727777/150‑1 (MQ=255)
                                                                             ccAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGGCCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCt                                                                    >  2:232157/1‑150 (MQ=255)
                                                                             ccAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCgg                                                                                                                                <  2:669575/90‑1 (MQ=255)
                                                                             ccAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCgg                                                                                                                                >  1:669575/1‑90 (MQ=255)
                                                                               aaaCCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTgg                                                                  <  2:926907/150‑1 (MQ=255)
                                                                               aaaCCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTgg                                                                  >  1:1603749/1‑150 (MQ=255)
                                                                                     ccATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGCCGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGc                                                                                  <  1:448504/128‑1 (MQ=255)
                                                                                     ccATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGc                                                                                  >  2:448504/1‑128 (MQ=255)
                                                                                       aTCAACCCGCGTGCCGGGCTGACTTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCg                                                          >  2:1450892/1‑150 (MQ=255)
                                                                                       aTCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTc                                                           >  2:987475/1‑149 (MQ=255)
                                                                                       aTCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTc                                                           <  1:987475/149‑1 (MQ=255)
                                                                                       aTCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCg                                                          >  2:372186/1‑150 (MQ=255)
                                                                                           aCCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAg                                                                                                    >  2:865140/1‑104 (MQ=255)
                                                                                           aCCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAg                                                                                                    <  1:865140/104‑1 (MQ=255)
                                                                                           aCCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCt                                                      >  1:1448694/1‑150 (MQ=255)
                                                                                               gcgTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCAcc                                                                            >  1:1034547/1‑124 (MQ=255)
                                                                                               gcgTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCAcc                                                                            <  2:1034547/124‑1 (MQ=255)
                                                                                                   gCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGAc                                                        >  2:342708/1‑140 (MQ=255)
                                                                                                   gCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGAc                                                        <  1:342708/140‑1 (MQ=255)
                                                                                                   gCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACAcc                                              <  2:1292113/150‑1 (MQ=255)
                                                                                                     cGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACcaa                                            <  1:1460344/150‑1 (MQ=255)
                                                                                                      gggCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATc                                                                                                                                        <  1:708994/57‑1 (MQ=255)
                                                                                                      gggCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATc                                                                                                                                        >  2:708994/1‑57 (MQ=255)
                                                                                                       ggCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGc                                                                                  <  2:594616/110‑1 (MQ=255)
                                                                                                       ggCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGc                                                                                  >  1:594616/1‑110 (MQ=255)
                                                                                                            aCGTTTGCCAGTGGCTTGCGCGCCATCCTGCGCCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCg                                     <  2:1145625/150‑1 (MQ=255)
                                                                                                              gTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGc                                                                                  >  2:492226/1‑103 (MQ=255)
                                                                                                              gTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGc                                                                                  <  1:492226/103‑1 (MQ=255)
                                                                                                              gTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGt                                   >  1:625560/1‑150 (MQ=255)
                                                                                                              gTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGt                                   >  1:543535/1‑150 (MQ=255)
                                                                                                              gTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGAGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGt                                   >  1:800137/1‑150 (MQ=255)
                                                                                                                     aGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATAATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCAAGCCTCACCGGCCATCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCa                            <  2:171145/150‑1 (MQ=255)
                                                                                                                      gTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAg                           <  1:49548/150‑1 (MQ=255)
                                                                                                                      gTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAg                           <  2:575744/150‑1 (MQ=255)
                                                                                                                      gTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAg                           <  2:1603749/150‑1 (MQ=255)
                                                                                                                         gCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGGTCGGTACAAGTCGCAGATCGAGGTGCAGGCCAGACACCCCAAAACAATGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGa                        <  1:372186/150‑1 (MQ=255)
                                                                                                                            tGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTg                                                                <  2:1173307/107‑1 (MQ=255)
                                                                                                                            tGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTg                                                                >  1:1173307/1‑107 (MQ=255)
                                                                                                                                cgcCATCCTGCGTCAGGGCCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCAc                                                                             >  2:60694/1‑90 (MQ=255)
                                                                                                                                cgcCATCCTGCGTCAGGGCCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCAc                                                                             <  1:60694/90‑1 (MQ=255)
                                                                                                                                cgcCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGAGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCt                 <  2:800137/150‑1 (MQ=255)
                                                                                                                                 gcCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCTg                <  2:1448694/150‑1 (MQ=255)
                                                                                                                                      ccTGCGTCAGGACCCTGACGTGAACATGCCCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCTGCGCGa           >  2:1177001/1‑150 (MQ=255)
                                                                                                                                      atcacGTCAGGGTCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCAcc                                                                            >  2:1522648/5‑85 (MQ=255)
                                                                                                                                      atcacGTCAGGGTCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCAcc                                                                            <  1:1522648/81‑1 (MQ=255)
                                                                                                                                          cGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCTGCGCGACATg       <  2:1273335/150‑1 (MQ=255)
                                                                                                                                            tCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCTGCGCGACATggg     >  1:1487307/1‑150 (MQ=255)
                                                                                                                                               ggACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCTGCGCGACATGGGCAt  <  2:1131944/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGAACAGCCTCAACGATGGCAGCCGCAACATCCTCACCGTCGAGGACCCGGTGGAGTACGCCATCGCCGGCATTGGCCAAACCGCCATCAACCCGCGTGCCGGGCTGACGTTTGCCAGTGGCCTGCGCGCCATCCTGCGTCAGGACCCTGACGTGATCATGCTCGGCGAAATCCGCGATCAGGAGACTGCGCAGATCGCCGTGCAGGCCAGCCTCACCGGCCACCTGGTGCTGTCGACCTTGCACACCAACAATGCGGTAGGCGCAGTGACCCGCCTGCGCGACATGGGCAT  >  NC_002947/1198477‑1198769

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