Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,661,519 C→A 10.8% A50D (GCC→GAC)  PP_1456 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,661,5190CA10.8% 270.1 / 15.2 92A50D (GCC→GAC) PP_1456membrane protein
Reads supporting (aligned to +/- strand):  ref base C (46/36);  new base A (6/4);  total (52/40)
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.09e-01

TGTTCTATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGT  >  NC_002947/1661372‑1661668
                                                                                                                                                   |                                                                                                                                                     
tGTTCTATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGccc                                                                                                                                                     <  1:628639/150‑1 (MQ=255)
          gATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCAt                                                                                                                                           >  2:574220/1‑150 (MQ=255)
          gATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACACCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGACTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCAt                                                                                                                                           >  1:599841/1‑150 (MQ=255)
             tGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCAtaac                                                                                                                                        <  2:1627858/150‑4 (MQ=255)
               ccTACATCAGCAGCGTAGTGCTGATCAGCTACGCCTTCTCCAGTGTCCCGCATCTGGACATCATCTGGTCTGCCTGGGGTGGCCTGGTGTTCCTCCTGCGCGACATGGTGCAGAGCCGTTTTGGTCATGGCGCCCTGGTCGCTATGCtgg                                                                                                                                      <  2:925664/150‑1 (MQ=255)
               ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtgg                                                                                                                                      <  2:304974/150‑1 (MQ=255)
               ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtgg                                                                                                                                      >  2:667779/1‑150 (MQ=255)
               ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtgg                                                                                                                                      >  1:1385894/1‑150 (MQ=255)
               ccTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtgg                                                                                                                                      >  1:1916305/1‑150 (MQ=255)
                cTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCTACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggt                                                                                                                                     <  1:123365/150‑1 (MQ=255)
                cTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCTGACCCGTTTGGGTCATGGCGCCCTGGTCGCCATGCtggt                                                                                                                                     >  1:1645348/1‑150 (MQ=255)
                 tACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggtg                                                                                                                                    <  2:507113/150‑1 (MQ=255)
                   cATCAGCAGCGTAGTGCTGATCAACTACACCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGc                                                                                                                                  <  2:599841/150‑1 (MQ=255)
                     tCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCAc                                                                                                                                <  1:1174425/150‑1 (MQ=255)
                     tCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCAc                                                                                                                                <  2:1639406/150‑1 (MQ=255)
                       agcagcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTTATCGCCATGCTGGTCGCAcag                                                                                                                              <  2:805341/150‑3 (MQ=255)
                       agcagcGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTg                                                                                                                              <  1:667779/150‑1 (MQ=255)
                               aGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTc                                                                                                                      <  2:147503/150‑1 (MQ=255)
                                   cTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTAt                                                                                                                  >  1:365988/1‑150 (MQ=255)
                                     gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGccc                                                                                                                                                     <  1:1917547/113‑1 (MQ=255)
                                     gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGccc                                                                                                                                                     >  2:1917547/1‑113 (MQ=255)
                                     gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGt                                                                                                                >  1:1623624/1‑150 (MQ=255)
                                     gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGt                                                                                                                >  2:304911/1‑150 (MQ=255)
                                     gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATgac                                                                                                                                        <  2:442508/126‑3 (MQ=255)
                                     gATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGACCAGGGCGCCATgac                                                                                                                                        >  1:442508/1‑124 (MQ=255)
                                      aTCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTTGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTc                                                                                                               >  1:1888760/1‑150 (MQ=255)
                                              cGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCgcc                                                                                                       >  2:865293/1‑148 (MQ=255)
                                                 cTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGcc                                                                                                    >  2:1134678/1‑150 (MQ=255)
                                                  ttCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCg                                                                                                   >  2:392075/1‑150 (MQ=255)
                                                  ttCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCg                                                                                                   >  1:1091883/1‑150 (MQ=255)
                                                  ttCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCg                                                                                                   >  1:1092114/1‑150 (MQ=255)
                                                      ccAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCa                                                                                               >  1:938881/1‑150 (MQ=255)
                                                        aGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATc                                                                                             >  2:639896/1‑150 (MQ=255)
                                                            ccccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGccc                                                                                         >  2:878197/1‑150 (MQ=255)
                                                             cccGCATCTGGATATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTc                                                                                                                      <  2:1104893/120‑1 (MQ=255)
                                                             cccGCATCTGGATATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTc                                                                                                                      >  1:1104893/1‑120 (MQ=255)
                                                              ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTg                                                                                       >  1:558619/1‑150 (MQ=255)
                                                              ccGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTg                                                                                       >  1:1151750/1‑150 (MQ=255)
                                                                gCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGc                                                                                     <  2:1408222/150‑1 (MQ=255)
                                                                  aTCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCTACATGGTGCAGACCCGGTTTGGTCATGGCGCCCTGGTCGCCATGTTGGTAGCACGAGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGccc                                                                                   >  2:1485960/1‑149 (MQ=255)
                                                                      ggACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAgcgc                                                                               >  2:816837/1‑150 (MQ=255)
                                                                        aCATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTAt                                                                                                                  >  2:1584637/1‑113 (MQ=255)
                                                                        aCATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTAt                                                                                                                  <  1:1584637/113‑1 (MQ=255)
                                                                        aCATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCa                                                                             >  2:173675/1‑150 (MQ=255)
                                                                          atcatcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACt                                                                           >  2:1215023/1‑150 (MQ=255)
                                                                           tcatcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGTAGCCGGCCATCGCCCTGGCCAGCGCCACTg                                                                          <  1:112919/150‑1 (MQ=255)
                                                                             atcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCg                                                                        >  1:256803/1‑150 (MQ=255)
                                                                              tcTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGt                                                                       >  2:54783/1‑150 (MQ=255)
                                                                                tGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggtg                                                                                                                                    >  2:1385674/1‑87 (MQ=255)
                                                                                tGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCtggtg                                                                                                                                    <  1:1385674/87‑1 (MQ=255)
                                                                                tGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttc                                                                     <  1:574220/150‑1 (MQ=255)
                                                                                tGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGctc                                                                     >  2:772433/1‑150 (MQ=255)
                                                                                     tGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTg                                                                          <  2:1702060/140‑1 (MQ=255)
                                                                                     tGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTg                                                                          >  1:1702060/1‑140 (MQ=255)
                                                                                       ccTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATg                                                                                                                 >  2:487045/1‑99 (MQ=255)
                                                                                       ccTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATg                                                                                                                 <  1:487045/99‑1 (MQ=255)
                                                                                       ccTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCAttt                                                              >  1:1191797/1‑150 (MQ=255)
                                                                                             gcggcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATACTGGTGGCACTGGTACTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGc                                                        >  2:1038036/1‑150 (MQ=255)
                                                                                                gcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGa                                                     >  1:296782/1‑150 (MQ=255)
                                                                                                gcCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGa                                                     <  1:1419148/150‑1 (MQ=255)
                                                                                                      tgtTCACCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGAct                                               <  1:639896/150‑1 (MQ=255)
                                                                                                       gtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGc                                                                         >  1:60393/1‑123 (MQ=255)
                                                                                                       gtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGc                                                                         <  2:60393/123‑1 (MQ=255)
                                                                                                       gtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttctt                                                                   <  1:7172/129‑1 (MQ=255)
                                                                                                       gtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGttctt                                                                   >  2:7172/1‑129 (MQ=255)
                                                                                                       gtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActg                                              <  2:938881/150‑1 (MQ=255)
                                                                                                       gtTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATACTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActg                                              <  1:1038036/150‑1 (MQ=255)
                                                                                                        ttCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActgg                                             >  1:565310/1‑150 (MQ=255)
                                                                                                           aTCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                          <  1:1520613/150‑1 (MQ=255)
                                                                                                            tCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctgg                                         >  2:470061/1‑150 (MQ=255)
                                                                                                             ccggcgcGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgt                                        >  2:1716463/4‑150 (MQ=255)
                                                                                                             ccTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGgt                                        <  2:1888760/150‑1 (MQ=255)
                                                                                                              cTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTg                                                                                                                              >  1:1239774/1‑63 (MQ=255)
                                                                                                              cTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTg                                                                                                                              <  2:1239774/63‑1 (MQ=255)
                                                                                                                gcgcGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGtt                                     <  2:1916305/150‑1 (MQ=255)
                                                                                                                gcgcGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGtt                                     <  2:1385894/150‑1 (MQ=255)
                                                                                                                 cgcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTACGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTc                                    <  2:365988/150‑1 (MQ=255)
                                                                                                                  gcgACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCa                                   <  1:1162858/150‑1 (MQ=255)
                                                                                                                    gACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGc                                 >  2:1170785/1‑150 (MQ=255)
                                                                                                                     aCATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCa                                >  2:566536/1‑150 (MQ=255)
                                                                                                                      cATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGa                                                     <  2:1599142/128‑1 (MQ=255)
                                                                                                                      cATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGa                                                     >  1:1599142/1‑128 (MQ=255)
                                                                                                                      cATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCa                                >  1:193252/1‑149 (MQ=255)
                                                                                                                      cATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCa                                <  2:193252/149‑1 (MQ=255)
                                                                                                                          gTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCTAGCTGATCGACTGGCTGGTGTTCAGCATCAcc                           <  1:392075/150‑1 (MQ=255)
                                                                                                                                       ttGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTgc              <  2:256803/150‑1 (MQ=255)
                                                                                                                                       accagcATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGc                                                        >  2:1038658/6‑108 (MQ=255)
                                                                                                                                       accagcATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGc                                                        <  1:1038658/103‑1 (MQ=255)
                                                                                                                                         ggTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTgcgc            <  1:173675/150‑1 (MQ=255)
                                                                                                                                         cagcATGTCGACCAGGTCTCCATGCTGGTTGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGATCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCATCCGCTgcgc            >  2:910401/4‑150 (MQ=255)
                                                                                                                                            cATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                          <  2:1265704/117‑1 (MQ=255)
                                                                                                                                            cATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggctg                                          >  1:1265704/1‑117 (MQ=255)
                                                                                                                                            cATGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAc         >  2:1666567/1‑150 (MQ=255)
                                                                                                                                             aTGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCTGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAcc        >  1:189117/1‑150 (MQ=255)
                                                                                                                                             aTGGCGACCAGGGCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGAcc        >  1:559180/1‑150 (MQ=255)
                                                                                                                                                gcgcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggc                                            >  2:8401/1‑111 (MQ=255)
                                                                                                                                                gcgcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGActggc                                            <  1:8401/111‑1 (MQ=255)
                                                                                                                                                 cgcCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCt    >  2:838014/1‑150 (MQ=255)
                                                                                                                                                   cccTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCtgt  <  2:1623624/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGTTCTATCTGATTGCCTACATCAGCAGCGTAGTGCTGATCAACTACGCCTTCTCCAGTGCCCCGCATCTGGACATCATCTGGTCTGCCTGGGGCGGCCTGGTGTTCATCCTGCGCGACATGGTGCAGACCCGTTTTGGTCATGGCGCCCTGGTCGCCATGCTGGTGGCACTGGTGCTGTCCTATGTCACCTCGGAGCCGGCCATCGCCCTGGCCAGCGCCACTGCGTTCTTCATTTCCGAGCTGATCGACTGGCTGGTGTTCAGCATCACCCGTCGGCCGCTGCGCGACCGCCTGT  >  NC_002947/1661372‑1661668

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: