Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,408,085 T→C 6.6% N168D (AAC→GAC)  PP_3021 ← LysE family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,408,0850TC 6.6% 216.7 / 4.7 76N168D (AAC→GAC) PP_3021LysE family transporter
Reads supporting (aligned to +/- strand):  ref base T (32/39);  new base C (2/3);  total (34/42)
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.75e-01

ATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAG  >  NC_002947/3407936‑3408226
                                                                                                                                                     |                                                                                                                                              
aTTTCACCTTGGTCGACATGCACACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGtt                                                                                                                                                >  1:1248676/1‑150 (MQ=255)
           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:1155633/1‑150 (MQ=255)
                              aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCa                                                                                                                 <  1:584833/150‑1 (MQ=255)
                              aacaaGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCa                                                                                                                 >  1:1454845/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                                                     >  2:493267/1‑150 (MQ=255)
                                                                                          ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                     >  1:1063528/1‑150 (MQ=255)
                                                                                          ccccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCt                                                     >  2:974777/1‑150 (MQ=255)
                                                                                           cccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                         <  1:343789/129‑1 (MQ=255)
                                                                                           cccGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGtt                                                                         >  2:343789/1‑129 (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                                        <  2:555333/150‑1 (MQ=255)
                                                                                                       tCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                        >  1:555333/1‑150 (MQ=255)
                                                                                                         cAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAAGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgg                                                               >  2:1347169/1‑125 (MQ=255)
                                                                                                         cAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAAGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCgg                                                               <  1:1347169/125‑1 (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                             <  1:394486/145‑1 (MQ=255)
                                                                                                                       tGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGc                             >  2:394486/1‑145 (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                                                    <  2:525599/113‑1 (MQ=255)
                                                                                                                                aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTg                                                    >  1:525599/1‑113 (MQ=255)
                                                                                                                                aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                        <  2:1076616/125‑1 (MQ=255)
                                                                                                                                aGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTg                                        >  1:1076616/1‑125 (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:1408005/150‑1 (MQ=255)
                                                                                                                                            cccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGa   <  1:689558/150‑1 (MQ=255)
                                                                                                                                             ccGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAg  >  1:126341/1‑150 (MQ=255)
                                                                                                                                              cGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGt                                                                          >  1:1567760/1‑77 (MQ=255)
                                                                                                                                              cGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGt                                                                          <  2:1567760/77‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
ATTTCACCTTGGTCGACATGCGCACCTTGCAACAAGTCAGCCCCGGGCCAGATGCAGCAACAACGCCCCCGCCAACATCACCAGCATGACCCCCGAGGCCCGCTCCAGCCACGCCAGGCTGCGGGCGAAGCGCCGCAGTACCCGCTGGTTGCCGATGGCCA‑CGGCGACCAGCAAGTCCCAGAGCAGGACGATGCTGAACATCCATAGCGCGTAGGCCGTTTTCCAGCCGGGGCTGACGCTGGCGACCATGCTGGCCAGGCTGGCATAGAACAGCGCGTTCTTGGGGTTGAG  >  NC_002947/3407936‑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: ‑