Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,489,780 A→C 6.3% M88L (ATG→CTG) ‡ PP_4825 → MarC family membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,489,7800AC 6.3% 240.2 / 6.8 79M88L (ATG→CTG) ‡PP_4825MarC family membrane protein
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base C (3/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

GCAACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAAC  >  NC_002947/5489632‑5489927
                                                                                                                                                         |                                                                                                                                                    
gcaaCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCa                                                                                                                                                      <  4:312742/149‑1 (MQ=255)
  aaCTGGCCTGGAAGGTCGCGTTCACAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATg                                                                                                                                                    >  2:40678/1‑149 (MQ=255)
   aCTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATgg                                                                                                                                                   <  2:374786/149‑1 (MQ=255)
       gCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGccc                                                                                                                                                <  3:147311/148‑1 (MQ=255)
        ccTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCTGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAg                                                                                                                                              <  3:151528/149‑1 (MQ=255)
         cTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAgg                                                                                                                                             <  8:197011/149‑1 (MQ=255)
            gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCCATGCCCAggg                                                                                                                                            <  1:342576/147‑1 (MQ=255)
            gAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCCATGCCCAGGGcc                                                                                                                                          <  7:261483/149‑2 (MQ=255)
                     gTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCcggcgg                                                                                                                                 >  5:352963/1‑149 (MQ=255)
                     gTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCcggcgg                                                                                                                                 >  2:126241/1‑149 (MQ=255)
                         gCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACa                                                                                                                             <  1:104816/149‑1 (MQ=255)
                          cAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAg                                                                                                                            <  7:25358/149‑1 (MQ=255)
                           aGCCGTGGGCGCCAGCGGGTTGCGGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCTATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAgg                                                                                                                           <  4:315855/149‑1 (MQ=255)
                              cTTGGGCGCCAGCGTGTTGTTGTACCTGTTCGGGCGGGCTATCTTCGTAATTTTCGGCAGCACCGCCTATGCCTTCCGCAGCGTTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAg                                                                                                                        <  1:336305/149‑1 (MQ=255)
                                     gCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACg                                                                                                                  <  2:238863/148‑1 (MQ=255)
                                     cccAGCGTGTTGCTGTACCTGTTCGGGCGGGTGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGt                                                                                                                 >  7:164444/2‑149 (MQ=255)
                                       cAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGtt                                                                                                                >  2:344809/1‑148 (MQ=255)
                                       cAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGtt                                                                                                                >  3:85216/1‑148 (MQ=255)
                                       cAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTc                                                                                                               <  5:81312/149‑1 (MQ=255)
                                        aGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGCCATCACCGCCGATGCCTGCCGCATGGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGCCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTca                                                                                                              <  3:214034/149‑1 (MQ=255)
                                           gtgtTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcagca                                                                                                           <  8:156409/149‑1 (MQ=255)
                                                gCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATTGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGAt                                                                                                       >  1:284373/1‑149 (MQ=255)
                                                    tACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCTATGCCTTCCGAATGGGTGCTGGCAGCGGGCTTTTCATTTC‑‑‑‑GGCCCTGGGCAGGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGa                                                                                                  <  8:309091/148‑1 (MQ=255)
                                                      ccTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCa                                                                                               >  2:237381/1‑149 (MQ=255)
                                                            cgggcgggCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAc                                                                                         >  4:32101/1‑149 (MQ=255)
                                                             gggcgggcGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAcc                                                                                        >  6:134571/1‑149 (MQ=255)
                                                               gcgggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCACGGCCCAGGGCAAGCCGGCGGt                                                                                                                                <  7:115598/108‑1 (MQ=255)
                                                               gcgggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCAAGGCCCAGGGCAAGCCGGCGGt                                                                                                                                >  8:115598/1‑108 (MQ=255)
                                                                  ggcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGa                                                                                   <  6:96776/149‑1 (MQ=255)
                                                                  ggcgAGCTGCGTAATTTGCGGCACCCCCGCCGATGCCTTCCGAAGGGGTGCGGGCCGCGTGCGGTTCATTTC‑‑‑‑GGCCCGGGGCTTGGCCCCGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGa                                                                                   <  2:152331/149‑1 (MQ=255)
                                                                   gcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGtt                                                                                                                <  5:375422/120‑1 (MQ=255)
                                                                   gcgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGtt                                                                                                                >  6:375422/1‑120 (MQ=255)
                                                                    cgATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGAc                                                                                  >  4:142335/1‑148 (MQ=255)
                                                                     gATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCa                                                                                >  6:150075/1‑149 (MQ=255)
                                                                     gATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCa                                                                                >  5:282718/1‑149 (MQ=255)
                                                                        cTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAc                                                                                         >  1:210133/1‑137 (MQ=255)
                                                                        cTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTAc                                                                                         <  2:210133/137‑1 (MQ=255)
                                                                         ttCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTcc                                                                            <  5:11932/149‑1 (MQ=255)
                                                                         ttCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTcc                                                                            >  5:137648/1‑149 (MQ=255)
                                                                            ggAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                         <  1:344497/149‑1 (MQ=255)
                                                                            ggAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCt                                                                         <  6:80166/149‑1 (MQ=255)
                                                                                ttttCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATCGTGACCATCGTACCGCTGACCATTCCGCTGAc                                                                      <  1:237381/149‑1 (MQ=255)
                                                                                 tttCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCg                                                                    >  8:54001/1‑149 (MQ=255)
                                                                                    cGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGtt                                                                  <  1:344809/148‑1 (MQ=255)
                                                                                      gCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGc                                                               <  2:284373/149‑1 (MQ=255)
                                                                                          cACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCg                                                           >  2:21536/1‑149 (MQ=255)
                                                                                           aCCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCgg                                                          <  7:54001/149‑1 (MQ=255)
                                                                                            ccgccgATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGc                                                         >  2:115281/1‑149 (MQ=255)
                                                                                              gccgATGCCTTCCGCATCGGGGCTGGCAGCGTGCTGTTCATTTCCTTTACCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCg                                                           <  1:126241/149‑1 (MQ=255)
                                                                                                      ccTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGt                                                >  2:182354/1‑149 (MQ=255)
                                                                                                       cTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTg                                               <  6:352963/149‑1 (MQ=255)
                                                                                                        ttCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTcag                                                                                                             >  5:41007/1‑87 (MQ=255)
                                                                                                            gCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGt                                                                   <  6:156434/124‑1 (MQ=255)
                                                                                                            gCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGt                                                                   >  5:156434/1‑124 (MQ=255)
                                                                                                             cATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCAGGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgc                                         <  4:85216/149‑1 (MQ=255)
                                                                                                               tCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCctgctg                                       <  3:32101/149‑1 (MQ=255)
                                                                                                                  gTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGt                                     >  2:70916/1‑148 (MQ=255)
                                                                                                                      tGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgg                                <  3:142335/149‑1 (MQ=255)
                                                                                                                        gCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgggg                              <  5:134571/149‑1 (MQ=255)
                                                                                                                        gCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgg                                <  5:150075/147‑1 (MQ=255)
                                                                                                                             gTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGag                         <  1:70916/149‑1 (MQ=255)
                                                                                                                                 tGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCGCTGCTGGTGATGGGGGTGagccag                     >  5:253325/1‑149 (MQ=255)
                                                                                                                                 tGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGagccag                     <  6:137648/149‑1 (MQ=255)
                                                                                                                                    tCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCt                                                                                     >  6:133437/1‑82 (MQ=255)
                                                                                                                                    tCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCt                                                                                     <  5:133437/82‑1 (MQ=255)
                                                                                                                                      aTTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCa                >  6:244372/1‑149 (MQ=255)
                                                                                                                                      aTTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCa                >  6:249812/1‑149 (MQ=255)
                                                                                                                                       tttC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgg                                <  8:87170/132‑1 (MQ=255)
                                                                                                                                       tttC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATgg                                >  7:87170/1‑132 (MQ=255)
                                                                                                                                       tttC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCAc               >  3:363564/1‑149 (MQ=255)
                                                                                                                                       tttC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCAc               >  3:66871/1‑149 (MQ=255)
                                                                                                                                            gcttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACt              >  6:248737/5‑149 (MQ=255)
                                                                                                                                              ttgCCCTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTgg            >  3:128660/3‑149 (MQ=255)
                                                                                                                                               ggCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTggg           <  7:207844/149‑1 (MQ=255)
                                                                                                                                                gCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACCTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGa          >  1:266653/1‑149 (MQ=255)
                                                                                                                                                 cccTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt         <  8:164444/149‑1 (MQ=255)
                                                                                                                                                  ccTGGGCCATGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGAt         >  8:34845/1‑148 (MQ=255)
                                                                                                                                                      ggCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaaa     <  4:66871/148‑1 (MQ=255)
                                                                                                                                                       gCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGaaaa    >  2:105421/1‑148 (MQ=255)
                                                                                                                                                        cATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAAc  >  1:189083/1‑149 (MQ=255)
                                                                                                                                                         |                                                                                                                                                    
GCAACTGGCCTGGAAGGTCGCGTTCGCAGCCTTGGTCGCCAGCGTGTTGCTGTACCTGTTCGGGCGGGCGATCTTCGGAATTTTCGGCATCACCGCCGATGCCTTCCGCATCGGTGCTGGCAGCGTGCTGTTCATTTC‑‑‑‑GGCCCTGGGCATGGCCCAGGGCAAGCCGGCGGTACAGGCAGACAACGTTCAGCAGGATGTGACCATCGTACCGCTGACCATTCCGCTGACCGTTGGCCCCGGCACCATCGGTGCCCTGCTGGTGATGGGGGTGAGCCAGCCGCACTGGGATGAAAAAC  >  NC_002947/5489632‑5489927

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: