Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,560,505 T→A 7.9% W142R (TGG→AGG) ‡ PP_2243 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,560,5050TA 7.9% 196.0 / 7.0 63W142R (TGG→AGG) ‡PP_2243hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/26);  new base A (2/3);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.30e-01

CCGATTGCGACCTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGCAATTGAGCTTT  >  NC_002947/2560356‑2560654
                                                                                                                                                     |                                                                                                                                                      
ccGATTGCGACCTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCct                                                                                                                                                        <  2:507484/150‑1 (MQ=255)
       cGACCTTCTGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCtggctggc                                                                                                                                                 <  1:61645/150‑1 (MQ=255)
        gACCTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTggcg                                                                                                                                                >  1:914894/1‑150 (MQ=255)
        gACCTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTggcg                                                                                                                                                >  1:549354/1‑150 (MQ=255)
           cTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGcc                                                                                                                                             >  1:933641/1‑150 (MQ=255)
                gCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGa                                                                                                                                        <  1:715676/150‑1 (MQ=255)
                 cGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGAc                                                                                                                                       <  2:914894/150‑1 (MQ=255)
                 cGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCAAGGCCGCCAGCCAGGCGGCCtggcgc                                                                                                                                       <  1:440693/150‑6 (MQ=255)
                  gTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGGATGGATGCAGGATCCGGCCAGTGCCGTGCAGGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCGTTCTACa                                                                                                                                      >  2:813128/1‑150 (MQ=255)
                      gCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACCGGcc                                                                                                                                  >  2:421655/1‑150 (MQ=255)
                      gCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGcc                                                                                                                                  >  2:463888/1‑150 (MQ=255)
                            tgtgtCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGa                                                                                                                            <  2:810667/150‑1 (MQ=255)
                                 cGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGActggc                                                                                                                      <  1:463888/150‑1 (MQ=255)
                                        gCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGCGTCGTTGGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGActggctggggtt                                                                                                               >  1:41613/1‑146 (MQ=255)
                                           aCTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGac                                                                                                            <  2:516057/150‑1 (MQ=255)
                                               cAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGacccac                                                                                                        <  2:1281985/150‑1 (MQ=255)
                                                  gCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCg                                                                                                     >  2:992288/1‑150 (MQ=255)
                                                     gAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTggcg                                                                                                                                                <  2:237292/105‑1 (MQ=255)
                                                     gAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTggcg                                                                                                                                                >  1:237292/1‑105 (MQ=255)
                                                     gAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGa                                                                                                    <  1:499443/148‑1 (MQ=255)
                                                     gAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGa                                                                                                    >  2:499443/1‑148 (MQ=255)
                                                               aGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATgcgcg                                                                                        >  1:826396/1‑150 (MQ=255)
                                                                  cGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCa                                                                                     >  2:63036/1‑150 (MQ=255)
                                                                     tGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGc                                                                                  <  2:933641/150‑1 (MQ=255)
                                                                       gCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                >  1:1399523/1‑150 (MQ=255)
                                                                        ccGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCAGCCAGGCGGCCt                                                                                                                                            <  2:434507/90‑1 (MQ=255)
                                                                        ccGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCAGCCAGGCGGCCt                                                                                                                                            >  1:434507/1‑90 (MQ=255)
                                                                         cGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                <  1:911955/148‑1 (MQ=255)
                                                                         cGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                >  2:911955/1‑148 (MQ=255)
                                                                         cGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGt                                                                              >  2:316647/1‑150 (MQ=255)
                                                                          gAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTa                                                                             <  1:366189/150‑1 (MQ=255)
                                                                             tGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTc                                                                                               <  2:472117/129‑1 (MQ=255)
                                                                             tGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTc                                                                                               >  1:472117/1‑129 (MQ=255)
                                                                             tGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAg                                                                          <  2:826396/150‑1 (MQ=255)
                                                                                  gCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCAt                                                                                             >  1:353851/1‑126 (MQ=255)
                                                                                  gCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCAt                                                                                             <  2:353851/126‑1 (MQ=255)
                                                                                    aGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTg                                                                   <  2:909338/150‑1 (MQ=255)
                                                                                     ggATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTgg                                                                  <  2:549354/150‑1 (MQ=255)
                                                                                     ggATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTgg                                                                  <  1:63036/150‑1 (MQ=255)
                                                                                        tCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCCTGCTGGAGCAGGTGcccggcgg                                                               >  2:659448/1‑143 (MQ=255)
                                                                                          cGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTc                                                             >  2:702373/1‑150 (MQ=255)
                                                                                            gCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCg                                                           >  2:1297417/1‑150 (MQ=255)
                                                                                              cAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTgg                                                                               <  1:116024/128‑1 (MQ=255)
                                                                                              cAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTgg                                                                               >  2:116024/1‑128 (MQ=255)
                                                                                              cAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGAt                                                         <  2:1027912/150‑1 (MQ=255)
                                                                                                 tGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAg                                                      >  2:239006/1‑150 (MQ=255)
                                                                                                       gCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCgg                                                <  1:1297417/150‑1 (MQ=255)
                                                                                                       gCATGGTGGTGAGTCGTTCGCTGCGATATGCCAGCGCCAGGCCGCCTGGCTGGCGCCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCgg                                                <  2:41613/150‑1 (MQ=255)
                                                                                                                          gCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGtgctg                             <  2:1389733/150‑1 (MQ=255)
                                                                                                                            tgagatATGCCAGTGCCAGGCCGCTTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc                           >  2:1380258/6‑150 (MQ=255)
                                                                                                                            tGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAg                                                                          <  2:947249/103‑1 (MQ=255)
                                                                                                                            tGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAg                                                                          >  1:947249/1‑103 (MQ=255)
                                                                                                                            tGCGCTATGCCAGCGCCAGACCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc                           >  2:814943/1‑150 (MQ=255)
                                                                                                                               gcTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc                           <  2:1419142/147‑1 (MQ=255)
                                                                                                                               gcTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc                           >  1:1419142/1‑147 (MQ=255)
                                                                                                                                cTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGCCGGCCCGGCCACTTCTTTGGCGGTGACCCACCCGAGGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTg                       >  1:603037/1‑150 (MQ=255)
                                                                                                                                    gCCAGCGCCAGGCCGCCTGGCTGGCGGCCTCCGACAGGCCGGGTGACTGTCTGGCTGTGACGCACCCGCTCGTCATGCGCGCCATGCTGGTACAGAGGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCTACGTGCTGCCGCTGTCAc                   <  1:813128/150‑1 (MQ=255)
                                                                                                                                     ccAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACg                  >  2:23935/1‑150 (MQ=255)
                                                                                                                                      cAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                <  1:589359/87‑1 (MQ=255)
                                                                                                                                      cAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                >  2:589359/1‑87 (MQ=255)
                                                                                                                                          gcCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGc             >  2:1052554/1‑150 (MQ=255)
                                                                                                                                          cgaaGGCCGCCAGCCAGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAg                                                      >  1:718554/4‑109 (MQ=255)
                                                                                                                                          cgaaGGCCGCCAGCCAGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAg                                                      <  2:718554/106‑1 (MQ=255)
                                                                                                                                           ccAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCGTGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGCa            >  1:502692/1‑150 (MQ=255)
                                                                                                                                             aGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGCAAt          >  1:126200/1‑150 (MQ=255)
                                                                                                                                                     tggctggcGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAg                                                                          >  2:1243129/1‑78 (MQ=255)
                                                                                                                                                     tggctggcGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAg                                                                          <  1:1243129/78‑1 (MQ=255)
                                                                                                                                                     tggctggcGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGCAATTGAGCttt  >  2:752242/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCGATTGCGACCTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGCAATTGAGCTTT  >  NC_002947/2560356‑2560654

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: