Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,788,401 G→C 7.2% R78R (CGC→CGG)  PP_3350 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,788,4010GC 7.2% 212.3 / 3.9 66R78R (CGC→CGG) PP_3350hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base C (2/3);  total (33/33)
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

GGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCT  >  NC_002947/3788257‑3788549
                                                                                                                                                |                                                                                                                                                    
ggTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTTCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCggg                                                                                                                                                  <  3:275457/149‑1 (MQ=255)
 gTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGc                                                                                                                                                 <  3:131330/149‑1 (MQ=255)
    ttCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTc                                                                                                                                              <  5:219095/149‑1 (MQ=255)
     tctcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGACGCATCGGCCCGGTAGCGGGCCt                                                                                                                                               >  5:128623/1‑147 (MQ=255)
       tcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTACTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGc                                                                                                                                           >  5:305358/1‑149 (MQ=255)
         gACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGccc                                                                                                                                         >  2:108005/1‑149 (MQ=255)
         gACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccgatgc                                                                                                                                           <  4:124113/147‑6 (MQ=255)
          aCCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGcccc                                                                                                                                        <  7:169557/149‑1 (MQ=255)
            cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCt                                                                                                                                               <  5:326679/140‑1 (MQ=255)
            cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCt                                                                                                                                               >  6:326679/1‑140 (MQ=255)
            cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGa                                                                                                                                      <  7:40319/149‑1 (MQ=255)
            cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccga                                                                                                                                              >  6:112844/1‑139 (MQ=255)
            cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGCTACCGGGccga                                                                                                                                              <  5:112844/141‑3 (MQ=255)
               gACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGtt                                                                                                                                   >  6:68687/1‑149 (MQ=255)
               gACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGtt                                                                                                                                   >  4:224211/1‑149 (MQ=255)
                 cGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGc                                                                                                                                 <  3:196618/149‑1 (MQ=255)
                      gTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCgg                                                                                                                                                   <  1:43168/126‑1 (MQ=255)
                      gTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGACCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCgg                                                                                                                                                   >  2:43168/1‑126 (MQ=255)
                          gCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGcc                                                                                                                        >  7:125487/1‑149 (MQ=255)
                               gCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGActc                                                                                                                   <  1:129348/149‑1 (MQ=255)
                                       gAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCTGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGcc                                                                                                           >  4:105030/1‑149 (MQ=255)
                                             gTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATgg                                                                                                      >  2:214393/1‑148 (MQ=255)
                                             gTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATgg                                                                                                      >  3:196475/1‑148 (MQ=255)
                                             gTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCAGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATgg                                                                                                      >  6:168375/1‑148 (MQ=255)
                                              tGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTa                                                                                                    <  8:125487/149‑1 (MQ=255)
                                                      ggCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAgg                                                                                             <  5:168375/148‑1 (MQ=255)
                                                       gCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGc                                                                                                                                                 >  6:67427/1‑95 (MQ=255)
                                                       gCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGc                                                                                                                                                 <  5:67427/95‑1 (MQ=255)
                                                            ggCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCg                                                                                       >  6:24851/1‑148 (MQ=255)
                                                              ccTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGaaa                                                                                    <  2:109302/149‑1 (MQ=255)
                                                               cTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAg                                                                                   <  5:68687/149‑1 (MQ=255)
                                                               cTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAg                                                                                   <  6:124495/149‑1 (MQ=255)
                                                                  cAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAgg                                                                                  <  4:196475/147‑1 (MQ=255)
                                                                            tgcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCtt                                                                      >  8:103439/1‑149 (MQ=255)
                                                                            tgcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCtt                                                                      >  8:267855/1‑149 (MQ=255)
                                                                             gcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGa                                                                                                                                      >  4:284853/1‑84 (MQ=255)
                                                                             gcAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGa                                                                                                                                      <  3:284853/84‑1 (MQ=255)
                                                                              cAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTgg                                                                    >  5:214859/1‑149 (MQ=255)
                                                                                    ggAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGTGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCgggg                                                               <  1:75356/148‑1 (MQ=255)
                                                                                     gAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtc                                                             >  2:68623/1‑149 (MQ=255)
                                                                                        tcgtcgTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                          >  8:313079/1‑149 (MQ=255)
                                                                                            cgtTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCgaga                                                      <  7:103439/149‑1 (MQ=255)
                                                                                                  gCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                                <  5:24851/149‑1 (MQ=255)
                                                                                                   cccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAggg                                                                                 >  8:160367/1‑115 (MQ=255)
                                                                                                   cccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAggg                                                                                 <  7:160367/115‑1 (MQ=255)
                                                                                                   cccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                                >  2:306498/1‑148 (MQ=255)
                                                                                                    ccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTa                                              <  7:267855/149‑1 (MQ=255)
                                                                                                       cgcAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTc                                           <  3:224211/149‑1 (MQ=255)
                                                                                                            gCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGa                                      >  3:75510/1‑149 (MQ=255)
                                                                                                              cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAg                                                                                                                                            >  7:311359/1‑45 (MQ=255)
                                                                                                              cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAg                                                                                                                                            <  8:311359/45‑1 (MQ=255)
                                                                                                              cAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATc                                    >  2:265460/1‑149 (MQ=255)
                                                                                                                   aagGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGcc                                                                        >  4:20517/1‑108 (MQ=255)
                                                                                                                      ggTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAgg                            <  6:305358/149‑1 (MQ=255)
                                                                                                                           gTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGccc                       >  6:244572/1‑149 (MQ=255)
                                                                                                                               gAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGGAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGa                                      <  8:112548/130‑1 (MQ=255)
                                                                                                                               gAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGa                                      >  7:112548/1‑130 (MQ=255)
                                                                                                                                aaGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                                <  6:48017/119‑1 (MQ=255)
                                                                                                                                    tgaggCCCGCTACCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcc              >  7:346415/4‑149 (MQ=255)
                                                                                                                                    aTCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcc              <  6:214859/149‑1 (MQ=255)
                                                                                                                                    aTCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcc              >  3:115061/1‑149 (MQ=255)
                                                                                                                                      aggCCCGCTACCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTcc              <  8:346415/146‑1 (MQ=255)
                                                                                                                                       ggCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCccg           >  7:156952/1‑149 (MQ=255)
                                                                                                                                            ggTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCccggccgg      <  5:244572/149‑1 (MQ=255)
                                                                                                                                             gTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCt                                         >  4:91585/1‑113 (MQ=255)
                                                                                                                                             gTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCt                                         <  3:91585/113‑1 (MQ=255)
                                                                                                                                                gCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGctct  >  7:44604/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGTTCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATCGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGCCGGCTCT  >  NC_002947/3788257‑3788549

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