Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,083,532 C→T 5.9% A48V (GCG→GTG)  PP_1862 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,083,5320CT 5.9% 251.5 / 3.4 84A48V (GCG→GTG) PP_1862membrane protein
Reads supporting (aligned to +/- strand):  ref base C (40/39);  new base T (2/3);  total (42/42)
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.98e-01

CTGATGATCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGATGCGC  >  NC_002947/2083387‑2083681
                                                                                                                                                 |                                                                                                                                                     
cTGATGATCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTg                                                                                                                                                   <  1:1085035/150‑1 (MQ=255)
cTGATGATCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTg                                                                                                                                                   <  1:250697/150‑1 (MQ=255)
  gattatCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGt                                                                                                                                                 <  1:1052094/150‑1 (MQ=255)
   atgatCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTg                                                                                                                                                >  2:1263418/1‑150 (MQ=255)
    tgatCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGa                                                                                                                                               <  2:1635856/150‑1 (MQ=255)
     gatCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAt                                                                                                                                              <  2:994215/150‑1 (MQ=255)
     gatCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAt                                                                                                                                              <  1:1263418/150‑1 (MQ=255)
      atCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAtg                                                                                                                                             >  2:974554/1‑150 (MQ=255)
         gAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAtggtg                                                                                                                                          <  2:1674812/150‑1 (MQ=255)
         gAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAtggtg                                                                                                                                          <  1:375449/150‑1 (MQ=255)
           gTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATggtggt                                                                                                                                        <  2:480805/150‑1 (MQ=255)
            tGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGgtggtg                                                                                                                                       <  1:681354/150‑1 (MQ=255)
              gttgttGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGcc                                                                                                                                     >  1:214967/1‑150 (MQ=255)
              gttgttGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGcc                                                                                                                                     <  1:974554/150‑1 (MQ=255)
              gttgttGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGcc                                                                                                                                     >  2:1009440/1‑150 (MQ=255)
               ttgttgTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCg                                                                                                                                    >  1:2403301/1‑150 (MQ=255)
               ttgttgTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATAGCGCCGCTGGGCGTGCTCTTCGGCCTGGAGCCGAAGTCGGATGCGGGGTCCGCGTTGGTGATGGTGGTGCCg                                                                                                                                    >  1:1622063/1‑150 (MQ=255)
                 gttgtATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACGAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGaa                                                                                                                                  >  1:1360150/1‑150 (MQ=255)
                  ttgtATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGc                                                                                                                                      >  1:1857994/1‑145 (MQ=255)
                  ttgtATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGc                                                                                                                                      <  2:1857994/145‑1 (MQ=255)
                  ttgtATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAAc                                                                                                                                 >  1:590788/1‑150 (MQ=255)
                  ttgtATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAAc                                                                                                                                 >  2:2345595/1‑150 (MQ=255)
                     tataTCGTGCTGGGCGAAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGAGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGtg                                                                                                                              >  1:55254/1‑150 (MQ=255)
                       tatCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGtgct                                                                                                                            >  1:1563227/1‑150 (MQ=255)
                         tCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGtgctgc                                                                                                                          >  2:1437375/1‑150 (MQ=255)
                           gTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGctgctg                                                                                                                        >  1:972176/1‑150 (MQ=255)
                            tGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAtggtg                                                                                                                                          >  1:1360496/1‑131 (MQ=255)
                            tGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAtggtg                                                                                                                                          <  2:1360496/131‑1 (MQ=255)
                                     cAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGc                                                                                                              <  1:1357863/150‑1 (MQ=255)
                                     cAGCCCTAGGCACCATGGCGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGc                                                                                                              <  1:211171/150‑1 (MQ=255)
                                        cccTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGt                                                                                                           >  1:2087751/1‑150 (MQ=255)
                                          cTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATggtggt                                                                                                                                        <  1:1605108/119‑1 (MQ=255)
                                          cTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATggtggt                                                                                                                                        >  2:1605108/1‑119 (MQ=255)
                                               cACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGTACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCa                                                                                                    <  2:1563227/150‑1 (MQ=255)
                                                aCCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAg                                                                                                   >  1:1430441/1‑150 (MQ=255)
                                                 ccATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCAGGCGTTTTGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTCCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAgc                                                                                                  <  1:1955164/150‑1 (MQ=255)
                                                  cATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAgcg                                                                                                 >  1:217573/1‑150 (MQ=255)
                                                     ggggggCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCaa                                                                                              >  1:1821970/1‑150 (MQ=255)
                                                       ggggCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAAcc                                                                                            >  1:1237882/1‑150 (MQ=255)
                                                        gggCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGctgctg                                                                                                                        <  2:1441340/121‑1 (MQ=255)
                                                        gggCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGctgctg                                                                                                                        >  1:1441340/1‑121 (MQ=255)
                                                              gTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATcgc                                                                                     <  2:1622063/150‑1 (MQ=255)
                                                                   gCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCtgtg                                                                                                                >  2:1924124/1‑118 (MQ=255)
                                                                   gCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCtgtg                                                                                                                <  1:1924124/118‑1 (MQ=255)
                                                                   gCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTc                                                                                >  1:598407/1‑150 (MQ=255)
                                                                     aTTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGc                                                                              >  1:412926/1‑150 (MQ=255)
                                                                     aTTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGc                                                                              <  2:2087751/150‑1 (MQ=255)
                                                                          tggCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGc                                                                         >  2:766360/1‑150 (MQ=255)
                                                                          tggCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGc                                                                         <  2:590788/150‑1 (MQ=255)
                                                                          tggCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGc                                                                         <  1:1009440/150‑1 (MQ=255)
                                                                              gggCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTg                                                                     >  2:1207617/1‑150 (MQ=255)
                                                                              gggCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTg                                                                     >  2:1249910/1‑150 (MQ=255)
                                                                              gggCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTg                                                                     >  1:1443068/1‑150 (MQ=255)
                                                                                gCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAt                                                                                                                                              >  1:1954450/1‑75 (MQ=255)
                                                                                gCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGAt                                                                                                                                              <  2:1954450/75‑1 (MQ=255)
                                                                                    aTTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTg                                                               >  1:336106/1‑150 (MQ=255)
                                                                                      tGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTc                                                             <  1:2345595/150‑1 (MQ=255)
                                                                                        cGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCgg                                                           >  2:803786/1‑150 (MQ=255)
                                                                                         gATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCt                                                                                                                   >  1:1421168/1‑93 (MQ=255)
                                                                                         gATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCt                                                                                                                   <  2:1421168/93‑1 (MQ=255)
                                                                                         gATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATAGCGTTTCGCCATGCGTTGCCGCTGTCGGc                                                          >  1:927377/1‑150 (MQ=255)
                                                                                           ttCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGt                                                        <  2:1430441/150‑1 (MQ=255)
                                                                                              cGGCGCTGGGCGTGCTCTTCTGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGGGGTGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCa                                                     >  2:426387/1‑150 (MQ=255)
                                                                                                gcgcTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGc                                                   <  2:217573/150‑1 (MQ=255)
                                                                                                   cTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGcttc                                                <  2:336106/150‑1 (MQ=255)
                                                                                                      ggCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTg                                             <  2:1821970/150‑1 (MQ=255)
                                                                                                        cGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCACAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGtt                                           <  2:2403301/150‑1 (MQ=255)
                                                                                                               ttCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCctgg                                    <  2:412926/150‑1 (MQ=255)
                                                                                                                  ggCCTGGACCAGCAGTTGGCGCTGGGGCCCGTGTTGGTGATGGTGGTGCATCACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCctggctg                                 <  2:1330437/150‑1 (MQ=255)
                                                                                                                  ggCCTGGACCAGCAGTTGGCGCAGGGGACCGTGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCAc                                                                                                      >  2:1090170/1‑81 (MQ=255)
                                                                                                                  ggCCTGGACCAGCAGTTGGCGCAGGGGACCGTGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCAc                                                                                                      <  1:1090289/81‑1 (MQ=255)
                                                                                                                  ggCCTGGACCAGCAGTTGGCGCAGGGGACCGTGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCAc                                                                                                      >  2:1090289/1‑81 (MQ=255)
                                                                                                                  ggCCTGGACCAGCAGTTGGCGCAGGGGACCGTGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCAc                                                                                                      <  1:1090170/81‑1 (MQ=255)
                                                                                                                   gCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCctggctgg                                >  1:1564062/1‑150 (MQ=255)
                                                                                                                    ccTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGGGGCGGTATCACCAGCGCAACCGCATAGCGTTTCGCCATGCTGTGCTGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTggg                               <  2:539946/150‑1 (MQ=255)
                                                                                                                       ggACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTc                            <  2:1237882/150‑1 (MQ=255)
                                                                                                                                gTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCAcc                                                                                                     >  2:1201057/1‑68 (MQ=255)
                                                                                                                                gTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCAcc                                                                                                     <  1:1201057/68‑1 (MQ=255)
                                                                                                                                 ttGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGt                                                                                   >  1:697744/1‑85 (MQ=255)
                                                                                                                                 ttGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGt                                                                                   <  2:697744/85‑1 (MQ=255)
                                                                                                                                     cgcAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTgg              <  1:1437375/150‑1 (MQ=255)
                                                                                                                                         ggggACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGt                                                                       >  2:1631834/1‑89 (MQ=255)
                                                                                                                                         ggggACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGt                                                                       <  1:1631834/89‑1 (MQ=255)
                                                                                                                                         ggggACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCGTGGCGGGGGTCGATCTTGGCGGTGGGGCt          >  2:1229919/1‑150 (MQ=255)
                                                                                                                                               cgcgTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTTGCTGGGGTCGATCTGGGCGGTGGGGCTTGATgc    <  1:803786/150‑1 (MQ=255)
                                                                                                                                               cgcgTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGATgc    >  2:2033929/1‑150 (MQ=255)
                                                                                                                                                gcgTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGATgcg   >  2:775926/1‑150 (MQ=255)
                                                                                                                                                 cgTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCAATACAGCCTCGCGGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGATgcgc  >  1:151154/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CTGATGATCGAGTGGTTGTTGTATATCGTGCTGGGCGCAGCCCTAGGCACCATGGGGGGCTTGTTTGGCATTGGTGGCGGGCTGATTGCGATTCCGGCGCTGGGCGTGCTCTTCGGCCTGGACCAGCAGTTGGCGCAGGGGACCGCGTTGGTGATGGTGGTGCCGAACGTGCTGCTGGCGCTGTGGCGGTATCACCAGCGCAACCGCATCGCGTTTCGCCATGCGTTGCCGCTGTCGGCGTGCAGCTTCTTGTTCGCCTGGCTGGGGTCGATCTGGGCGGTGGGGCTTGATGCGC  >  NC_002947/2083387‑2083681

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: