Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,560,509 T→A 7.9% L143Q (CTG→CAG)  PP_2243 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,560,5090TA 7.9% 190.8 / 7.0 63L143Q (CTG→CAG) 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 = 5.85e-01

CGACCTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGCAATTGAGCTTT  >  NC_002947/2560363‑2560654
                                                                                                                                                  |                                                                                                                                                  
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                                                                                                                                                >  1:237292/1‑105 (MQ=255)
                                              gAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTggcg                                                                                                                                                <  2:237292/105‑1 (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                                                                                >  2:911955/1‑148 (MQ=255)
                                                                  cGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                <  1:911955/148‑1 (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                                                                          >  1:947249/1‑103 (MQ=255)
                                                                                                                     tGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAg                                                                          <  2:947249/103‑1 (MQ=255)
                                                                                                                     tGCGCTATGCCAGCGCCAGACCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc                           >  2:814943/1‑150 (MQ=255)
                                                                                                                        gcTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc                           >  1:1419142/1‑147 (MQ=255)
                                                                                                                        gcTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGcc                           <  2:1419142/147‑1 (MQ=255)
                                                                                                                         cTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGCCGGCCCGGCCACTTCTTTGGCGGTGACCCACCCGAGGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTg                       >  1:603037/1‑150 (MQ=255)
                                                                                                                             gCCAGCGCCAGGCCGCCTGGCTGGCGGCCTCCGACAGGCCGGGTGACTGTCTGGCTGTGACGCACCCGCTCGTCATGCGCGCCATGCTGGTACAGAGGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCTACGTGCTGCCGCTGTCAc                   <  1:813128/150‑1 (MQ=255)
                                                                                                                              ccAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACg                  >  2:23935/1‑150 (MQ=255)
                                                                                                                               cAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                >  2:589359/1‑87 (MQ=255)
                                                                                                                               cAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTg                                                                                <  1:589359/87‑1 (MQ=255)
                                                                                                                                   gcCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGc             >  2:1052554/1‑150 (MQ=255)
                                                                                                                                   cgaaGGCCGCCAGCCAGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAg                                                      <  2:718554/106‑1 (MQ=255)
                                                                                                                                   cgaaGGCCGCCAGCCAGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAg                                                      >  1:718554/4‑109 (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)
                                                                                                                                                  |                                                                                                                                                  
CGACCTTGGGCGTTGGCAAGGTGTGTCGTTGAAGCAACTACAAGCCGAGCAACCGCAGGCGTTGGCCGAATGGCTGCAGGATCCGGCCAGTGCCGTGCATGGTGGTGAGTCGTTCGCTGCGCTATGCCAGCGCCAGGCCGCCTGGCTGGCGGCCTTCGACAGGCCGGGTGACTGGCTGGCGGTGACCCACCCGATGGTCATGCGCGCCATGCTGGTACAGGTGCTTGGGTGTCCGATGAGTGCCGGGCAGCGTATCGACGTGCTGCCGCTGTCACGTCTGCAATTGAGCTTT  >  NC_002947/2560363‑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: