Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,081 G→C 6.2% A434G (GCC→GGC)  trmF ← tRNA uridine 5‑carboxymethylaminomethyl modification protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,0810GC 6.2% 201.3 / 3.9 64A434G (GCC→GGC) trmFtRNA uridine 5‑carboxymethylaminomethyl modification protein
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base C (2/2);  total (32/32)
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.49e-01

GATTTCAGGCGCTGTTCCTCACGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAG  >  NC_002947/2934‑3226
                                                                                                                                                   |                                                                                                                                                 
gATTTCAGGCGCTGTTCCTCACGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCAAGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCCCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCgg                                                                                                                                                   >  1:72446/1‑148 (MQ=255)
  tttCAGGCGCTGTTCCTCACGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCAc                                                                                                                                                >  5:199652/1‑149 (MQ=255)
      aGGCGCTGTTCCTCACGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGa                                                                                                                                            <  1:202503/149‑1 (MQ=255)
       ggCGCTGTTCCTCACGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAg                                                                                                                                           >  2:45303/1‑149 (MQ=255)
        gcgcTGTTCCTCACGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAgg                                                                                                                                          <  6:171073/149‑1 (MQ=255)
              ttCCTCACGCTCGATACCGTCGCGCTTTGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACa                                                                                                                                    <  6:276635/149‑1 (MQ=255)
                    aCGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCgg                                                                                                                               <  6:199652/148‑1 (MQ=255)
                    aCGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCg                                                                                                                                <  1:45303/147‑1 (MQ=255)
                      gCTCGATACCGTCGCGCTTGGCTCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTAc                                                                                                                            <  5:144429/149‑1 (MQ=255)
                      gCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTAc                                                                                                                            >  7:289013/1‑149 (MQ=255)
                      gCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTAc                                                                                                                            <  7:200870/149‑1 (MQ=255)
                        tCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACgg                                                                                                                          >  7:26382/1‑149 (MQ=255)
                           aTACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTc                                                                                                                       <  6:158207/149‑1 (MQ=255)
                            tACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCTGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTcc                                                                                                                      <  4:244567/149‑1 (MQ=255)
                            tACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTcc                                                                                                                      >  7:230669/1‑149 (MQ=255)
                                     gcTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTAccc                                                                                                             <  4:163210/149‑1 (MQ=255)
                                         ggCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAg                                                                                                          <  8:26382/148‑1 (MQ=255)
                                          gcgcAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAgg                                                                                                         <  6:97596/148‑1 (MQ=255)
                                                aaTGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAg                                                                                                  >  4:2087/1‑149 (MQ=255)
                                                    cGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGgtcg                                                                                              >  3:204706/1‑149 (MQ=255)
                                                     ggCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAgg                                                                                                         >  5:82095/1‑137 (MQ=255)
                                                     ggCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAgg                                                                                                         <  6:82095/137‑1 (MQ=255)
                                                      gCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCCCGCAGTATCAGCCGGTACTCTGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGgtcgt                                                                                             >  1:290595/1‑148 (MQ=255)
                                                        ccAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTcca                                                                                          <  8:289013/149‑1 (MQ=255)
                                                         cAGCTCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTccac                                                                                         <  6:241282/149‑1 (MQ=255)
                                                         cAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTccac                                                                                         >  5:34721/1‑149 (MQ=255)
                                                          aGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTccac                                                                                         >  3:148098/1‑148 (MQ=255)
                                                             gcTGGTCATCGATCAGGCCCAGTTAACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAg                                                                                      >  5:237358/1‑148 (MQ=255)
                                                                 gTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGc                                                                                 <  6:237358/149‑1 (MQ=255)
                                                                     tCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGAt                                                                             >  8:284172/1‑149 (MQ=255)
                                                                       gATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATg                                                                            >  8:155274/1‑148 (MQ=255)
                                                                        aTCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGt                                                                                                                              <  5:190783/97‑1 (MQ=255)
                                                                        aTCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGt                                                                                                                              >  6:190783/1‑97 (MQ=255)
                                                                             gCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGAt                                                                             >  8:286110/1‑141 (MQ=255)
                                                                             gCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGAt                                                                             <  7:286110/141‑1 (MQ=255)
                                                                             gCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCt                                                                     >  7:72439/1‑149 (MQ=255)
                                                                             gCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCt                                                                     >  2:111176/1‑149 (MQ=255)
                                                                                     cACGGGCCTTTTCGGTCAGGCCGAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCTGTTCTCGGCACTCGAGGTAAACTTGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcg                                                              <  8:265835/148‑1 (MQ=255)
                                                                                     cACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcg                                                              <  7:155274/148‑1 (MQ=255)
                                                                                     cACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCCGTACGTCTCCTGATTCCCCAAGGTAATAACGTCGACCACCAGTAAGCCGATGTATGCCTCGTcccg                                                              >  7:265835/1‑148 (MQ=255)
                                                                                      aCGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggc                                                            >  1:102619/1‑149 (MQ=255)
                                                                                      aCGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggc                                                            >  5:281208/1‑149 (MQ=255)
                                                                                        ggCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgc                                                          <  2:290595/149‑1 (MQ=255)
                                                                                          cccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcgg                                                        >  3:282374/1‑149 (MQ=255)
                                                                                          cccTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCACGTCGTCCACCAGTACGCCGATGTATGCCTCGTcgcggcgcgg                                                        >  8:89244/1‑149 (MQ=255)
                                                                                           cgttttCGGTCAGGCGCAGGTCGGCGTTGTCTTCTCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGAACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCggg                                                       <  2:72446/147‑1 (MQ=255)
                                                                                            cTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCggg                                                       >  3:311400/1‑148 (MQ=255)
                                                                                                  ggTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAg                                                 <  3:2087/148‑1 (MQ=255)
                                                                                                  ggTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGc                                                >  3:27039/1‑149 (MQ=255)
                                                                                                    tCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTg                                              <  8:230669/149‑1 (MQ=255)
                                                                                                          gcAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACgg                                        >  7:118008/1‑149 (MQ=255)
                                                                                                           cAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGCCCCGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGc                                       >  4:231238/1‑149 (MQ=255)
                                                                                                           cAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGCCACGCGAGGTAAACATGCGGTACGGATCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGc                                       >  4:226978/1‑149 (MQ=255)
                                                                                                            aGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGcc                                      <  8:72439/149‑1 (MQ=255)
                                                                                                                        tCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGACTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGaa                           <  3:35808/148‑1 (MQ=255)
                                                                                                                        tCTTCGCGAAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTTATCAGGTCTTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACAAg                                                 >  7:133260/1‑126 (MQ=255)
                                                                                                                        tCTACGCGCAGTATCAGCCGGTACTCGGCACGCGATGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAg                                                 <  8:133260/126‑1 (MQ=255)
                                                                                                                             gcgcAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTgcg                      <  6:281208/148‑1 (MQ=255)
                                                                                                                                    aTCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTa                >  7:244120/1‑147 (MQ=255)
                                                                                                                                    aTCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTa                <  8:244120/147‑1 (MQ=255)
                                                                                                                                       aGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTccc           <  4:204706/149‑1 (MQ=255)
                                                                                                                                             tACTCGCCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCcag     <  3:226978/149‑1 (MQ=255)
                                                                                                                                             tACTCGCCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCcag     <  3:231238/149‑1 (MQ=255)
                                                                                                                                              aCTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCcag     >  4:116058/1‑148 (MQ=255)
                                                                                                                                                tCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTAGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCcagcag  >  7:14220/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATTTCAGGCGCTGTTCCTCACGCTCGATACCGTCGCGCTTGGCGCAGAATGCGGCCCAGCGCTGGTCATCGATCAGGCCCAGTTCACGGCCCTTTTCGGTCAGGCGCAGGTCGGCGTTGTCTTCGCGCAGTATCAGCCGGTACTCGGCACGCGAGGTAAACATGCGGTACGGCTCCTGAGTACCCAAGGTAATCAGGTCGTCCACCAGTACGCCGATGTATGCCTCGTCGCGGCGCGGGCACCAGCTGTCACGGCCCTGCGCCCGAAGTGCGGCGTTAGTCCCGGCCAGCAG  >  NC_002947/2934‑3226

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: