Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,139,680 C→A 8.5% Q162H (CAG→CAT)  PP_5481 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,139,6800CA 8.5% 251.4 / 8.6 82Q162H (CAG→CAT) PP_5481hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/41);  new base A (3/4);  total (37/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

TAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCTCGTC  >  NC_002947/2139534‑2139826
                                                                                                                                                  |                                                                                                                                                   
tAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGc                                                                                                                                                  <  1:1415426/150‑1 (MQ=255)
tAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGc                                                                                                                                                  <  2:815372/150‑1 (MQ=255)
  ggggTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCg                                                                                                                                                <  2:1227502/150‑1 (MQ=255)
     gtgtCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcc                                                                                                                                             >  1:1363727/1‑150 (MQ=255)
        tCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCAGTGGCGTCGTCGCCAATGCCggtgcaa                                                                                                                                          <  2:1110002/150‑6 (MQ=255)
          ggATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCGCTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccc                                                                                                                                        >  2:1548730/1‑150 (MQ=255)
          ggATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccc                                                                                                                                        >  2:323556/1‑150 (MQ=255)
           gagatCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTAAGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGcccaccca                                                                                                                                       >  1:1056480/4‑150 (MQ=255)
            atatCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCccacccac                                                                                                                                      <  2:574379/150‑1 (MQ=255)
                cAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAg‑c                                                                                                                                 <  1:1236743/150‑1 (MQ=255)
                        gCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACCCCACTGCCGGCCCACCCACCAG‑CGCCAGTAg                                                                                                                         >  2:773336/1‑150 (MQ=255)
                          aCGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGaa                                                                                                                       >  2:214338/1‑150 (MQ=255)
                              cGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACAtg                                                                                                                   <  1:323556/150‑1 (MQ=255)
                              cGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACAtg                                                                                                                   <  1:1168265/150‑1 (MQ=255)
                                  cACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCCCCAG‑CGCCAGTAGAACAtggtgg                                                                                                               <  1:1124733/150‑1 (MQ=255)
                                     cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGcc                                                                                                                              <  1:1505130/132‑1 (MQ=255)
                                     cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGcc                                                                                                                              >  2:1505130/1‑132 (MQ=255)
                                     cTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGgaga                                                                                                            >  2:1000383/1‑150 (MQ=255)
                                       tCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAg                                                                                                          <  2:342474/150‑1 (MQ=255)
                                             tggtggCCACTGCGGTGCTCTACGTTGTNTGCGTTCGCTTGCTGNTNANCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCa                                                                                                    >  1:906213/1‑150 (MQ=255)
                                              ggtggCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                                   <  1:810898/150‑1 (MQ=255)
                                                         cGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGc                                                                                        >  2:1376281/1‑150 (MQ=255)
                                                          ggTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGca                                                                                       >  2:283197/1‑150 (MQ=255)
                                                                 tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                                   >  2:592843/1‑131 (MQ=255)
                                                                 tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                                   >  1:1378628/1‑131 (MQ=255)
                                                                 tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                                   <  1:592843/131‑1 (MQ=255)
                                                                 tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCaa                                                                                                   <  2:1378628/131‑1 (MQ=255)
                                                                 tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAg                                                                                                  <  2:1004036/132‑1 (MQ=255)
                                                                 tACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAg                                                                                                  >  1:1004036/1‑132 (MQ=255)
                                                                  aCGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTa                                                                                                        <  1:81809/125‑1 (MQ=255)
                                                                  aCGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTa                                                                                                        >  2:81809/1‑125 (MQ=255)
                                                                   cGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAg‑c                                                                                                                                 <  2:172835/99‑1 (MQ=255)
                                                                   cGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAg‑c                                                                                                                                 >  1:172835/1‑99 (MQ=255)
                                                                    gttgttTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGca                                                                             >  2:614836/1‑150 (MQ=255)
                                                                         ttGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCAcacca                                                                        <  2:246279/150‑1 (MQ=255)
                                                                             gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGccc                                                                                           <  2:1608997/127‑1 (MQ=255)
                                                                             gTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGccc                                                                                           >  1:1608997/1‑127 (MQ=255)
                                                                              ttCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccgag                                                                   >  2:614073/1‑150 (MQ=255)
                                                                                      cTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGAt                                                                                               >  2:1000424/1‑114 (MQ=255)
                                                                                      cTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGAt                                                                                               <  1:1000424/114‑1 (MQ=255)
                                                                                       tGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAg‑cgc                                                                                                                               <  1:420799/81‑1 (MQ=255)
                                                                                       tGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAg‑cgc                                                                                                                               >  2:420799/1‑81 (MQ=255)
                                                                                          tCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCCGTAGAAGATGGTGGAGCAGTCAGCCCTGGCGCGCACCACCACGCCGCACACCACCCACCCGCgcaggcag                                                       >  2:430298/1‑150 (MQ=255)
                                                                                           cAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAgg                                                      <  1:1121609/150‑1 (MQ=255)
                                                                                           cAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAgg                                                      <  1:614073/150‑1 (MQ=255)
                                                                                            aGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAggg                                                     >  1:40086/1‑150 (MQ=255)
                                                                                              cGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCATCACGCCGCACACCACCGAGCCGAGCAGGAAGGGtt                                                   <  1:1548730/150‑1 (MQ=255)
                                                                                              cGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGtt                                                   <  2:1056480/150‑1 (MQ=255)
                                                                                                gCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGc                                                 <  1:776505/150‑1 (MQ=255)
                                                                                                  cgcgGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGgag                                                                                                             <  1:824231/88‑1 (MQ=255)
                                                                                                  cgcgGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGgag                                                                                                             >  2:824231/1‑88 (MQ=255)
                                                                                                   gcgGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCcac                                              <  2:906213/150‑1 (MQ=255)
                                                                                                         tGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGc                                                                                        <  1:641668/102‑1 (MQ=255)
                                                                                                         tGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGc                                                                                        >  2:641668/1‑102 (MQ=255)
                                                                                                                cTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCt                                 <  1:1378797/150‑1 (MQ=255)
                                                                                                                     gCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCg                            <  1:1000383/150‑1 (MQ=255)
                                                                                                                        ttGCACTGGCATTGGCGACGACGCCACTGACGGCCC‑CCCCCCAGCCGTTAGGAGAACATGGTGGAGAAGCCAAGGTTGATGCCCAGCATCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTg                         <  1:773336/150‑1 (MQ=255)
                                                                                                                        ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccgagccga                                                               <  2:50777/112‑1 (MQ=255)
                                                                                                                        ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCAccgagccga                                                               >  1:50777/1‑112 (MQ=255)
                                                                                                                        ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCt                                 <  1:485038/142‑1 (MQ=255)
                                                                                                                        ttGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCt                                 >  2:485038/1‑142 (MQ=255)
                                                                                                                         tgggccGGCAGTGGCGTCGTCGCCAATGCCGGTCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGa                        >  2:1513341/5‑150 (MQ=255)
                                                                                                                            aCCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATcc                     >  2:812581/1‑150 (MQ=255)
                                                                                                                             ccGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCcaca                                             <  2:1558803/125‑1 (MQ=255)
                                                                                                                             ccGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCcaca                                             >  1:1558803/1‑125 (MQ=255)
                                                                                                                              cGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTg                         <  1:1332819/144‑1 (MQ=255)
                                                                                                                              cGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTg                         >  2:1332819/1‑144 (MQ=255)
                                                                                                                               ggCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACt                  <  1:812581/150‑1 (MQ=255)
                                                                                                                               ggCATTGGCGAAGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTATAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGACACAGCGTCGACAGCTGGTCGGTGATCCACt                  <  1:712469/150‑1 (MQ=255)
                                                                                                                                gcagtGGCGTCGTCGCCAATGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACt                  <  1:550415/145‑1 (MQ=255)
                                                                                                                                gCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACt                  >  2:550415/1‑149 (MQ=255)
                                                                                                                                gCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACt                  >  2:551471/1‑149 (MQ=255)
                                                                                                                                gCAGTGGCGTCGTCGCCAATGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACt                  <  1:551471/149‑1 (MQ=255)
                                                                                                                                    tGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAg             <  1:151486/150‑1 (MQ=255)
                                                                                                                                     ggCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCCGTTGGGCATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAgg            <  1:430298/150‑1 (MQ=255)
                                                                                                                                         acgacgCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAg        <  2:1363727/150‑1 (MQ=255)
                                                                                                                                              gCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACCTGGGGGAGAAGTAGGCAAGGATGCCCAGCACCAAGCCGCACACCACCGAGCCGAGCAGGAAGGGTGGGCACAGCGACGACAGCTGGTCGGTGATCCACTGGAAGGACAGCgcgt   >  1:460597/1‑150 (MQ=255)
                                                                                                                                              gCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGGAAGGGGGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCtcgt   >  1:116171/1‑150 (MQ=255)
                                                                                                                                              gCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCtcgt   >  1:1504149/1‑150 (MQ=255)
                                                                                                                                               ccACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCtcgtc  <  2:925272/150‑1 (MQ=255)
                                                                                                                                                cACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCtc     >  1:1035085/1‑146 (MQ=255)
                                                                                                                                                cACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCtc     <  2:1035085/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TAGGGGTGTCGGATATCAGCATTGGCACGGCGATCACCTTCATGCTGGTGGCCACTGCGGTGCTCTACGTTGTTTGCGTTCGCTTGCTGGTCAGCGGCCGCGGCATGCGGGCCTGAAGCCTTGCACCGGCATTGGCGACGACGCCACTGCCGGCCCACCCACCAG‑CGCCAGTAGAACATGGTGGAGAAGTAGGCAAGGATGCCCAGCACCACGCCGCACACCACCGAGCCGAGCAGGAAGGGTTGCCACAGCGTCGACAGCTGGTCGGTGATCCACTGGAAGGTCAGCTCGTC  >  NC_002947/2139534‑2139826

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: ‑