Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 966,599 A→T 5.3% L41M (TTG→ATG) ‡ PP_0828 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947966,5990AT 5.3% 292.7 / 5.1 93L41M (TTG→ATG) ‡PP_0828membrane protein
Reads supporting (aligned to +/- strand):  ref base A (45/43);  new base T (3/2);  total (48/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.95e-01

AACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATG  >  NC_002947/966458‑966747
                                                                                                                                             |                                                                                                                                                    
aaCTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAACATGTTCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCa                                                                                                                                              >  1:18409/1‑150 (MQ=255)
 aCTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCAt                                                                                                                                             >  2:1403942/1‑150 (MQ=255)
     ttCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaag                                                                                                                                            <  2:1857402/147‑3 (MQ=255)
     ttCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaag                                                                                                                                            >  1:1857402/1‑145 (MQ=255)
     ttCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGaaa                                                                                                                                         >  2:1927052/1‑150 (MQ=255)
      tCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAAt                                                                                                                                        <  2:1426828/150‑1 (MQ=255)
      tCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAAt                                                                                                                                        <  1:31978/150‑1 (MQ=255)
        aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCg                                                                                                                                      >  1:1441672/1‑150 (MQ=255)
         ccAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCCTGAAATCGa                                                                                                                                     >  1:404715/1‑150 (MQ=255)
         ccAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGa                                                                                                                                     <  1:494531/150‑1 (MQ=255)
         ccAGCGCGCCCCACAGGGTGCGAGCCTCCAACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGa                                                                                                                                     <  2:347175/150‑1 (MQ=255)
          cAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGAt                                                                                                                                    >  2:996606/1‑150 (MQ=255)
           aGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATg                                                                                                                                   <  1:1468687/150‑1 (MQ=255)
            gcgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGa                                                                                                                                     <  1:1566373/147‑1 (MQ=255)
            gcgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGa                                                                                                                                     >  2:1566373/1‑147 (MQ=255)
              gcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTg                                                                                                                                <  1:874443/150‑1 (MQ=255)
                  cccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATAGATGGTGACCa                                                                                                                            >  1:1174298/1‑150 (MQ=255)
                   ccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaa                                                                                                                                             >  1:1528685/1‑131 (MQ=255)
                   ccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCTTGAACAGCaa                                                                                                                                             <  2:1528685/132‑2 (MQ=255)
                   ccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAg                                                                                                                           >  1:1057336/1‑150 (MQ=255)
                     acaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGc                                                                                                                                               <  2:1816341/128‑1 (MQ=255)
                     acaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGc                                                                                                                                               >  1:1816341/1‑128 (MQ=255)
                             cGAGCCTCCCACCAGCGGCTGTACACGCTGATCCTACGAAAATCCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAg                                                                                                                 <  2:18409/150‑1 (MQ=255)
                                 ccTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATc                                                                                                                                       <  2:232189/124‑1 (MQ=255)
                                 ccTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATc                                                                                                                                       >  1:232189/1‑124 (MQ=255)
                                   tCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACCCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagc                                                                                                           <  1:1788769/150‑1 (MQ=255)
                                     ccaccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATca                                                                                                             <  2:383705/146‑1 (MQ=255)
                                     ccaccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATca                                                                                                             >  1:383705/1‑146 (MQ=255)
                                      caccaGCGGCTGTAGGCGCTGCTGTTGTGAAAACTCACCAGCACCACCAGCGCTGAGCTCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                        >  1:1323056/1‑150 (MQ=255)
                                           gCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGaa                                                                                                   >  1:200787/1‑150 (MQ=255)
                                               cTGTAGGCGCTGCTGTTGCGAACACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGGCAGCGGGATCAGTGGCAGGGTGAACTTGCTGCTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGc                                                                                               <  1:218991/150‑1 (MQ=255)
                                                tgtacccgcTTCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGAAGGGAGAACTTGCTGTTCCAGAACAGCATAAAATCGACGGTGACCAGCACATCCAAGATAAGCAGCCAGAACAACg                                                                                              >  2:862986/7‑150 (MQ=255)
                                                  tAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGAc                                                                                            <  2:1708046/150‑1 (MQ=255)
                                                  tAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGAc                                                                                            >  1:1644378/1‑150 (MQ=255)
                                                        ctgctgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCg                                                                                      <  2:924069/150‑1 (MQ=255)
                                                              ttGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCa                                                                                                                            >  2:1835021/1‑106 (MQ=255)
                                                              ttGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCa                                                                                                                            <  1:1835021/106‑1 (MQ=255)
                                                              ttGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCGGCCAGAACAGCGACCAGCCGATATAg                                                                                <  2:404715/150‑1 (MQ=255)
                                                              ttGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAg                                                                                >  2:1587225/1‑150 (MQ=255)
                                                               tGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGc                                                                               <  2:1644378/150‑1 (MQ=255)
                                                               tGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGc                                                                               <  2:683652/150‑1 (MQ=255)
                                                                 cGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCg                                                                             >  1:1196120/1‑150 (MQ=255)
                                                                 cGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCg                                                                             >  2:314239/1‑150 (MQ=255)
                                                                  gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            <  2:200787/150‑1 (MQ=255)
                                                                  gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            <  2:881548/150‑1 (MQ=255)
                                                                  gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                            <  2:29890/150‑1 (MQ=255)
                                                                   aaaaCTCACCAGCACCACCAGCGCTGAGCCCGGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGAt                                                                                                               >  1:280111/1‑114 (MQ=255)
                                                                   aaaaCTCACCAGCACCACCAGCGCTGAGCCCGGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGAt                                                                                                               <  2:280111/114‑1 (MQ=255)
                                                                       cTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTc                                                                       <  2:1067456/150‑1 (MQ=255)
                                                                        tCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGCGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCt                                                                      >  2:1206718/1‑150 (MQ=255)
                                                                          aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTc                                                                    >  2:730254/1‑150 (MQ=255)
                                                                          aCCAGAACCACGAGCGATGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTc                                                                    <  2:1807196/150‑1 (MQ=255)
                                                                            cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcac                                                                  >  1:129029/1‑150 (MQ=255)
                                                                            cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcac                                                                  <  1:1105763/150‑1 (MQ=255)
                                                                            cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCCATGGTCTTcac                                                                  >  2:1315891/1‑150 (MQ=255)
                                                                               caccacCAGCGCTGAGCCCATCAAGGTCAGCCGCATCAGTGCCAGGGTTAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCATCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcaccac                                                               <  1:739084/150‑1 (MQ=255)
                                                                                  caccaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGc                                                                                               >  1:1765071/1‑115 (MQ=255)
                                                                                  caccaGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGc                                                                                               <  2:1765071/115‑1 (MQ=255)
                                                                                        cgcTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAGCTTGGTGTTTAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTa                                                      <  2:1323056/150‑1 (MQ=255)
                                                                                         gcTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTAc                                                     <  1:314239/150‑1 (MQ=255)
                                                                                              gCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCg                                                >  1:677907/1‑150 (MQ=255)
                                                                                                ccAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACtt                                                   >  2:592719/1‑145 (MQ=255)
                                                                                                ccAGCAGGGTCAGAGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACtt                                                   <  1:592719/145‑1 (MQ=255)
                                                                                                 cagcagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGc                                             >  2:855030/1‑150 (MQ=255)
                                                                                                        gTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGAc                                       <  1:1590295/149‑1 (MQ=255)
                                                                                                        gTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGAc                                       >  2:1590295/1‑149 (MQ=255)
                                                                                                          cAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          <  2:1307167/144‑1 (MQ=255)
                                                                                                          cAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTc                                    <  1:730254/150‑1 (MQ=255)
                                                                                                          cAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTc                                    >  2:239480/1‑150 (MQ=255)
                                                                                                          cAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAAAAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGAt                                          >  1:1307167/1‑144 (MQ=255)
                                                                                                           aGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCt                                   <  2:129029/150‑1 (MQ=255)
                                                                                                            gCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTc                                                                    >  1:464314/1‑116 (MQ=255)
                                                                                                            gCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTc                                                                    <  2:464314/116‑1 (MQ=255)
                                                                                                                    aGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGa                          <  1:996606/150‑1 (MQ=255)
                                                                                                                     gTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                                   <  2:791393/108‑1 (MQ=255)
                                                                                                                     gTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                                   >  1:791393/1‑108 (MQ=255)
                                                                                                                      tGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACAt                                                                                                                      >  1:681939/1‑56 (MQ=255)
                                                                                                                      tGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACAt                                                                                                                      <  2:681939/56‑1 (MQ=255)
                                                                                                                         cAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTa                     >  1:1064477/1‑150 (MQ=255)
                                                                                                                            ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGc                                                          >  1:1058290/1‑110 (MQ=255)
                                                                                                                            ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGc                                                          <  2:1058290/110‑1 (MQ=255)
                                                                                                                            ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTc                                >  1:1150219/1‑136 (MQ=255)
                                                                                                                            ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTc                                <  2:1150219/136‑1 (MQ=255)
                                                                                                                            ggTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCt                  >  1:405402/1‑150 (MQ=255)
                                                                                                                             gTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTc                 <  1:1927052/150‑1 (MQ=255)
                                                                                                                              tGAACTTGCTGTTCAAGAAGAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCg                >  1:1425626/1‑150 (MQ=255)
                                                                                                                               tttcatGCTGTTCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGa               >  2:1751070/6‑150 (MQ=255)
                                                                                                                                  cTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAAc            >  2:165470/1‑150 (MQ=255)
                                                                                                                                       tGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggg       >  1:979021/1‑150 (MQ=255)
                                                                                                                                         ttCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggggg     >  2:253271/1‑150 (MQ=255)
                                                                                                                                         ttCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggggg     >  2:593983/1‑150 (MQ=255)
                                                                                                                                         ttCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggg       <  1:1903862/148‑1 (MQ=255)
                                                                                                                                         ttCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCggg       >  2:1903862/1‑148 (MQ=255)
                                                                                                                                            aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATg  <  1:253271/150‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
AACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTTCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGGGATG  >  NC_002947/966458‑966747

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: