Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,408,090 Δ1 bp 6.4% coding (497/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,0900A. 6.4% 237.3 / 10.7 78coding (497/615 nt)PP_3021LysE family transporter
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base . (2/3);  total (35/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.05e-01

GTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAG  >  NC_002947/3407947‑3408226
                                                                                                                                               |                                                                                                                                         
gTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCa                                                                                                                                     >  1:1418022/1‑150 (MQ=255)
 tCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑c                                                                                                                                   <  2:1356377/150‑1 (MQ=255)
   gACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGACCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑Cgg                                                                                                                                 >  2:426895/1‑150 (MQ=255)
       tGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGAc                                                                                                                             <  1:686280/150‑1 (MQ=255)
         cgaaCCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCa                                                                                                                           <  1:1243619/147‑1 (MQ=255)
              cTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAg                                                                                                                      <  2:1248676/150‑1 (MQ=255)
                tGCAACAAGTCAGCCCCGGGCCCGCTGCAGCAACCACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTc                                                                                                                    <  2:1320224/150‑1 (MQ=255)
                tGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTCGCCG‑TGGCCATCGGCAACCAGcgggta                                                                                                                    <  2:678029/150‑6 (MQ=255)
                 gCAACAAGTCAGCCCCGGGACAGATGCAGCAACAACGCCACCACCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTcc                                                                                                                   <  2:1186132/150‑1 (MQ=255)
                   aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCa                                                                                                                 <  1:584833/150‑1 (MQ=255)
                   aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCa                                                                                                                 >  1:1454845/1‑150 (MQ=255)
                   aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCa                                                                                                                 >  1:1155633/1‑150 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCaa                                                                                                                       >  2:7962/1‑141 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCaa                                                                                                                       <  1:7962/141‑1 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCagag                                                                                                              >  1:706265/1‑150 (MQ=255)
                      aaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCagag                                                                                                              >  2:31445/1‑150 (MQ=255)
                           aGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCACGa                                                                                                         <  1:1643097/150‑1 (MQ=255)
                             ccccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACg                                                                                                       >  2:929493/1‑150 (MQ=255)
                              cccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGa                                                                                                      <  1:426895/150‑1 (MQ=255)
                              cccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTCGCCG‑TGGCCATCGGCAACCAGcg                                                                                                                        >  1:966898/1‑131 (MQ=255)
                              cccGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTCGCCG‑TGGCCATCGGCAACCAGcg                                                                                                                        <  2:966898/132‑2 (MQ=255)
                                     cAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGaa                                                                                               <  1:214762/150‑1 (MQ=255)
                                             gcaACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATa                                                                                       <  1:1635256/150‑1 (MQ=255)
                                             gcaACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGAAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATa                                                                                       >  1:1469135/1‑150 (MQ=255)
                                               aacaacGCCCCCGCCGACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAgc                                                                                     <  2:1146901/150‑1 (MQ=255)
                                               aacaacGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAgc                                                                                     >  2:121801/1‑150 (MQ=255)
                                                  aacGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGcgcg                                                                                  <  1:993490/150‑1 (MQ=255)
                                                   acGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGt                                                                                 >  1:943509/1‑150 (MQ=255)
                                                       ccccGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGc                                                                             >  1:1575626/1‑150 (MQ=255)
                                                             cAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtttt                                                                       >  2:266101/1‑150 (MQ=255)
                                                              aaCATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                      <  2:1155633/150‑1 (MQ=255)
                                                                   cACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGcc                                                                 <  1:266101/150‑1 (MQ=255)
                                                                    aCCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCg                                                                >  2:689558/1‑150 (MQ=255)
                                                                      cAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCggg                                                              <  1:121801/150‑1 (MQ=255)
                                                                        gCATGACCCCCGAGGCCCGCTCCACCCACGCCAGGCTGCGGGCGAAGCGCCGCACTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGACCATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGc                                                            >  2:385030/1‑150 (MQ=255)
                                                                            gACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGAc                                                        <  2:706265/150‑1 (MQ=255)
                                                                            gACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGAAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGAc                                                        <  2:1469135/150‑1 (MQ=255)
                                                                               ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                     >  1:1063528/1‑150 (MQ=255)
                                                                               ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                     >  2:974777/1‑150 (MQ=255)
                                                                               ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                     >  2:493267/1‑150 (MQ=255)
                                                                                cccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                         >  2:343789/1‑129 (MQ=255)
                                                                                cccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                         <  1:343789/129‑1 (MQ=255)
                                                                                 ccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGTTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCGTAGCGCGGGGGCTGTTTTCCAGCCGGAGCTGACGCTgg                                                   >  2:591136/1‑150 (MQ=255)
                                                                                 ccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTgg                                                   >  2:271732/1‑150 (MQ=255)
                                                                                 ccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTgg                                                   >  1:153938/1‑150 (MQ=255)
                                                                                  cGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGc                                                  <  1:929493/150‑1 (MQ=255)
                                                                                    aGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGa                                                <  2:1454845/150‑1 (MQ=255)
                                                                                            tCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                        >  1:555333/1‑150 (MQ=255)
                                                                                            tCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                        <  2:555333/150‑1 (MQ=255)
                                                                                              cAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAAGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgg                                                               <  1:1347169/125‑1 (MQ=255)
                                                                                              cAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAAGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgg                                                               >  2:1347169/1‑125 (MQ=255)
                                                                                                  cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCTGCTGGTCGCCG‑TGGCCATCGGCAACCAGcgg                                                                                                                       <  1:1530595/65‑3 (MQ=9)
                                                                                                  cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCTGCTGGTCGCCG‑TGGCCATCGGCAACCAGcgg                                                                                                                       >  2:1530595/1‑63 (MQ=9)
                                                                                                  cACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAgg                                  <  2:153938/150‑1 (MQ=255)
                                                                                                   aCGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGc                                 <  2:1227504/150‑1 (MQ=255)
                                                                                                       cAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                      >  2:1358540/1‑109 (MQ=255)
                                                                                                       cAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTc                                                                      <  1:1358540/109‑1 (MQ=255)
                                                                                                            tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGc                             >  2:394486/1‑145 (MQ=255)
                                                                                                            tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGc                             <  1:394486/145‑1 (MQ=255)
                                                                                                            tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCCGGCATAGa                        >  2:304601/1‑150 (MQ=255)
                                                                                                             gCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGaa                       <  2:1063528/150‑1 (MQ=255)
                                                                                                                 gCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCTCTAGCGCGCAGGCCGTTCACCAGCCGGGGCGGGCGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAgc                   <  1:31445/150‑1 (MQ=255)
                                                                                                                 gCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAgc                   <  1:493267/150‑1 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTg                                                    >  1:525599/1‑113 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTg                                                    <  2:525599/113‑1 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                        >  1:1076616/1‑125 (MQ=255)
                                                                                                                     aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                        <  2:1076616/125‑1 (MQ=255)
                                                                                                                         ccgcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgggg                                                             <  2:137887/100‑1 (MQ=255)
                                                                                                                         ccgcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgggg                                                             >  1:137887/1‑100 (MQ=255)
                                                                                                                          cgcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGTCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTg          <  1:385030/150‑1 (MQ=255)
                                                                                                                           gcAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTgg         <  2:1575626/150‑1 (MQ=255)
                                                                                                                                 cccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGa   <  1:689558/150‑1 (MQ=255)
                                                                                                                                 cccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGa   <  1:1408005/150‑1 (MQ=255)
                                                                                                                                  ccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAg  >  1:126341/1‑150 (MQ=255)
                                                                                                                                   cGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGt                                                                          <  2:1567760/77‑1 (MQ=255)
                                                                                                                                   cGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGt                                                                          >  1:1567760/1‑77 (MQ=255)
                                                                                                                                             cGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTgg                              >  1:702281/1‑111 (MQ=255)
                                                                                                                                             cGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTgg                              <  2:702281/111‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                         
GTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAG  >  NC_002947/3407947‑3408226

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: ‑