Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 274,145 A→T 5.9% L142Q (CTG→CAG)  PP_0221 ← methionine ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947274,1450AT 5.9% 208.3 / 2.7 68L142Q (CTG→CAG) PP_0221methionine ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (30/33);  new base T (2/2);  total (33/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 = 6.53e-01

CACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCC  >  NC_002947/274000‑274293
                                                                                                                                                   |                                                                                                                                                         
cACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagca                                                                                                                                                        <  4:194925/149‑1 (MQ=255)
cACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagca                                                                                                                                                        <  4:201672/149‑1 (MQ=255)
cACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagc                                                                                                                                                         <  4:311486/148‑1 (MQ=255)
     ggTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCTCCTCGACCTGAAGGATCTTAAGCTTTTTAGGGTTGGCGAGTATGTCAGCCTAGGTGTACTTGTAGCACACGCCGGGATTGATGGTGATCAACCCGGCCTTGGCcagcagcag                                                                                                                                                    >  8:203036/1‑148 (MQ=255)
       tCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagcagcag                                                                                                                                                 <  3:68953/149‑1 (MQ=255)
         ggCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGc                                                                                                                                               >  8:329053/1‑149 (MQ=255)
            gtcgtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGATGTATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCaa                                                                                                                                            <  2:314053/149‑2 (MQ=255)
               gtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagca                                                                                                                                                     <  3:219749/137‑1 (MQ=255)
               gtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCcagcagca                                                                                                                                                     >  4:219749/1‑137 (MQ=255)
               gtcgtAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGc                                                                                                                                         >  4:413188/1‑149 (MQ=255)
                     ggCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACAGCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCggg                                                                                                                                     <  6:182733/148‑3 (MQ=37)
                     ggCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGt                                                                                                                                  >  2:270691/1‑149 (MQ=255)
                      gcgcGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGtt                                                                                                                                 >  3:267145/1‑149 (MQ=255)
                           tACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACAGCGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCGGGTTGATca                                                                                                                            <  8:242476/149‑2 (MQ=25)
                             ccAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCggg                                                                                                                          >  2:270644/1‑149 (MQ=255)
                                    cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTgg                                                                                                                   <  3:77184/149‑1 (MQ=255)
                                    cACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCCGGTCGTTgg                                                                                                                   >  7:228909/1‑149 (MQ=255)
                                        gCCTCGACCTGAAGGATCATCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGAt                                                                                                               <  5:12758/149‑1 (MQ=255)
                                           tCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGc                                                                                                            >  5:241528/1‑149 (MQ=255)
                                            cGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGc                                                                                                            >  2:25797/1‑148 (MQ=255)
                                             gACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGc                                                                                                            <  1:270691/147‑1 (MQ=255)
                                                ccgAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATc                                                                                                                             <  1:332645/125‑1 (MQ=255)
                                                ccgAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATc                                                                                                                             >  2:332645/3‑127 (MQ=255)
                                                   aaGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTgg                                                                                                    >  7:391253/1‑149 (MQ=255)
                                                   aaGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTgg                                                                                                    >  1:48334/1‑149 (MQ=255)
                                                    aGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGgc                                                                                                   <  2:135035/149‑1 (MQ=255)
                                                        tCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCt                                                                                               >  6:139130/1‑149 (MQ=255)
                                                          ttCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCAGGGCTTGAGGGTGATCAACACGGCCTTGGCCAGCATCAGCAGGCCACGGCAGATGTTGATCTGGTCGTTGGCGATGGCCACCCTGGCGCCtt                                                                                              >  5:338347/1‑148 (MQ=255)
                                                           tCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCtt                                                                                              >  4:240069/1‑147 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGTGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGTCGCCtt                                                                                              <  7:203036/146‑1 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGATTAAAAAACCCGCCTTTGGCAAGAGGAGCCCGCCACGGGCGGGGTTTGAACGGTTCTTTGGCGATGGCCACGATgccccat                                                                                               >  1:131063/1‑141 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGcc                                                                                                <  6:199104/144‑1 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGcc                                                                                                <  4:267145/144‑1 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCtt                                                                                              <  4:205913/146‑1 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCt                                                                                               <  7:163102/145‑1 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTc                                                                                             <  3:426178/147‑1 (MQ=255)
                                                            cAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCa                                                                                            <  3:240069/148‑1 (MQ=255)
                                                             aGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAAGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAg                                                                                           <  2:423271/148‑1 (MQ=255)
                                                              gCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCATTGTAGCCCACGCCGGGCTTGAGGGTGATCAAACCGGCCTTGGCCAGCAGAAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAgg                                                                                          >  7:293551/1‑148 (MQ=255)
                                                                      gggTTGGCGAGGATGTCAGCCTCGGTGGGCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTCGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCgg                                                                              <  7:80512/148‑1 (MQ=255)
                                                                       ggTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGc                                                                             >  2:348727/1‑148 (MQ=255)
                                                                         ttGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCa                                                                            >  4:62431/1‑147 (MQ=255)
                                                                               aGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCtttt                                                                    <  5:127518/149‑1 (MQ=255)
                                                                                ggATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTg                                                                   <  3:399424/149‑1 (MQ=255)
                                                                                   tggcAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTAc                                                                <  3:246977/146‑1 (MQ=255)
                                                                                       aGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACtttt                                                            <  7:329053/149‑1 (MQ=255)
                                                                                         ccTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTc                                                          >  7:290557/1‑149 (MQ=255)
                                                                                              gTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                     >  6:249400/1‑149 (MQ=255)
                                                                                               tGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                     >  6:151034/1‑148 (MQ=255)
                                                                                               tGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                     >  7:212331/1‑148 (MQ=255)
                                                                                                ggCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTa                                                     >  5:192648/1‑147 (MQ=255)
                                                                                                 gCCTTGTAGCCGACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAg                                                  <  5:160233/149‑1 (MQ=255)
                                                                                                 gCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCCGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAg                                                  >  6:266241/1‑149 (MQ=255)
                                                                                                  ccTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAgg                                                 >  5:178949/1‑149 (MQ=255)
                                                                                                   cTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGc                                                <  3:413188/149‑1 (MQ=255)
                                                                                                      gTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGcc                                               <  5:247121/147‑1 (MQ=255)
                                                                                                           ccACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttg                                        >  7:29753/1‑149 (MQ=255)
                                                                                                            gaCGCCTGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCATCAGGCCACGGACGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACgttgt                                       <  8:293551/148‑1 (MQ=255)
                                                                                                                ccGGGCTTGAGGGGGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgat                                   <  6:241528/149‑1 (MQ=255)
                                                                                                                 cGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTtgatg                                  <  4:9791/149‑1 (MQ=255)
                                                                                                                    gCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGatgatg                               >  2:43160/1‑149 (MQ=255)
                                                                                                                          gggTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                         <  6:192648/149‑1 (MQ=255)
                                                                                                                            gTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGGCTCA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCggg                         <  1:343246/149‑1 (MQ=255)
                                                                                                                                aTCAACCCGGGCTTTGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTTATGATGCCCGGGGCGt                     <  6:338347/148‑1 (MQ=255)
                                                                                                                                 tCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGt                     >  8:228068/1‑147 (MQ=255)
                                                                                                                                   aaCCCGGCCTTGGCCAGCCGCAGCAGGCCACGGCCGCTTGTGATCGGGTCGTTGGCGTTGGCCACGCTGGCGCCTTCAGG‑‑CACTGCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTaa                   <  1:348727/149‑1 (MQ=255)
                                                                                                                                     cccGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGca               >  3:37774/1‑149 (MQ=255)
                                                                                                                                       gcGGCCGTGGCCTGCTGCTGCTGGCCAAGGCCGGGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcac              >  8:343721/2‑148 (MQ=38)
                                                                                                                                        cGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcaca             <  6:178949/149‑1 (MQ=255)
                                                                                                                                        cGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGcac              >  1:212053/1‑148 (MQ=255)
                                                                                                                                            cgtGGCCTGCTGCTGCTGGCCAAGGCCGGGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGtt        >  7:288013/3‑149 (MQ=38)
                                                                                                                                               ggCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTg       <  7:228068/148‑1 (MQ=255)
                                                                                                                                                gCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGaa     <  8:29753/149‑1 (MQ=255)
                                                                                                                                                 ccAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAgg                 <  2:128553/136‑1 (MQ=255)
                                                                                                                                                 ccAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAgg                 >  1:128553/1‑136 (MQ=255)
                                                                                                                                                  cagcagcagcagGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGc   <  8:212331/149‑1 (MQ=255)
                                                                                                                                                  cagcagcagcagGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGc   >  8:53940/1‑149 (MQ=255)
                                                                                                                                                   agcagcagcagGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGcc  <  2:131063/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
CACCAGGTCGGCGTCGTCGTAGGCGCGTACCAGTTGCACCGCCTCGACCTGAAGGATCTTCAGCTTTTTAGGGTTGGCGAGGATGTCAGCCTCGGTGGCCTTGTAGCCCACGCCGGGCTTGAGGGTGATCAACCCGGCCTTGGCCAGCAGCAGCAGGCCACGGCCGCTGTTGATCGGGTCGTTGGCGATGGCCACGCTGGCGCCTTCAGG‑‑CA‑‑GCTCGGCAAAGCTTTTGTACTTTTTCGAGTACAGGCCAACGTTGTTGATGATGCCCGGGGCGTAAGGCACAAGGTTGAAGCC  >  NC_002947/274000‑274293

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: