Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,149,444 C→T 8.5% intergenic (‑430/+18) PP_mr51 ← / ← PP_4535 C4 ncRNA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,149,4440CT 8.5% 196.4 / 6.2 70intergenic (‑430/+18)PP_mr51/PP_4535C4 ncRNA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (34/30);  new base T (2/4);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01

TGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTTCCTT  >  NC_002947/5149305‑5149590
                                                                                                                                           |                                                                                                                                                     
tGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTg                                                                                                                                               >  8:71827/1‑148 (MQ=255)
 gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGtt                                                                                                                                                <  4:239740/146‑1 (MQ=255)
 gAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGTTg                                                                                                                                               <  3:131906/147‑1 (MQ=255)
 gAGTCGCTGAGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGtt                                                                                                                                                <  4:120924/146‑1 (MQ=255)
  aGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGaa                                                                                                                                            <  7:329460/149‑1 (MQ=255)
  aGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGaa                                                                                                                                            <  4:347707/149‑1 (MQ=255)
   gTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACATCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAAc                                                                                                                                           <  7:71827/149‑1 (MQ=255)
         gtgCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGt                                                                                                                                     >  2:298019/1‑149 (MQ=255)
              ttCTGAATGACCACGAGCAATAATACTCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCGCATGCTACTTCCAACCAGTAGGCAATCCCCCCCTGGTTGGAACTTCAGTTACTc                                                                                                                                <  6:147301/149‑1 (MQ=255)
               tCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGTTg                                                                                                                                               <  7:216056/133‑1 (MQ=39)
               tCTGAATGACCACGAGCAATAAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGGGGGGATAGCCTCCTGGTTg                                                                                                                                               >  8:216056/1‑133 (MQ=38)
                     tGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAAt                                                                                                                          >  5:193586/1‑148 (MQ=255)
                       aCCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCa                                                                                                                        <  1:298019/148‑1 (MQ=255)
                        ccACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGa                                                                                                                      >  7:19415/1‑149 (MQ=255)
                         cACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGAc                                                                                                                     >  5:348388/1‑149 (MQ=255)
                                  tCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGtt                                                                                                            >  3:226484/1‑148 (MQ=255)
                                  tCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTtc                                                                                                           <  6:348388/149‑1 (MQ=255)
                                       cGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCAc                                                                                                      <  8:19415/149‑1 (MQ=255)
                                        gCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACcgcgc                                                                                                                <  6:53096/139‑1 (MQ=255)
                                        gCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACcgcgc                                                                                                                >  5:53096/1‑139 (MQ=255)
                                          ataataCATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATc                                                                                                                         >  7:223883/1‑128 (MQ=255)
                                          ataataCATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATc                                                                                                                         <  8:223883/128‑1 (MQ=255)
                                             ataCATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCCCCGAgcg                                                                                                >  1:330453/1‑149 (MQ=255)
                                              taCATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAgcgc                                                                                               >  1:241051/1‑149 (MQ=255)
                                                  tAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGAc                                                                                           <  4:226484/149‑1 (MQ=255)
                                                   aaTGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGAc                                                                                           >  8:121418/1‑148 (MQ=255)
                                                        gtCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGt                                                                                      >  5:178081/1‑148 (MQ=255)
                                                           cACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTc                                                                                  <  8:319773/149‑1 (MQ=255)
                                                           cACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTc                                                                                  >  5:244793/1‑149 (MQ=255)
                                                           cACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTc                                                                                  >  7:319773/1‑149 (MQ=255)
                                                            aCGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGAc                                                                                           <  3:338384/139‑1 (MQ=255)
                                                            aCGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGAc                                                                                           >  4:338384/1‑139 (MQ=255)
                                                              gccgccCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCg                                                                                 >  1:224207/1‑147 (MQ=255)
                                                              gccgccCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCg                                                                                 >  7:164120/1‑147 (MQ=255)
                                                              gccgccCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGaa                                                                              <  2:224207/149‑1 (MQ=255)
                                                              gccgccCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGaa                                                                              <  8:164120/149‑1 (MQ=255)
                                                                     gccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACa                                                                          <  5:328682/146‑1 (MQ=255)
                                                                     gccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACa                                                                          >  6:328682/1‑146 (MQ=255)
                                                                     gccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAg                                                                       <  2:241051/149‑1 (MQ=255)
                                                                     gccATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAg                                                                       <  6:193586/149‑1 (MQ=255)
                                                                        atatCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTg                                                                    <  6:178081/149‑1 (MQ=255)
                                                                           tCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTa                                                                                                                                   <  8:189248/85‑1 (MQ=255)
                                                                           tCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTa                                                                                                                                   >  7:189248/1‑85 (MQ=255)
                                                                            cTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCt                                                                >  2:223138/1‑149 (MQ=255)
                                                                             ttGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGACGACa                                                                          <  1:236428/138‑1 (MQ=255)
                                                                             ttGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTg                                                               >  1:236587/1‑149 (MQ=255)
                                                                             ttGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTg                                                               >  6:184152/1‑149 (MQ=255)
                                                                                       cTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAgg                                                    <  4:146376/149‑1 (MQ=255)
                                                                                        ttCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGc                                                   >  2:213147/1‑149 (MQ=255)
                                                                                        ttCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGc                                                   >  4:245800/1‑149 (MQ=255)
                                                                                                 cAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                          >  1:38184/1‑149 (MQ=255)
                                                                                                        aTCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCa                                   <  6:244793/149‑1 (MQ=255)
                                                                                                          ccTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                                 >  4:331866/1‑149 (MQ=255)
                                                                                                           cTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                                 <  5:184152/148‑1 (MQ=255)
                                                                                                           cTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                                 <  3:331866/148‑1 (MQ=255)
                                                                                                           cTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGAc                                <  7:121418/149‑1 (MQ=255)
                                                                                                               tCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCaa                                                            >  4:14296/1‑118 (MQ=255)
                                                                                                               tCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCaa                                                            <  3:14296/118‑1 (MQ=255)
                                                                                                                ccAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGa                           >  8:157271/1‑149 (MQ=255)
                                                                                                                 cAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACAGCAGGCCATCCGCGTGATCGCAGACAATGAt                          >  3:211878/1‑149 (MQ=255)
                                                                                                                  aaGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACCGGCAGGCCATCCGCGTGATCGCAGACAATGAt                          <  1:223138/149‑1 (MQ=255)
                                                                                                                     ttCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                                 >  2:252361/1‑138 (MQ=255)
                                                                                                                     gtCCAACCATGATGCTATCCCCCCCTTGTTGGTAATTCAGTTACTCGGCAATCAGACCGCGCGTTGTCCACCGTGCGCGGTCGTAGTCCTCGCAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGa                                 <  1:252361/138‑1 (MQ=255)
                                                                                                                         aaCCAGGGGGGGATAGCCTCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                          >  2:127497/1‑125 (MQ=38)
                                                                                                                         aaCCAGGGGGGGATAGCCTCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGt                                          <  1:127497/125‑1 (MQ=38)
                                                                                                                         aaCCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGc                   <  7:167756/148‑1 (MQ=255)
                                                                                                                         aaCCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGc                   >  8:167756/1‑148 (MQ=255)
                                                                                                                         aaCCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt                  >  6:355896/1‑149 (MQ=255)
                                                                                                                          aCCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTAGCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCt                  >  1:274345/1‑149 (MQ=255)
                                                                                                                                    tATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGtt       >  1:234168/1‑149 (MQ=255)
                                                                                                                                        ccccccTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTTCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGAGGCTTGAAGCCGGTTTCCt   >  4:226087/1‑149 (MQ=255)
                                                                                                                                        ccccccTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTTCCt   <  2:295459/149‑1 (MQ=255)
                                                                                                                                        ccccccTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCGTGACGCCGGTTTCCt   <  3:226087/149‑1 (MQ=255)
                                                                                                                                          ccccTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTTCCtt  >  5:131768/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                                     
TGAGTCGCTGTGCATTCTGAATGACCACGAGCAATCAGACGCATAATACATAATGTGTCCACGCCGCCCGCCATATCTTGGCTTCAACTTCAGTGTCCAAGCCAATCCTCATCCAAGTTCCAACCAGGAGGCTATCCCCCCCTGGTTGGAACTTCAGTTACTCGGCAATCAGACCGCGCGTTCTCCACCGAGCGCGGACGAAGTCCTCGAAGACATAGTTGCCCTGCAACGGCAGGCCATCCGCGTGATCGCAGACAATGATTTGCGGCTTGACGCCGGTTTCCTT  >  NC_002947/5149305‑5149590

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: