Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 754,763 C→A 8.3% intergenic (‑19/+149) PP_0643 ← / ← PP_0645 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947754,7630CA 8.3% 186.0 / 6.8 60intergenic (‑19/+149)PP_0643/PP_0645hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base A (2/3);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

GTGCAGCGGCAGCAGCCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAGCTCTGGCG  >  NC_002947/754616‑754906
                                                                                                                                                   |                                                                                                                                                
gTGCAGCGGCAGCAGCCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAattca                                                                                                                                                 <  2:106644/149‑1 (MQ=255)
         cagcagCCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATGAATTCGCCgcc                                                                                                                                        <  3:252967/149‑3 (MQ=255)
          agcagcCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGc                                                                                                                                       >  7:234120/1‑149 (MQ=255)
          agcagcCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATGAATTCGCCgccc                                                                                                                                       <  7:62484/149‑4 (MQ=255)
             agcCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCttt                                                                                                                                    <  4:107757/149‑1 (MQ=255)
             agcCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCttt                                                                                                                                    >  3:25010/1‑149 (MQ=255)
             agcCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCttt                                                                                                                                    <  4:30038/149‑1 (MQ=255)
                    gcgcgGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGtt                                                                                                                             <  8:68651/149‑1 (MQ=255)
                    gcgcgGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGtt                                                                                                                             >  7:68651/1‑149 (MQ=255)
                    gcgcgGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGtt                                                                                                                             >  6:191128/1‑149 (MQ=255)
                    gcgcgGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCACCGGGCTTTGCCTGtt                                                                                                                             <  6:74177/149‑1 (MQ=255)
                        ggTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGccc                                                                                                                         <  2:57391/149‑1 (MQ=255)
                             cgtccgtcACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCga                                                                                                                     <  3:32063/148‑1 (MQ=255)
                                   tcACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCcg                                                                                                              >  5:154137/1‑149 (MQ=255)
                                   tcACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCcg                                                                                                              >  7:5756/1‑149 (MQ=255)
                                     aCCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGTCCCCGATGATCcgcg                                                                                                            >  2:232207/1‑149 (MQ=255)
                                      ccTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGCGCTTTGCCTGTTGCCCCCGATGATCCgcgc                                                                                                           <  1:30038/149‑1 (MQ=255)
                                       ctggctggATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCa                                                                                                          >  5:330495/1‑149 (MQ=255)
                                       ctggctggATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCa                                                                                                          >  3:199959/1‑149 (MQ=255)
                                       ctggctggATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCa                                                                                                          >  4:118003/1‑149 (MQ=255)
                                            tggATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGtt                                                                                                     >  8:64476/1‑149 (MQ=255)
                                               atttCTTGGTGGCCAGCAAGTTCGGCACGCCCGCGTTGCTCAGGTGCCGCTCCTGCGCGCGGACGGTCGAGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGGTCATGGTTGGc                                                                                                  <  8:108178/146‑1 (MQ=255)
                                                tatCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCga                                                                                                                     >  1:189223/1‑129 (MQ=255)
                                                tatCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCga                                                                                                                     <  2:189223/129‑1 (MQ=255)
                                                tatCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGcc                                                                                                 >  5:74981/1‑149 (MQ=255)
                                                   cTTGGTGGCTAGCGACTTGGGCACGCCCTCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCa                                                                                                          <  8:132868/137‑1 (MQ=255)
                                                   cTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCAGCGCa                                                                                                          >  7:132868/1‑137 (MQ=255)
                                                      ggtggCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCgg                                                                                               <  8:5756/145‑1 (MQ=255)
                                                      ggtggCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCgg                                                                                               <  6:154137/145‑1 (MQ=255)
                                                          gCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAgg                                                                                        <  2:24452/148‑1 (MQ=255)
                                                             aGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGccc                                                                                     >  1:62346/1‑148 (MQ=255)
                                                               cGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAAc                                                                                  <  5:191128/149‑1 (MQ=255)
                                                                gACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACa                                                                                 <  6:197421/149‑1 (MQ=255)
                                                                   ttGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATc                                                                               >  4:45046/1‑148 (MQ=255)
                                                                     gggCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATc                                                                               <  3:45046/146‑1 (MQ=255)
                                                                     gggCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGa                                                                             <  8:131011/148‑1 (MQ=255)
                                                                      ggCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAg                                                                            <  3:118003/148‑1 (MQ=255)
                                                                           gCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCTTTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAg                                                                      <  3:73016/149‑1 (MQ=255)
                                                                           gCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAg                                                                      >  3:201791/1‑149 (MQ=255)
                                                                           gCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAg                                                                      >  3:94097/1‑149 (MQ=255)
                                                                                            gCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCa                                                      >  8:33116/1‑148 (MQ=255)
                                                                                              cGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAg                                                   >  6:24436/1‑149 (MQ=255)
                                                                                               gCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGc                                                  >  5:332774/1‑149 (MQ=255)
                                                                                                   cAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCt                                                <  4:25010/147‑1 (MQ=255)
                                                                                                    aGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑c                                              <  4:275478/147‑1 (MQ=255)
                                                                                                    aGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑c                                              >  3:275478/1‑147 (MQ=255)
                                                                                                        cgcgGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATtggt                                        <  7:64476/149‑1 (MQ=255)
                                                                                                                   gcgAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCggccagg                              <  6:74981/148‑1 (MQ=255)
                                                                                                                           aCGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGggccaggcca                     >  8:81508/1‑149 (MQ=255)
                                                                                                                           aCGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGggccaggcc                      <  2:62346/148‑1 (MQ=255)
                                                                                                                               aTTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCggccagggcca                          >  4:314276/1‑140 (MQ=255)
                                                                                                                               aTTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCggccagggcca                          <  3:314276/140‑1 (MQ=255)
                                                                                                                                     cagGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCtgtg             >  1:223407/1‑147 (MQ=255)
                                                                                                                                       gggCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGTCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCTACATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAg          <  1:232207/149‑1 (MQ=255)
                                                                                                                                        ggcggcGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAGct        <  6:208880/149‑1 (MQ=255)
                                                                                                                                        ggcggcGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAGct        >  3:282570/1‑149 (MQ=255)
                                                                                                                                        cccggcgAATGAATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCggccagggccag                         >  2:91145/3‑132 (MQ=255)
                                                                                                                                        cccggcgAATGAATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCggccagggccag                         <  1:91145/130‑1 (MQ=255)
                                                                                                                                         ccggcgAATGAATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAGctc       >  2:36770/2‑149 (MQ=255)
                                                                                                                                            gcgAATTCCTTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAGCTCTgg    <  4:201791/149‑1 (MQ=255)
                                                                                                                                              gAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAGCTCTGgcg  >  7:30414/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTGCAGCGGCAGCAGCCGGCGCGCGGTTTCGTCCGTCACCTGGCTGGATATCTTGGTGGCCAGCGACTTGGGCACGCCCGCGTTGCGCAGGTGCCGCTCCAGCGCGCGGACGGTCGCGAAATCACGCATTCAGCAGGGCGGCGAATTCATTCGCCGGGCTTTGCCTGTTGCCCCCGATGATCCGCGCATGGTTGGCCGGCTCCAGGCCCAACATCGAGCCGAAGCTGATCATCTGGCGCATAGCCT‑CATTGGTGGCGGCCAGGGCCAGGCCACCCCTGTGAAGCTCTGGCG  >  NC_002947/754616‑754906

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑