Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 743,420 C→T 8.2% intergenic (+99/+175) PP_0635 → / ← PP_0636 group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947743,4200CT 8.2% 180.3 / 4.2 61intergenic (+99/+175)PP_0635/PP_0636group II intron‑encoding maturase/cold shock DNA‑binding domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base T (2/3);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGG  >  NC_002947/743272‑743562
                                                                                                                                                     |                                                                                                                                               
ggCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGc                                                                                                                                                 >  6:144139/1‑149 (MQ=255)
      tctcGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCTCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTa                                                                                                                                           <  7:188241/149‑1 (MQ=255)
          gATACTGGATACGCTTAACCGGCTTAGCCGCGTAACCTCAACCGCCGTGTACGGAACCGAACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTAtg                                                                                                                                         <  6:250568/147‑1 (MQ=255)
               tGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGc                                                                                                                                   >  6:12524/1‑148 (MQ=255)
                 gATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAAc                                                                                                                                >  4:151947/1‑149 (MQ=255)
                 gATACGATTAACCGGCTTAGCCGAGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCaa                                                                                                                                 >  7:51269/1‑148 (MQ=255)
                   tACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCt                                                                                                                              >  2:127528/1‑149 (MQ=255)
                      gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAGCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                           >  1:236930/1‑149 (MQ=255)
                      gATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCg                                                                                                                           >  5:309895/1‑149 (MQ=255)
                       aTTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAAccgaagg                                                                                                                          >  7:106462/1‑144 (MQ=255)
                             cGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGt                                                                                                                    >  1:120799/1‑149 (MQ=255)
                               gCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTa                                                                                                                 <  2:58513/149‑1 (MQ=255)
                                cTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTa                                                                                                                 >  2:204006/1‑148 (MQ=255)
                                   aGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACa                                                                                                               <  4:92867/147‑1 (MQ=255)
                                   aGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACAt                                                                                                              <  7:26473/148‑1 (MQ=255)
                                    gCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAg                                                                                                            >  6:276730/1‑149 (MQ=255)
                                     ccGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGc                                                                                                           <  7:132851/149‑1 (MQ=18)
                                     ccGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGc                                                                                                           <  8:106462/149‑1 (MQ=18)
                                            aCCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGacaca                                                                                                    >  6:15987/1‑149 (MQ=255)
                                              cTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGgcacac                                                                                                   >  4:247184/1‑148 (MQ=18)
                                               tGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAg                                                                                                 >  3:129229/1‑149 (MQ=255)
                                                 aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGTCCGCTATGTAACGCAACCGAAGGTTGCTTCACATAGCGGGcacacaa                                                                                                 <  3:247184/147‑2 (MQ=18)
                                                 aaCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGc                                                                                                <  2:236930/148‑1 (MQ=255)
                                                  aCCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgt                                                                                              >  3:261780/1‑149 (MQ=255)
                                                     gccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                           <  3:151947/149‑1 (MQ=255)
                                                     gccgTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtcgt                                                                                           <  5:12524/149‑1 (MQ=255)
                                                           tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtc                                                                                             <  2:211991/141‑1 (MQ=255)
                                                           tACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGcgtc                                                                                             >  1:211991/1‑141 (MQ=255)
                                                             cGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTAACATAGCGGACACACAGCGTCGTGACTTGc                                                                                    <  1:204006/149‑1 (MQ=255)
                                                              accaCCGTGCGTACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCg                                                                                  >  2:280500/4‑149 (MQ=255)
                                                              accaCCGTGCGTACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCg                                                                                  >  4:189171/4‑149 (MQ=255)
                                                              accaCCGTGCGTACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCg                                                                                  >  4:327533/4‑149 (MQ=255)
                                                                aaCCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCggg                                                                                <  8:51269/149‑1 (MQ=255)
                                                                  ccGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGa                                                                              <  1:127528/149‑1 (MQ=255)
                                                                    gttcGCACGGTGGTGTTAGAGGACGGAGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGaaa                                                                            <  4:61456/146‑1 (MQ=255)
                                                                    gTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGaaa                                                                            >  3:210263/1‑149 (MQ=255)
                                                                          aCGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACtt                                                                      >  7:26730/1‑149 (MQ=255)
                                                                           cGGTTGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTg                                                                     <  5:221558/149‑1 (MQ=255)
                                                                                   tgAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCgcggc                                                             <  6:128607/149‑1 (MQ=255)
                                                                                     agagGACGGCGGGTGTGAATCCGCCTACTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTAAAAAGCGGACACACAGCGTCGTGAATTGCCGGGCAAACTACTTGGTAGCGGCgt                                                           >  6:110194/1‑149 (MQ=39)
                                                                                        ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg                                                        <  2:120799/149‑1 (MQ=255)
                                                                                        ggACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGtgtg                                                        <  4:129229/149‑1 (MQ=255)
                                                                                           cggcggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAg                                                     <  5:144139/149‑1 (MQ=255)
                                                                                              cggGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAg                                                     >  4:262164/1‑146 (MQ=255)
                                                                                               gggTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCTTCTATGTGTTGTTTTCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                 <  1:280500/149‑1 (MQ=255)
                                                                                               gggTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                 <  3:327533/149‑1 (MQ=255)
                                                                                               gggTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCaga                                                 <  6:309895/149‑1 (MQ=255)
                                                                                                 gtgtGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTTGTTTCCCCGATATGTAAAGCAACATACGGTTGCTTAACTTACCGGAAACACACCTTCTTGACTTACAGGGAAAAATAATTTGTAGAGGCGtttgataca                                                   >  1:237442/1‑141 (MQ=255)
                                                                                                       aTCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAg                                         <  5:15987/149‑1 (MQ=255)
                                                                                                           gCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGa                                     <  5:276730/149‑1 (MQ=255)
                                                                                                                 tACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGtt                                >  5:237532/1‑148 (MQ=255)
                                                                                                                  aCTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                              >  1:225249/1‑149 (MQ=255)
                                                                                                                  aCTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                              >  1:323965/1‑149 (MQ=255)
                                                                                                                  aCTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGa                              >  1:329104/1‑149 (MQ=255)
                                                                                                                      gATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGt                          >  7:133752/1‑149 (MQ=255)
                                                                                                                        ttGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGTCTGTAGGGGAGCGGTGACAgagt                        >  8:310293/1‑146 (MQ=255)
                                                                                                                            gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCgt                                                           <  4:289728/110‑1 (MQ=255)
                                                                                                                            gcgAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGt                    >  5:308642/1‑149 (MQ=255)
                                                                                                                              gAAATCAGACCTGCCTTGTGTTCCCGCTATGTGAAGCAAACTTCGGTTGCGTTAAATAGCGGACACACAGCGTAGTGAATTGCAGGGGAAACTAATTGGTCGCGGCGTGTGGAGAAGATCATGTAGGGGAGCGTTGAAAGTTGCGGtt                   >  2:158652/1‑148 (MQ=255)
                                                                                                                                aaTCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCa                <  3:189171/149‑1 (MQ=255)
                                                                                                                                           ccTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTgg     <  8:133752/149‑1 (MQ=255)
                                                                                                                                               tgtgtgCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCgg  >  4:327128/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GGCTGGTCTCGATACTGGATACGATTAACCGGCTTAGCCGCGTAACCTGAACCGCCGTGTACGGAACCGTACGCACGGTGGTGTGAGAGGACGGCGGGTGTGAATCCGCCTCCTACTCGATTGCGCGAAATCAGACCTGCCTTGTGTGCCCGCTATGTGAAGCAACCTTCGGTTGCGTTACATAGCGGACACACAGCGTCGTGACTTGCCGGGGAAACTACTTGGTCGCGGCGTGTGGAGCAGAGACTGTAGGGGAGCGTTGACAGTGGCGGTTTCAAGCTAGCTTGGCGG  >  NC_002947/743272‑743562

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: