Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,605,714 G→T 6.8% H179Q (CAC→CAA)  PP_1407 ← sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,605,7140GT 6.8% 181.9 / 2.9 59H179Q (CAC→CAA) PP_1407sulfate transporter
Reads supporting (aligned to +/- strand):  ref base G (28/27);  new base T (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

CGGCAACGCCGAAGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAG  >  NC_002947/1605566‑1605861
                                                                                                                                                    |                                                                                                                                                    
cGGCAACGCCGAAGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAg                                                                                                                                                      >  5:203104/1‑149 (MQ=255)
    aaCGCCGAAGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGaa                                                                                                                                                  >  6:320256/1‑149 (MQ=255)
      cGCCGAAGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAAtt                                                                                                                                                >  3:167098/1‑149 (MQ=255)
      cGCCGAAGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAAtt                                                                                                                                                <  4:99670/149‑1 (MQ=255)
         cGAAGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTcc                                                                                                                                             <  4:167098/149‑1 (MQ=255)
           aaGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTcccc                                                                                                                                           >  1:25505/1‑149 (MQ=255)
              cGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAg                                                                                                                                        <  8:205767/149‑1 (MQ=255)
                gTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGgc                                                                                                                                      >  8:225282/1‑149 (MQ=255)
                 tCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGgcg                                                                                                                                     >  1:169284/1‑149 (MQ=255)
                        cGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑gg                                                                                                                                >  1:192920/1‑146 (MQ=255)
                            cAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGt                                                                                                                         <  2:137597/149‑1 (MQ=255)
                              ggCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTg                                                                                                                       <  2:25505/149‑1 (MQ=255)
                                 cAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGAt                                                                                                                     >  2:24325/1‑148 (MQ=255)
                                  aaCCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCa                                                                                                                   <  3:137004/149‑1 (MQ=255)
                                    ccAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCGAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAgg                                                                                                                 >  4:77337/1‑149 (MQ=255)
                                      agcagcATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCt                                                                                                               <  7:225282/149‑1 (MQ=255)
                                         agcaTGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCGCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAg                                                                                                            <  5:154414/149‑1 (MQ=255)
                                         agcaTGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAg                                                                                                            <  2:169284/149‑1 (MQ=255)
                                          gcaTGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGgcgc                                                                                                                                    <  5:282205/125‑1 (MQ=255)
                                          gcaTGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGgcgc                                                                                                                                    >  6:282205/1‑125 (MQ=255)
                                           caTGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGa                                                                                                          >  2:22670/1‑149 (MQ=255)
                                           caTGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGa                                                                                                          >  5:296193/1‑149 (MQ=255)
                                           caTGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGa                                                                                                          >  2:235622/1‑149 (MQ=255)
                                                    cgcCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTCACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGCATGAAAACGt                                                                                                 >  4:248604/1‑149 (MQ=255)
                                                     gcCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTc                                                                                                <  6:138773/149‑1 (MQ=255)
                                                           aCCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTc                                                                                          >  5:24041/1‑149 (MQ=255)
                                                              gTGAGGGCTGCGGGCAGGAGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGAt                                                                                       >  6:201478/1‑149 (MQ=255)
                                                                 aGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCtt                                                                                    <  1:22670/149‑1 (MQ=255)
                                                                  gggCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCTCGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCtt                                                                                    <  1:235622/149‑1 (MQ=255)
                                                                    gCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGaa                                                                                 <  1:24325/149‑1 (MQ=255)
                                                                        cGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGccc                                                                             <  5:320256/149‑1 (MQ=255)
                                                                               gcgcGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCATTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCt                                                                      >  3:6976/1‑149 (MQ=255)
                                                                               gcgcGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCATTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCt                                                                      >  3:2480/1‑149 (MQ=255)
                                                                               gcgcGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCt                                                                      >  5:160862/1‑149 (MQ=255)
                                                                                 gcgGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTg                                                                    <  6:182347/149‑1 (MQ=255)
                                                                                     gTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGa                                                                  >  3:272474/1‑147 (MQ=255)
                                                                                        gcgcCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCATTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCt                                                             <  4:2480/149‑1 (MQ=255)
                                                                                        gcgcCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCATTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCt                                                             <  4:6976/149‑1 (MQ=255)
                                                                                           ccGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGcc                                                          >  4:309710/1‑149 (MQ=255)
                                                                                                gCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGc                                                      >  5:109415/1‑148 (MQ=255)
                                                                                                  tgCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGAt                                                   <  1:279003/148‑1 (MQ=255)
                                                                                                  aGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCt                                                             >  4:243287/1‑139 (MQ=255)
                                                                                                  aGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCt                                                             <  3:243287/139‑1 (MQ=255)
                                                                                                    ccAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAg                                                                                <  3:224517/118‑1 (MQ=255)
                                                                                                    ccAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAg                                                                                >  4:224517/1‑118 (MQ=255)
                                                                                                     cAGATGACCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTACCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCCCCTTGGACAGCTGACAGGCGAtcag                                                >  4:281009/1‑149 (MQ=255)
                                                                                                      aGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAACCCCACCACCTTGGACAGCTGCCCGGCGAtcagg                                               >  2:321080/1‑149 (MQ=255)
                                                                                                         tGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCGTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGAGAtcaggct                                             >  8:180937/1‑148 (MQ=255)
                                                                                                           aGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGAtcaggctcag                                          >  3:177121/1‑149 (MQ=255)
                                                                                                             cAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGAtcaggctcaggc                                        <  4:272474/149‑1 (MQ=255)
                                                                                                             cAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGAtcaggctcaggc                                        >  1:186339/1‑149 (MQ=255)
                                                                                                                   cAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATAc                                  >  3:96304/1‑149 (MQ=255)
                                                                                                                            cAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCTGGCTCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGt                         <  5:13451/149‑1 (MQ=255)
                                                                                                                            cAGGCCGATGATCAATGTGACCCAGTGAATTACCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGt                         <  3:281009/149‑1 (MQ=255)
                                                                                                                                      atCAATGTGACCCAGTGAATTTCCCCCAGCCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCa                <  6:109415/148‑1 (MQ=255)
                                                                                                                                          aTGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATc           <  4:177121/149‑1 (MQ=255)
                                                                                                                                           tgtgACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCg          <  6:203104/149‑1 (MQ=255)
                                                                                                                                             tgACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGc        <  3:77337/149‑1 (MQ=255)
                                                                                                                                                 ccAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCcgc     >  1:292609/1‑148 (MQ=255)
                                                                                                                                                  cAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCga   <  6:296193/149‑1 (MQ=255)
                                                                                                                                                    gTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCgag  >  2:66764/1‑148 (MQ=255)
                                                                                                                                                    gTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCgag  >  7:256323/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGCAACGCCGAAGCGGTCGAGGCCGAACAGGCCAACCAGCAGCATGAACAGCGCCACTACCGTGAGGGCTGCGGGCAGGCGCGGGTAGCGCCGTGGCAGCCAGATGAGCAGGCCCAGGGCGGCCAGGCCGATGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCT‑GGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAG  >  NC_002947/1605566‑1605861

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