Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,379,582 T→C 6.6% intergenic (‑244/+179) PP_3855 ← / ← PP_3856 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,379,5820TC 6.6% 173.1 / 2.2 60intergenic (‑244/+179)PP_3855/PP_3856hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/33);  new base C (2/2);  total (25/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 = 8.59e-01

GGCGGCTATGAGCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCGAAGAATCGCATCGGCCTCCTCCAGACG  >  NC_002947/4379435‑4379727
                                                                                                                                                   |                                                                                                                                                 
ggcggcTATGAGCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCtata                                                                                                                                                  <  4:262260/149‑1 (MQ=255)
ggcggcTATGAGCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCtat                                                                                                                                                   >  7:222855/1‑148 (MQ=255)
 gcggcTATGAGCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATaa                                                                                                                                                 <  5:97775/149‑1 (MQ=255)
           gCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCTCAc                                                                                                                                       >  7:199669/1‑149 (MQ=255)
           gCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCAc                                                                                                                                       >  7:177162/1‑149 (MQ=255)
           gCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCAc                                                                                                                                       <  5:157047/149‑1 (MQ=255)
            cAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCAc                                                                                                                                       >  2:112159/1‑148 (MQ=255)
               aGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTACAACGACGCGCAcc                                                                                                                                      >  2:63562/1‑146 (MQ=255)
                 gACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAagag                                                                                                                                 <  2:61279/149‑1 (MQ=255)
                      ttGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGc                                                                                                                             <  7:87126/148‑1 (MQ=255)
                        ggCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCtt                                                                                                                           <  5:267178/148‑1 (MQ=255)
                             gggtAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGTCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAAc                                                                                                                     <  8:130732/146‑1 (MQ=255)
                              gCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACa                                                                                                                    <  8:222855/149‑1 (MQ=255)
                                     gTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCtttt                                                                                                             >  1:160068/1‑149 (MQ=255)
                                     gTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCtttt                                                                                                             <  1:112159/149‑1 (MQ=255)
                                     gTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTACAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCttt                                                                                                              <  1:63562/148‑1 (MQ=255)
                                      tCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCttttt                                                                                                            <  6:276102/149‑1 (MQ=255)
                                       cGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCttttt                                                                                                            >  3:150646/1‑148 (MQ=255)
                                       cGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATAGGCGAAAAATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCGATAAACAAGCAACCCAACCGGGGGTTATGAACAGGCTTTTtg                                                                                                           >  2:310690/1‑149 (MQ=255)
                                        gTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTtga                                                                                                          >  4:2959/1‑149 (MQ=255)
                                        gTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTtg                                                                                                           <  2:160068/148‑1 (MQ=255)
                                          tcttctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTtgat                                                                                                         >  3:185108/1‑148 (MQ=255)
                                            ttctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACgcgc                                                                                                                                         <  8:183356/114‑1 (MQ=255)
                                            ttctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACgcgc                                                                                                                                         >  7:183356/1‑114 (MQ=255)
                                            atctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGAGATGCCACTCCTCCTAACCAGATAGGGAGCTACAACGACGCGCACCAAGAGTGGCTTCTGAACAGGCTTTTtgatgat                                                                                                      >  6:323871/2‑149 (MQ=255)
                                             tctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCTAACCAGCTAGGGTGCTATAATGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTtgatgat                                                                                                      <  7:155687/148‑1 (MQ=255)
                                             tctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTtgatgat                                                                                                      <  4:153979/148‑1 (MQ=255)
                                              ctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTtgatgat                                                                                                      <  3:151394/147‑1 (MQ=255)
                                              ctATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTtgatgat                                                                                                      <  8:31030/147‑1 (MQ=255)
                                                 tCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACt                                                                                                  >  2:87638/1‑148 (MQ=255)
                                                    ggCTGCAATCACTACTTCTATAGCCCGTATGGTAGTCGTTGGGATCAGCGAAACATACCCTTCTATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGc                                                                                              <  7:109778/149‑1 (MQ=255)
                                                         cAATCTCTACTTCAATAGCCCGTATGGGAGTCTTTTGGATCAGCGAAACATTCCGTGCGATGACTCTCCTCCGAACGAGCTAGGGTGCTACAACGACGCGCAGCAAGAGGGGCTTCTGAACATGCTTTTTGATGATGACTAAGCAAt                                                                                           <  5:323871/147‑1 (MQ=255)
                                                         cAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAAt                                                                                           <  3:2959/147‑1 (MQ=255)
                                                                tactTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAgc                                                                                  >  7:213933/1‑149 (MQ=255)
                                                                          gCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGc                                                                         <  8:177162/148‑1 (MQ=255)
                                                                              gTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCt                                                                    <  8:213933/149‑1 (MQ=255)
                                                                                   ggAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTACGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGt                                                                <  8:199669/148‑1 (MQ=255)
                                                                                            gggATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCtt                                                      >  7:23585/1‑149 (MQ=255)
                                                                                                tCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTttcat                                                  <  4:150646/149‑1 (MQ=255)
                                                                                                tCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTttcat                                                  >  2:97258/1‑149 (MQ=255)
                                                                                                tCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTttcat                                                  >  6:49169/1‑149 (MQ=255)
                                                                                                    cGAAACATACCCTGCGATGCCACTCCTCCTAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTtcattcat                                              <  1:310690/149‑1 (MQ=255)
                                                                                                      aaaCATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCt                                                       >  2:208325/1‑138 (MQ=255)
                                                                                                      aaaCATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCt                                                       <  1:208325/138‑1 (MQ=255)
                                                                                                       aaCATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAagag                                                                                                                                 <  4:41132/63‑1 (MQ=255)
                                                                                                       aaCATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAagag                                                                                                                                 >  3:41132/1‑63 (MQ=255)
                                                                                                          aTACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACTTTTCATTCATTGCTCa                                        <  1:87638/149‑1 (MQ=255)
                                                                                                          aTACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTttca                                                   <  6:258512/138‑1 (MQ=255)
                                                                                                          aTACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTttca                                                   >  5:258512/1‑138 (MQ=255)
                                                                                                             cccTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAgg                                      <  8:23585/148‑1 (MQ=255)
                                                                                                               cTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTATCATTCATTCCTAAGGtt                                    >  5:310686/1‑148 (MQ=255)
                                                                                                               cTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTGCTGATCAGGCTTTTTGATGATGACTACGCACTCAGTCCCGCTTGCTTGGCGGGGTTCGTTGGAGAGCTTTCATTCATTGCTCAGGtg                                    >  4:63535/1‑147 (MQ=255)
                                                                                                                  cGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTAt                                 <  4:185108/148‑1 (MQ=255)
                                                                                                                       ccACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCga                            <  7:288966/148‑1 (MQ=255)
                                                                                                                       ccACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCga                            <  5:89223/148‑1 (MQ=255)
                                                                                                                        cACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCgaa                           <  5:49169/148‑1 (MQ=255)
                                                                                                                        cACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCg                             >  5:290506/1‑146 (MQ=255)
                                                                                                                             ctccGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCGAAGAATc                      <  6:290506/148‑1 (MQ=255)
                                                                                                                                 gAACCAGCTAGGTTTCTATAACGACGCGCACCAAGAGTGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTccattcat                                              <  4:17279/120‑1 (MQ=255)
                                                                                                                                 gAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCGAAGAATCGCATc                 >  4:94393/1‑149 (MQ=255)
                                                                                                                                 gAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGAGACCTTccattcat                                              >  3:17279/1‑120 (MQ=255)
                                                                                                                                     cAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCGAAGAATCGCATCGGc              >  8:279640/1‑148 (MQ=255)
                                                                                                                                                 tataACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGATTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCGAAGAATCGCATCGGCCTCCTCCAGACg  >  5:309548/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGCGGCTATGAGCAAAGGACTTTTGGCAGGGCTAGTGGTCGTTCTTCTATCGGGCTGCAATCACTACTTCAATAGCCCGTATGGGAGTCGTTGGGATCAGCGAAACATACCCTGCGATGCCACTCCTCCGAACCAGCTAGGGTGCTATAACGACGCGCACCAAGAGGGGCTTCTGAACAGGCTTTTTGATGATGACTAAGCAATCAGTCCAGCTTGCTTGGCGGGCTTCGTTGGCGACCTTTCATTCATTGCTCAGGTTTATGGCGAAGAATCGCATCGGCCTCCTCCAGACG  >  NC_002947/4379435‑4379727

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