Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,592,302 T→A 6.5% G69G (GGT→GGA PP_1396 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,592,3020TA 6.5% 181.2 / 3.0 61G69G (GGT→GGAPP_1396hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/29);  new base A (2/2);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

GCAGAGCTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGTTCAACAACCCGGACAGCAACATGGCACGTAG  >  NC_002947/1592156‑1592449
                                                                                                                                                   |                                                                                                                                                    
gCAGAGCTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTa                                                                                                                                                    <  5:186898/149‑1 (MQ=255)
  agagCTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGtt                                                                                                                                                     >  4:319818/1‑146 (MQ=255)
  agagCTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGtt                                                                                                                                                     <  3:319818/146‑1 (MQ=255)
  agagCTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACt                                                                                                                                                  >  8:128012/1‑149 (MQ=255)
      cTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACtt                                                                                                                                                 >  4:15655/1‑146 (MQ=255)
      cTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACtt                                                                                                                                                 <  3:15655/146‑1 (MQ=255)
                    gCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTc                                                                                                                                           <  5:198353/138‑1 (MQ=255)
                    gCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTc                                                                                                                                           >  6:198353/1‑138 (MQ=255)
                     cGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTg                                                                                                                                <  7:128012/147‑1 (MQ=255)
                      gTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAg                                                                                                                             <  1:130534/149‑1 (MQ=255)
                      gTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAg                                                                                                                             <  1:72821/149‑1 (MQ=255)
                        gACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAg                                                                                                                             <  5:255181/147‑1 (MQ=255)
                            tGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGt                                                                                                                                                      >  4:276757/1‑119 (MQ=255)
                            tGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGt                                                                                                                                                      <  3:276757/119‑1 (MQ=255)
                            tGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTg                                                                                                                                               <  8:153268/126‑1 (MQ=255)
                            tGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTg                                                                                                                                               >  7:153268/1‑126 (MQ=255)
                              cTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCg                                                                                                                      <  6:83996/148‑1 (MQ=255)
                                   gCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCg                                                                                                                      <  7:123170/143‑1 (MQ=255)
                                   gCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCg                                                                                                                      >  8:123170/1‑143 (MQ=255)
                                      cTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCATCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAg                                                                                                                             <  4:51884/133‑1 (MQ=255)
                                      cTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGcc                                                                                                                                        <  6:238836/123‑1 (MQ=255)
                                      cTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAg                                                                                                                             >  3:51884/1‑133 (MQ=255)
                                       tCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGt                                                                                                            >  2:270675/1‑149 (MQ=255)
                                                 cacCGGGTTGCCTGAGTTTTACCGTGCCGCGTTCTCAGGCGGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAg                                                                                                  <  5:330873/149‑1 (MQ=255)
                                                   ccGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCGTAAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAgg                                                                                                 >  1:274134/1‑149 (MQ=255)
                                                           ccTGAGTTTTAGCGTGCAGCGTTCTCTGGCTGATGCTGACCTGTGGTTGCTCAGCGGATGCTGGCAGGACCAGCATGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTcagc                                                                                         <  5:73481/148‑1 (MQ=255)
                                                           ccTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGTAGCTGGCAGGACCAGCAGGCGATGAGTGGTCACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGc                                                                                                      <  6:31282/135‑1 (MQ=255)
                                                           ccTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGTAGCTGGCAGGACCAGCAGGCGATGAGTGGTCACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGc                                                                                                      >  5:31282/1‑135 (MQ=255)
                                                           ccTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTcag                                                                                          >  2:236839/1‑147 (MQ=255)
                                                             tGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCagcagc                                                                                      >  2:327729/1‑149 (MQ=255)
                                                                        gTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACCTCTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGc                                                                            <  1:270675/149‑1 (MQ=255)
                                                                        gTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGATACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCa                                                                           >  5:7771/1‑149 (MQ=255)
                                                                        gTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGATACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCa                                                                           >  5:2778/1‑149 (MQ=255)
                                                                                  ctcAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAACCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGAc                                                                 >  7:10292/1‑149 (MQ=255)
                                                                                       gCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTg                                                               >  1:312547/1‑146 (MQ=255)
                                                                                       gCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTg                                                               <  2:312547/146‑1 (MQ=255)
                                                                                         tGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGt                                                          >  8:190990/1‑149 (MQ=255)
                                                                                           aTGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCTATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTg                                                         <  3:40924/148‑1 (MQ=255)
                                                                                           aTGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTg                                                         <  1:327729/148‑1 (MQ=255)
                                                                                             gCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAg                                                                                                  >  7:153136/1‑105 (MQ=255)
                                                                                              cTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGAt                                                      >  1:191680/1‑148 (MQ=255)
                                                                                                gACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTg                                                                             >  1:339208/1‑123 (MQ=255)
                                                                                                gACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTg                                                                             <  2:339208/123‑1 (MQ=255)
                                                                                                    tgtgGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTa                                                >  3:282057/1‑148 (MQ=255)
                                                                                                      tgGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAg                                             <  7:190990/149‑1 (MQ=255)
                                                                                                      tgGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGATACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTaaa                                              <  6:7771/148‑1 (MQ=255)
                                                                                                      tgGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGATACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTaaa                                              <  6:2778/148‑1 (MQ=255)
                                                                                                              tCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTgg                                                        >  3:47262/1‑131 (MQ=255)
                                                                                                              tCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTgg                                                        <  4:47262/131‑1 (MQ=255)
                                                                                                              tCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGct                                      >  6:310228/1‑149 (MQ=255)
                                                                                                                aGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAg                                                                                                                             >  6:172275/1‑60 (MQ=255)
                                                                                                                aGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACg                                        >  5:286503/1‑145 (MQ=255)
                                                                                                                aGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACg                                        <  6:286503/145‑1 (MQ=255)
                                                                                                                   ggAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGt                                 <  1:236839/149‑1 (MQ=255)
                                                                                                                    gAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGtt                                >  7:174204/1‑149 (MQ=255)
                                                                                                                     aaGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCATCGACTGAGCGTGGATGAGGTAAAGt                                            >  2:217208/1‑136 (MQ=255)
                                                                                                                     aaGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCATCGACTGAGCGTGGATGAGGTAAAGt                                            <  1:217208/136‑1 (MQ=255)
                                                                                                                             aGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACt                                                                <  4:129216/108‑1 (MQ=255)
                                                                                                                             aGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACt                                                                >  3:129216/1‑108 (MQ=255)
                                                                                                                               gACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGTTCAACAACCCg                      >  1:49600/1‑148 (MQ=255)
                                                                                                                                 ccAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGTTCAACAACCCGGa                    >  3:22036/1‑148 (MQ=255)
                                                                                                                                   agcagGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGACGAGGTAAAGTGACGCTCTTGTTCAACAAGCCGGTCAg                 <  8:10292/149‑1 (MQ=255)
                                                                                                                                    gcagGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAgg                                                  >  6:345928/1‑115 (MQ=255)
                                                                                                                                    gcagGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAgg                                                  <  5:345928/115‑1 (MQ=255)
                                                                                                                                                  gTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGTTCAACAACCCGGACAGCAACATGGCACGTAg  <  8:174204/149‑1 (MQ=255)
                                                                                                                                                   ttACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGTTCAACAACCCGGACAGCAACATg          <  1:54466/140‑1 (MQ=255)
                                                                                                                                                   ttACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGTTCAACAACCCGGACAGCAACATg          >  2:54466/1‑140 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGAGCTGGGCGCGCGCCTGCGTGACCTGCTCGAGCCCTCGCTGCGTACACCGGGTTGCCTGAGTTTTACCGTGCAGCGTTCTCAGGCTGATGCTGACCTGTGGTTGCTCAGCGGAAGCTGGCAGGACCAGCAGGCGATGAGTGGTTACTTTGCTTCGCCGACACTGGAGGTGTTCGGTGAACTGGTGCAGGCCCAGGTGGTCAGCAGCCTGGACTTGCATACCTTCGACTGAGCGTGGATGAGGTAAAGTGACGCTCTTGTTCAACAACCCGGACAGCAACATGGCACGTAG  >  NC_002947/1592156‑1592449

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: