Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,255,852 A→G 6.4% intergenic (‑115/+21) PP_2852 ← / ← PP_2853 sulfatase domain‑containing protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,255,8520AG 6.4% 185.2 / 2.4 64intergenic (‑115/+21)PP_2852/PP_2853sulfatase domain‑containing protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (2/2);  total (30/34)
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.16e-01

GGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAA  >  NC_002947/3255716‑3255995
                                                                                                                                         |                                                                                                                                               
ggtggtCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCTGTGCGAACGGGCCCGGctttg                                                                                                                                     <  2:205471/149‑5 (MQ=255)
ggtggtCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGAt                                                                                                                                     <  2:326635/149‑1 (MQ=255)
 gtggtCGGACACGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATAAGGACAGCGCAGGGAAGGCGGCCAGTTACAAAGCTTGTAACCAAAGCCGGGCCCCTTCCCACAGGCCCGGCCGATg                                                                                                                                    >  8:301095/1‑149 (MQ=255)
         aCAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTAt                                                                                                                            >  8:176228/1‑149 (MQ=255)
             gTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAg                                                                                                                         >  2:170309/1‑148 (MQ=255)
             gTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATAAGGCCAGGGCAGGGAAGGAGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAg                                                                                                                         >  7:320983/1‑148 (MQ=255)
              ttcttcACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGaa                                                                                                                       <  4:27186/149‑1 (MQ=255)
                   cacaAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGt                                                                                                                     <  7:355783/146‑1 (MQ=255)
                    acaAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaat                                                                                                                  >  5:227083/1‑148 (MQ=255)
                     caAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaataa                                                                                                                <  2:344022/149‑1 (MQ=255)
                     caAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGtaat                                                                                                                  >  1:344022/1‑147 (MQ=255)
                      aaaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGATGGCGGCCAGTTACAAGGATTGTAAGCAAATCCGGGCCCGTTCGCTATGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                               <  7:301095/149‑1 (MQ=255)
                      aaaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                               >  5:127806/1‑149 (MQ=255)
                      aaaaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGACAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAAATATCAGAAGTAATAAg                                                                                                               >  5:183837/1‑149 (MQ=255)
                        aaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTAt                                                                                                                            >  3:163172/1‑134 (MQ=255)
                        aaaCTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTAt                                                                                                                            <  4:163172/134‑1 (MQ=255)
                               ccTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAAttt                                                                                                        >  3:319539/1‑147 (MQ=255)
                                     gCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTg                                                                                                 >  5:25981/1‑148 (MQ=255)
                                     gCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTg                                                                                                 >  8:88732/1‑148 (MQ=255)
                                       cATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                              <  7:109309/149‑1 (MQ=255)
                                       cATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                              <  5:201475/149‑1 (MQ=255)
                                        aTTAAATCGTCATGCATCTGTGCTTGATGATATCGAGGCCAGGGCAGGTAAGGCGCCCAGTTACAAGGATTGTAACAAAAGCCGGGCCCTTTCGCACAGGCCCGCCCGATGGAAACTATCAGAAGAAATAAGGGAATTACAGGCTgaa                                                                                               >  1:243046/1‑148 (MQ=255)
                                              tCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                               <  6:104425/125‑1 (MQ=255)
                                              tCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAg                                                                                                               >  5:104425/1‑125 (MQ=255)
                                                gTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCa                                                                                     <  4:319539/149‑1 (MQ=255)
                                                   aTGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAgg                                                                                   <  5:228413/148‑1 (MQ=255)
                                                    tGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTg                                                                                 <  7:88732/149‑1 (MQ=255)
                                                      cATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCg                                                                               <  7:176228/149‑1 (MQ=255)
                                                       aTCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTTCAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGt                                                                              <  6:127806/149‑1 (MQ=255)
                                                       aTCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGt                                                                              >  5:361689/1‑149 (MQ=255)
                                                        tCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGt                                                                                  <  5:29211/144‑1 (MQ=255)
                                                        tCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGt                                                                                  >  6:29211/1‑144 (MQ=255)
                                                                 gatgatATCTGGCCAGGGCATGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTc                                                                      <  8:320983/147‑1 (MQ=255)
                                                                  atgatATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGa                                                                   >  5:195576/1‑149 (MQ=255)
                                                                       atCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAg                                                                                    >  6:60453/1‑127 (MQ=255)
                                                                       atCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAg                                                                                    <  5:60453/127‑1 (MQ=255)
                                                                        tCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                              <  8:1971/116‑1 (MQ=255)
                                                                        tCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTgaaa                                                                                              >  7:1971/1‑116 (MQ=255)
                                                                        tCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATg                                                             <  6:25981/149‑1 (MQ=255)
                                                                           aGGCCAGGGCAGTGAAGGCTGCCAGTTACTAGGCTTGTAACCAAAGCAGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                                    <  2:158170/108‑1 (MQ=255)
                                                                           aGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAgg                                                                                                    >  1:158170/1‑108 (MQ=255)
                                                                               cagggcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAgg                                                       >  4:374349/1‑149 (MQ=255)
                                                                                  ggcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTg                                                    >  8:215600/1‑149 (MQ=255)
                                                                                   gcagggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTgg                                                   <  6:227083/149‑1 (MQ=255)
                                                                                       ggAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCAc                                               >  6:131853/1‑149 (MQ=255)
                                                                                        gAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCAc                                               <  5:131853/148‑1 (MQ=255)
                                                                                         aaGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGa                                             <  1:170309/149‑1 (MQ=255)
                                                                                           ggcggcCAGTTACAAGGGTTGTAACCAAAGCCGGGCGCGTTAGCTCATGCCCGGCCGATGGCAACTATCAGAAGTAATATGGTAATTTCAGGCTGAAAGAAAAGTCAGGGGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGAtt                                           <  6:183837/149‑2 (MQ=255)
                                                                                           ggcggcCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATg                                           <  1:63667/149‑1 (MQ=255)
                                                                                               tttggTTACAAGCCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGt                                         >  1:378742/5‑147 (MQ=255)
                                                                                                 cAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAgg                                     <  6:361689/149‑1 (MQ=255)
                                                                                                 cAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGATGGGCACGATGGTCAgg                                     >  6:208700/1‑149 (MQ=255)
                                                                                                         aGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGt                             >  2:201812/1‑149 (MQ=255)
                                                                                                           ccTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTgg                           <  2:378742/148‑1 (MQ=255)
                                                                                                            cTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTggc                          >  1:161204/1‑149 (MQ=255)
                                                                                                               gTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTggcggc                       >  8:330170/1‑149 (MQ=255)
                                                                                                                   ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAg                   <  7:330170/149‑1 (MQ=255)
                                                                                                                   ccAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAg                   <  5:208700/149‑1 (MQ=255)
                                                                                                                      aaGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTc                <  6:195576/149‑1 (MQ=255)
                                                                                                                          cGGGCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGt                                    <  2:14971/125‑1 (MQ=255)
                                                                                                                          cGGGCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGt                                    >  1:14971/1‑125 (MQ=255)
                                                                                                                            ggCCTGTGCGAACGGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGcc           >  6:319661/1‑148 (MQ=255)
                                                                                                                               ccGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAata       <  2:243046/149‑1 (MQ=255)
                                                                                                                                     gCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGaa  <  1:201812/148‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                               
GGTGGTCGGACAGGTTCTTCACAAAAACTTCCCTTTTGCCATTAAATCGTCATGCATCTGTGCTTGATGATATCAGGCCAGGGCAGGGAAGGCGGCCAGTTACAAGGCTTGTAACCAAAGCCGGGCCCGTTCGCACAGGCCCGGCCGATGGCAACTATCAGAAGTAATAAGGGAATTTCAGGCTGAAAGAAAAGTCAGGTGCGTCATCGGTCAGACCAATGGACAGGTTGGGCACGATGGTCAGGTTGTCGGTGGCGGCGAAGGTCAGGCCAATATTGAA  >  NC_002947/3255716‑3255995

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