Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,213,230 G→C 8.0% G19A (GGA→GCA) ‡ PP_1956 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,213,2300GC 8.0% 190.0 / 4.5 62G19A (GGA→GCA) ‡PP_1956hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (25/32);  new base C (2/3);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

TCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCCTACCTACTGGACAGAGGAC  >  NC_002947/2213084‑2213375
                                                                                                                                                  |                                                                                                                                                 
tctTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATTTGGTTTTTgg                                                                                                                                                   <  5:76337/147‑1 (MQ=255)
  tttGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAAt                                                                                                                                               <  4:238057/149‑1 (MQ=255)
    tGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGaaa                                                                                                                                                <  5:119120/146‑1 (MQ=255)
      cAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATAcc                                                                                                                                            <  5:330765/148‑1 (MQ=255)
         gTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACAACAt                                                                                                                                        >  6:279209/1‑149 (MQ=255)
           aaTCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACAtt                                                                                                                                       <  2:55242/148‑1 (MQ=255)
            aTCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTtg                                                                                                                                     <  4:238291/149‑1 (MQ=255)
             tCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgc                                                                                                                                      >  2:227193/1‑145 (MQ=255)
             tCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgc                                                                                                                                      <  1:227193/147‑3 (MQ=255)
                ggTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTtgctgc                                                                                                                                 >  6:139252/1‑149 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGcc                                                                                                                                >  3:261929/1‑149 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGcc                                                                                                                                >  4:128468/1‑149 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgcag                                                                                                                                    <  5:208999/145‑5 (MQ=255)
                 gTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTATTTCCAAAAACCACAtgcag                                                                                                                                    <  7:272388/145‑5 (MQ=255)
                     aaaCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGaaa                                                                                                                            <  5:139252/149‑1 (MQ=255)
                      aaCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGaaaa                                                                                                                           <  6:128929/149‑1 (MQ=255)
                       aCCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGaaaa                                                                                                                           >  6:328907/1‑148 (MQ=255)
                         cGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAatat                                                                                                                        <  7:43992/149‑1 (MQ=255)
                           gAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCa                                                                                                                      <  1:199518/149‑1 (MQ=255)
                            accAACCGCCTCCAAAAATAAAAAATTCCTGCCTCAAAAGGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCaa                                                                                                                     <  1:111099/147‑1 (MQ=255)
                                aCCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAagcca                                                                                                                 >  8:106317/1‑149 (MQ=255)
                                      tCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGc                                                                                                           >  6:207735/1‑149 (MQ=255)
                                        cAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGccc                                                                                                         <  4:261929/149‑1 (MQ=255)
                                         aataatAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGccc                                                                                                         >  1:215907/1‑148 (MQ=255)
                                              taaCAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTgg                                                                                                   >  7:93297/1‑149 (MQ=255)
                                                aCAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggcg                                                                                                 <  5:207735/149‑1 (MQ=255)
                                                 cAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAAtt                                                                                                     <  6:113089/144‑1 (MQ=255)
                                                 cAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAAtt                                                                                                     >  5:113089/1‑144 (MQ=255)
                                                  aaTATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggcgg                                                                                                >  1:245427/1‑148 (MQ=255)
                                                    tatGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCAc                                                                                             <  3:128468/149‑1 (MQ=255)
                                                           ccGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGata                                                                                        <  6:327092/147‑1 (MQ=255)
                                                           ccGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATAtt                                                                                      >  6:154838/1‑149 (MQ=255)
                                                            cGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTc                                                                                     <  6:115353/149‑1 (MQ=255)
                                                               cATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGAc                                                                                  <  5:301074/149‑1 (MQ=255)
                                                                aTTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACa                                                                                 <  5:154838/149‑1 (MQ=255)
                                                                     cgaagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGa                                                                            <  8:93297/149‑1 (MQ=255)
                                                                      gaagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTgg                                                                                                   <  7:262570/125‑1 (MQ=255)
                                                                      gaagcgaaTGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTgg                                                                                                   >  8:262570/1‑125 (MQ=255)
                                                                                 tGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCt                                                                        <  5:175207/141‑1 (MQ=255)
                                                                                 tGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCt                                                                        >  6:175207/1‑141 (MQ=255)
                                                                                 tGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTccac                                                                >  4:71176/1‑149 (MQ=255)
                                                                                        tttATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGTGGCACTGATATTCGATAAGCGATTCTACCTCCACCAGgtgt                                                         >  6:189654/1‑149 (MQ=255)
                                                                                               ataataAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATgg                                                   <  5:328907/148‑1 (MQ=255)
                                                                                                   taaGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTgg                                              >  4:278829/1‑149 (MQ=255)
                                                                                                         cgcacgcTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTg                                        <  2:215907/149‑1 (MQ=255)
                                                                                                               cTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAg                                  <  7:106317/149‑1 (MQ=255)
                                                                                                                tATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTgg                                              >  7:216600/1‑136 (MQ=255)
                                                                                                                tATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTgg                                              <  8:216600/136‑1 (MQ=255)
                                                                                                                    tcGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCtt                             >  6:309300/1‑149 (MQ=255)
                                                                                                                      gcgGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTg                            >  7:201926/1‑148 (MQ=255)
                                                                                                                            aTTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggc                                                                                                  <  8:146140/72‑1 (MQ=255)
                                                                                                                            aTTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTggc                                                                                                  >  7:146140/1‑72 (MQ=255)
                                                                                                                             ttGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGTGGCACTGATATTCGATAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTcct                    <  5:189654/149‑1 (MQ=255)
                                                                                                                                cagcaaaTGTGGTATTTCCAAAAACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTcctacc                 >  7:180202/7‑149 (MQ=255)
                                                                                                                                 ccTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTcctacc                 >  1:81674/1‑148 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                              <  6:327446/84‑1 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                              <  6:114772/84‑1 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                              >  5:327446/1‑84 (MQ=255)
                                                                                                                                    gCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGc                                                                              >  5:114772/1‑84 (MQ=255)
                                                                                                                                     cATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCctacctac              <  5:279209/147‑1 (MQ=255)
                                                                                                                                      aTGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGTCCCAATTGCCGGCACTGATCTTCGACAAGCGATTCTACCTCCACCAGTTGATTATGGCCTGGGCGTTGTTGTAGGGCTTGGATTTCCTACCTACTg            <  8:97069/148‑1 (MQ=255)
                                                                                                                                            tttttGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCCTACCTACTGGACagag     <  3:71176/149‑1 (MQ=255)
                                                                                                                                               ttGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGAGGAGACCTCCACCAGGAGTTATTTGAATGGGCGTTGGTGTAGGGCTTGGATTTACTACCTACTGGACAGCGGAc  >  1:271652/1‑149 (MQ=255)
                                                                                                                                                  gAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCtt                             <  2:172278/119‑1 (MQ=255)
                                                                                                                                                  gAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCtt                             >  1:172278/1‑119 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TCTTTGCAAGTAATCTGGTAAAAACCGGAGCAACCGGCTCCAATAATAACAATATGCTGCCGACATTCGCGAAGCGAATGGTGAGGTTTTTATGAATAATAAGTACGCACGCTATCTCGCGGTAATTGCCCTGCATGTGGTTTTTGGAAATACCACATTTGCTGCCGAAAATATCAAGCCAGCCGGCCCAATTGGCGGCACTGATATTCGACAAGCGATTCTACCTCCACCAGGTGTTTATGGCCTGGGCGTTGGTGTAGGGCTTGGATTTCCTACCTACTGGACAGAGGAC  >  NC_002947/2213084‑2213375

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