Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,408,096:1 +T 6.3% coding (491/615 nt) PP_3021 ← LysE family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,408,0961.T 6.3% 318.8 / 15.3 79coding (491/615 nt)PP_3021LysE family transporter
Reads supporting (aligned to +/- strand):  ref base . (34/40);  new base T (2/3);  total (36/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

GTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGC  >  NC_002947/3407947‑3408243
                                                                                                                                                      |                                                                                                                                                   
gTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCa                                                                                                                                                      >  1:1418022/1‑150 (MQ=255)
 tCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCAc                                                                                                                                                    <  2:1356377/150‑1 (MQ=255)
   gACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGACCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACgg                                                                                                                                                  >  2:426895/1‑150 (MQ=255)
       tGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGAc                                                                                                                                              <  1:686280/150‑1 (MQ=255)
         cgaaCCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCa                                                                                                                                            <  1:1243619/147‑1 (MQ=255)
              cTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAg                                                                                                                                       <  2:1248676/150‑1 (MQ=255)
                tGCAACAAGTCAGCCCCGGGCCCGCTGCAGCAACCACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTc                                                                                                                                     <  2:1320224/150‑1 (MQ=255)
                tGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTCGCCGTGGCCATCGGCAACCAGcgggta                                                                                                                                     <  2:678029/150‑6 (MQ=255)
                 gCAACAAGTCAGCCCCGGGACAGATGCAGCAACAACGCCACCACCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTcc                                                                                                                                    <  2:1186132/150‑1 (MQ=255)
                   aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCa                                                                                                                                  >  1:1155633/1‑150 (MQ=255)
                   aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCa                                                                                                                                  >  1:1454845/1‑150 (MQ=255)
                   aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCa                                                                                                                                  <  1:584833/150‑1 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCaa                                                                                                                                        <  1:7962/141‑1 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCaa                                                                                                                                        >  2:7962/1‑141 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCagag                                                                                                                               >  2:31445/1‑150 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCagag                                                                                                                               >  1:706265/1‑150 (MQ=255)
                           aGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCACGa                                                                                                                          <  1:1643097/150‑1 (MQ=255)
                             ccccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACg                                                                                                                        >  2:929493/1‑150 (MQ=255)
                              cccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGa                                                                                                                       <  1:426895/150‑1 (MQ=255)
                              cccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTCGCCGTGGCCATCGGCAACCAGcg                                                                                                                                         <  2:966898/132‑2 (MQ=255)
                              cccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTCGCCGTGGCCATCGGCAACCAGcg                                                                                                                                         >  1:966898/1‑131 (MQ=255)
                                     cAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGaa                                                                                                                <  1:214762/150‑1 (MQ=255)
                                             gcaACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATa                                                                                                        <  1:1635256/150‑1 (MQ=255)
                                             gcaACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGAAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATa                                                                                                        >  1:1469135/1‑150 (MQ=255)
                                               aacaacGCCCCCGCCGACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAgc                                                                                                      <  2:1146901/150‑1 (MQ=255)
                                               aacaacGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAgc                                                                                                      >  2:121801/1‑150 (MQ=255)
                                                  aacGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGcgcg                                                                                                   <  1:993490/150‑1 (MQ=255)
                                                   acGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGt                                                                                                  >  1:943509/1‑150 (MQ=255)
                                                       ccccGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGc                                                                                              >  1:1575626/1‑150 (MQ=255)
                                                             cAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtttt                                                                                        >  2:266101/1‑150 (MQ=255)
                                                              aaCATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                                       <  2:1155633/150‑1 (MQ=255)
                                                                   cACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGcc                                                                                  <  1:266101/150‑1 (MQ=255)
                                                                    aCCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCg                                                                                 >  2:689558/1‑150 (MQ=255)
                                                                      cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCggg                                                                               <  1:121801/150‑1 (MQ=255)
                                                                        gCATGACCCCCGAGGCCCGCTCCACCCACGCCAGGCTGCGGGCGAAGCGCCGCACTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGACCATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGc                                                                             >  2:385030/1‑150 (MQ=255)
                                                                            gACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGAc                                                                         <  2:706265/150‑1 (MQ=255)
                                                                            gACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGAAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGAc                                                                         <  2:1469135/150‑1 (MQ=255)
                                                                               ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                                      >  2:493267/1‑150 (MQ=255)
                                                                               ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                                      >  1:1063528/1‑150 (MQ=255)
                                                                               ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                                      >  2:974777/1‑150 (MQ=255)
                                                                                cccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                                          >  2:343789/1‑129 (MQ=255)
                                                                                cccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                                          <  1:343789/129‑1 (MQ=255)
                                                                                 ccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGTTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCGTAGCGCGGGGGCTGTTTTCCAGCCGGAGCTGACGCTgg                                                                    >  2:591136/1‑150 (MQ=255)
                                                                                 ccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTgg                                                                    >  2:271732/1‑150 (MQ=255)
                                                                                 ccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTgg                                                                    >  1:153938/1‑150 (MQ=255)
                                                                                  cGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGc                                                                   <  1:929493/150‑1 (MQ=255)
                                                                                    aGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGa                                                                 <  2:1454845/150‑1 (MQ=255)
                                                                                            tCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                                         >  1:555333/1‑150 (MQ=255)
                                                                                            tCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                                         <  2:555333/150‑1 (MQ=255)
                                                                                              cAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAAGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgg                                                                                <  1:1347169/125‑1 (MQ=255)
                                                                                              cAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAAGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgg                                                                                >  2:1347169/1‑125 (MQ=255)
                                                                                                  cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCTGCTGGTCGCCGTGGCCATCGGCAACCAGcgg                                                                                                                                        <  1:1530595/65‑3 (MQ=9)
                                                                                                  cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCTGCTGGTCGCCGTGGCCATCGGCAACCAGcgg                                                                                                                                        >  2:1530595/1‑63 (MQ=9)
                                                                                                  cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAgg                                                   <  2:153938/150‑1 (MQ=255)
                                                                                                   aCGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGc                                                  <  2:1227504/150‑1 (MQ=255)
                                                                                                       cAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                                       >  2:1358540/1‑109 (MQ=255)
                                                                                                       cAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                                       <  1:1358540/109‑1 (MQ=255)
                                                                                                            tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGc                                              <  1:394486/145‑1 (MQ=255)
                                                                                                            tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGc                                              >  2:394486/1‑145 (MQ=255)
                                                                                                            tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCCGGCATAGa                                         >  2:304601/1‑150 (MQ=255)
                                                                                                             gCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGaa                                        <  2:1063528/150‑1 (MQ=255)
                                                                                                                 gCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCTCTAGCGCGCAGGCCGTTCACCAGCCGGGGCGGGCGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAgc                                    <  1:31445/150‑1 (MQ=255)
                                                                                                                 gCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAgc                                    <  1:493267/150‑1 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTg                                                                     <  2:525599/113‑1 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTg                                                                     >  1:525599/1‑113 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                                         >  1:1076616/1‑125 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                                         <  2:1076616/125‑1 (MQ=255)
                                                                                                                         ccgcAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgggg                                                                              <  2:137887/100‑1 (MQ=255)
                                                                                                                         ccgcAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgggg                                                                              >  1:137887/1‑100 (MQ=255)
                                                                                                                          cgcAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGTCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTg                           <  1:385030/150‑1 (MQ=255)
                                                                                                                           gcAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTgg                          <  2:1575626/150‑1 (MQ=255)
                                                                                                                                 cccGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGa                    <  1:689558/150‑1 (MQ=255)
                                                                                                                                 cccGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGa                    <  1:1408005/150‑1 (MQ=255)
                                                                                                                                  ccGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAg                   >  1:126341/1‑150 (MQ=255)
                                                                                                                                   cGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGt                                                                                           >  1:1567760/1‑77 (MQ=255)
                                                                                                                                   cGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGt                                                                                           <  2:1567760/77‑1 (MQ=255)
                                                                                                                                             cGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTgg                                               <  2:702281/111‑1 (MQ=255)
                                                                                                                                             cGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTgg                                               >  1:702281/1‑111 (MQ=255)
                                                                                                                                                   ccACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGtt                      >  1:240350/1‑130 (MQ=255)
                                                                                                                                                   ccACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGtt                      <  2:240350/130‑1 (MQ=255)
                                                                                                                                                   ccACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGc  >  2:1580074/1‑150 (MQ=255)
                                                                                                                                                     aCGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACaag      >  1:1027219/1‑144 (MQ=255)
                                                                                                                                                     aCGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACaag      <  2:1027219/144‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCACGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAGGATGCCAGACAAGAAGC  >  NC_002947/3407947‑3408243

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: