Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,072,218 T→A 6.5% I103F (ATC→TTC)  pstB‑II ← phosphate ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,072,2180TA 6.5% 180.4 / 2.4 62I103F (ATC→TTC) pstB‑IIphosphate ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (2/2);  total (33/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 = 9.79e-01

ATCGAGCACGCGCTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGTTCATCTGC  >  NC_002947/6072072‑6072366
                                                                                                                                                     |                                                                                                                                                    
aTCGAGCACGCGCTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGatat                                                                                                                                                    <  5:211600/149‑1 (MQ=255)
aTCGAGCACGCGCTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGatat                                                                                                                                                    >  7:111980/1‑149 (MQ=255)
       aCGCGCTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGcc                                                                                                                                             >  1:198453/1‑149 (MQ=255)
       aCGCGCTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGcc                                                                                                                                             <  2:198453/149‑1 (MQ=255)
       aCGCGCTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGcc                                                                                                                                             >  7:155120/1‑149 (MQ=255)
           gcTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCGGTa                                                                                                                                          >  4:22319/1‑148 (MQ=255)
                ttgttgATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGAAATTGTTGCCGTACAGGtt                                                                                                                                    >  6:38071/1‑149 (MQ=255)
                  gttgATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGtt                                                                                                                                    <  8:155120/147‑1 (MQ=255)
                     gATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACCCCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATga                                                                                                                               >  7:129569/1‑148 (MQ=255)
                     gATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGAt                                                                                                                                 <  6:92907/147‑1 (MQ=255)
                     gATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGAt                                                                                                                                 >  5:92907/1‑147 (MQ=255)
                     gATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGAt                                                                                                                                 <  8:129569/147‑1 (MQ=255)
                     gATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATgg                                                                                                                               >  2:61771/1‑149 (MQ=255)
                      aTGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGATAGATGGTCTTGGGGAACGGGTTCGGCTTATGGAACACCATGCCCACACGCCGGCGAATTTCGGCAACGTGTTCGCCGTTGCGGTAGATATTGTTGCCGTACAGGTTGATgg                                                                                                                               >  1:58373/1‑149 (MQ=255)
                        acccTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCAc                                                                                                                            >  6:288036/2‑149 (MQ=255)
                         cccTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCAc                                                                                                                            <  7:268279/148‑1 (MQ=255)
                              gATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTc                                                                                                                        <  2:58373/147‑1 (MQ=255)
                                  cgcAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCacac                                                                                                                   <  5:288036/148‑1 (MQ=255)
                                  cgcAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACg                                                                                                                  <  7:193950/149‑1 (MQ=255)
                                  cgcAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGACATTGTTGCCGTACAGGTTGATGGCACCTTCCACACg                                                                                                                  >  5:321569/1‑149 (MQ=255)
                                      ggCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGTAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGc                                                                                                                >  1:273452/1‑148 (MQ=255)
                                       gCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAAc                                                                                                             <  1:196876/149‑1 (MQ=255)
                                              aCCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAAc                                                                                                      <  1:61771/149‑1 (MQ=255)
                                              aCCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGAAATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAAc                                                                                                      <  5:38071/149‑1 (MQ=255)
                                                               tGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCgg                                                                                      >  7:257606/1‑149 (MQ=255)
                                                                   cTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTg                                                                                   >  3:114354/1‑148 (MQ=255)
                                                                    ttGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTg                                                                                   >  1:231596/1‑147 (MQ=255)
                                                                        ggAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGAATTTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGt                                                                                     <  1:286723/141‑1 (MQ=255)
                                                                        ggAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGAAATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGt                                                                                     >  2:286723/1‑141 (MQ=255)
                                                                            cGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCgca                                                                         <  8:111980/149‑1 (MQ=255)
                                                                               gTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCgcag                                                                        >  1:98705/1‑148 (MQ=255)
                                                                               gTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCgcag                                                                        >  8:187005/1‑147 (MQ=255)
                                                                               gTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCgcag                                                                        >  2:340921/1‑147 (MQ=255)
                                                                               ctatctcTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCgcagc                                                                       <  7:257888/142‑1 (MQ=255)
                                                                                   ggCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGt                                                                  <  7:187005/149‑1 (MQ=255)
                                                                                     cTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCg                                                                <  2:231596/149‑1 (MQ=255)
                                                                                      ttCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCgcagc                                                                       >  8:257888/1‑141 (MQ=255)
                                                                                             aaCACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACtt                                                            <  1:162073/146‑1 (MQ=255)
                                                                                             aaCACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACtt                                                            >  2:162073/1‑146 (MQ=255)
                                                                                                 ccATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAg                                                     >  3:239834/1‑149 (MQ=255)
                                                                                                  cATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGc                                                    <  4:239834/149‑1 (MQ=255)
                                                                                                       ccACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGACATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGAt                                               <  6:321569/149‑1 (MQ=255)
                                                                                                        cacaCGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAACGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATg                                              >  3:231549/1‑149 (MQ=255)
                                                                                                         acacGCCGGTGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGGCGACTTGCCGCAGCCGGATgg                                             >  8:60328/1‑149 (MQ=255)
                                                                                                         acacGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATgg                                             >  8:134335/1‑149 (MQ=255)
                                                                                                           acGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGc                                           >  5:61939/1‑149 (MQ=255)
                                                                                                           acGCCGGCGCAGTTCGGCAACGACTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGc                                           >  7:104570/1‑149 (MQ=255)
                                                                                                              ccGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGa                                        <  1:340921/149‑1 (MQ=255)
                                                                                                              ccGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGa                                        <  6:61939/149‑1 (MQ=255)
                                                                                                               cGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGAt                                       <  4:114354/149‑1 (MQ=255)
                                                                                                                     aGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGc                                 <  7:134335/149‑1 (MQ=255)
                                                                                                                      gTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCg                                >  5:120598/1‑149 (MQ=255)
                                                                                                                       ttCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCgg                               <  8:257606/149‑1 (MQ=255)
                                                                                                                           gCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCAc                           >  6:175331/1‑149 (MQ=255)
                                                                                                                                 tCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCgctg                     <  5:175331/149‑1 (MQ=255)
                                                                                                                                        cccTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCgg               >  4:131615/1‑149 (MQ=255)
                                                                                                                                        cccTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCgg               >  5:154453/1‑149 (MQ=255)
                                                                                                                                             gCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGAt            <  2:98705/147‑1 (MQ=255)
                                                                                                                                             gCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATg           <  2:273452/148‑1 (MQ=255)
                                                                                                                                               ggTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGtt         >  7:95288/1‑148 (MQ=255)
                                                                                                                                                 tAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAgg                                  <  6:28841/121‑1 (MQ=255)
                                                                                                                                                 tAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAgg                                  >  5:28841/1‑121 (MQ=255)
                                                                                                                                                 tAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGTTCAt      >  6:47695/1‑149 (MQ=255)
                                                                                                                                                 tAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGTTCAt      <  7:60328/149‑1 (MQ=255)
                                                                                                                                                    atatTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGTTCATCTg   >  8:67966/1‑149 (MQ=255)
                                                                                                                                                     tatTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGTTCATCTg   >  1:347141/1‑148 (MQ=255)
                                                                                                                                                     tatTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGTTCATCTGc  <  6:154453/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATCGAGCACGCGCTTCTTGTTGATGCCCTGGATGCGCAGGCCGTAGACCACGTTCTCGTAGATGGTCTTGGGGAACGGGTTCGGCTTCTGGAACACCATGCCCACACGCCGGCGCAGTTCGGCAACGTCTTCGCCCTTGCGGTAGATATTGTTGCCGTACAGGTTGATGGCACCTTCCACACGGCAACCGTCAACCAAGTCGTTCATCCGGTTGAAGGTGCGCAGCAGGGTCGACTTGCCGCAGCCGGATGGGCCGATGAAGGCGGTCACGCGCTGCTTCGGGATGTTCATCTGC  >  NC_002947/6072072‑6072366

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: