Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,368,199 G→C 7.7% intergenic (+28/‑155) PP_1193 → / → PP_1194 membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,368,1990GC 7.7% 193.1 / 3.0 63intergenic (+28/‑155)PP_1193/PP_1194membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base C (3/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGTGCCTCT  >  NC_002947/1368061‑1368345
                                                                                                                                           |                                                                                                                                                  
aaCCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGtttt                                                                                                                                          <  6:276819/149‑1 (MQ=255)
      gTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGa                                                                                                                                    <  8:190804/149‑1 (MQ=255)
       tGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAg                                                                                                                                   <  8:286814/149‑1 (MQ=255)
       tGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAg                                                                                                                                   >  6:147944/1‑149 (MQ=255)
         aCGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCa                                                                                                                                 <  8:98852/149‑1 (MQ=255)
           gATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAgg                                                                                                                               >  7:264973/1‑149 (MQ=255)
             tCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTa                                                                                                                                       >  8:259217/1‑139 (MQ=255)
             tCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTa                                                                                                                                       <  7:259217/139‑1 (MQ=255)
             tCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                             <  2:177298/149‑1 (MQ=255)
                    aaCCATGTGGGTGGCGAACCCGGTAAGTGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTATgg                                                                                                                                    >  1:3748/1‑135 (MQ=255)
                    aaCCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTATgg                                                                                                                                    <  2:3748/135‑2 (MQ=255)
                        aTGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCt                                                                                                                  <  8:264973/149‑1 (MQ=255)
                         tgtgGGTTGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCtt                                                                                                                 <  4:58795/149‑1 (MQ=255)
                         tgtgGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGCCCCGAAGTGGCCCATACAAACCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTGtt                                                                                                                 >  6:113420/1‑149 (MQ=255)
                          gtgggtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                             <  6:160821/136‑1 (MQ=255)
                          gtgggtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                             >  5:160821/1‑136 (MQ=255)
                              gtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCtt                                                                                                            <  2:293423/149‑1 (MQ=255)
                              gtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCtt                                                                                                            >  1:293423/1‑149 (MQ=255)
                               tggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCttt                                                                                                           <  1:60775/149‑1 (MQ=255)
                               tggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCtt                                                                                                            >  2:60775/1‑148 (MQ=255)
                                         ggTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGCCGGGGTTTGTATGAGCAGGGTCTGGCCTGGCTGGCCTTTATCACGGCag                                                                                                 >  2:128766/1‑148 (MQ=255)
                                         ggTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCAt                                                                                                 >  7:78749/1‑149 (MQ=255)
                                           tAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACTTCAttt                                                                                               >  4:246413/1‑149 (MQ=255)
                                                   cACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAAt                                                                                          >  1:150158/1‑147 (MQ=255)
                                                   cACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAAt                                                                                          >  8:175180/1‑147 (MQ=255)
                                                   cACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAAt                                                                                          >  3:133928/1‑147 (MQ=255)
                                                     cGGTGTTCATGCCGGCGGCGAGGCGGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACCAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAgc                                                                                      <  1:128766/149‑1 (MQ=255)
                                                     cGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTc                                                                                                                   >  2:299025/1‑120 (MQ=255)
                                                     cGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTc                                                                                                                   <  1:299025/120‑1 (MQ=255)
                                                             aTGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAg                                                                              >  5:313503/1‑149 (MQ=255)
                                                               gCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATc                                                                                                        <  8:126173/121‑1 (MQ=255)
                                                                          ggCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTc                                                                    <  3:85585/146‑1 (MQ=255)
                                                                          ggCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTc                                                                    >  4:85585/1‑146 (MQ=255)
                                                                            ccGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAgcg                                                                                     <  5:115375/127‑1 (MQ=255)
                                                                            ccGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAgcg                                                                                     >  6:115375/1‑127 (MQ=255)
                                                                             cgtggATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGTTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCAt                                                              <  5:249363/146‑1 (MQ=255)
                                                                                   tGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                       <  6:63639/134‑1 (MQ=255)
                                                                                   tGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                       >  5:63639/1‑134 (MQ=255)
                                                                                    gACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGc                                                                                                                        >  4:33055/1‑84 (MQ=255)
                                                                                    gACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGc                                                                                                                        <  3:33055/84‑1 (MQ=255)
                                                                                    gACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCggg                                                        <  7:175180/148‑1 (MQ=255)
                                                                                             agAATCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGc                                              <  2:78779/149‑1 (MQ=255)
                                                                                                 gCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGCCAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGtt                                          >  5:177337/1‑149 (MQ=255)
                                                                                                     aTTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCg                                      <  8:78749/149‑1 (MQ=255)
                                                                                                           gcTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGAt                                <  5:147944/149‑1 (MQ=255)
                                                                                                           gcTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGAt                                <  6:313503/149‑1 (MQ=255)
                                                                                                                 ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                       <  2:169968/104‑1 (MQ=255)
                                                                                                                 ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                       >  1:169968/1‑104 (MQ=255)
                                                                                                                 ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTgg                                   >  3:146470/1‑140 (MQ=255)
                                                                                                                 ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTgg                                   <  4:146470/140‑1 (MQ=255)
                                                                                                                 ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTc                          >  1:261553/1‑149 (MQ=255)
                                                                                                                 ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTc                          >  4:257460/1‑149 (MQ=255)
                                                                                                                          cccATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCAGTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATCATTGCGGTTTGCGTGGGATTAATTCCCGCGTCt                  >  5:225049/1‑148 (MQ=255)
                                                                                                                           ccATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTta                >  1:329125/1‑149 (MQ=255)
                                                                                                                             aTACAAAACCCGTCCGCTGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTtata              >  1:193510/1‑149 (MQ=255)
                                                                                                                             aTACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTta                >  3:196735/1‑147 (MQ=255)
                                                                                                                                 aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg          <  2:261553/149‑1 (MQ=255)
                                                                                                                                 aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg          >  5:133311/1‑149 (MQ=255)
                                                                                                                                 aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg          <  4:196735/149‑1 (MQ=255)
                                                                                                                                 aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGAGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg          >  6:281608/1‑149 (MQ=255)
                                                                                                                                   aaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGt        >  7:26897/1‑149 (MQ=255)
                                                                                                                                        aGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGTGCctc   >  4:167127/1‑149 (MQ=255)
                                                                                                                                         gCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCCTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGTGCctct  <  5:281608/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
AACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGTGCCTCT  >  NC_002947/1368061‑1368345

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