Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,285,434 A→T 6.5% V144D (GTC→GAC)  PP_2887 ← catalase‑like peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,285,4340AT 6.5% 180.5 / 2.4 61V144D (GTC→GAC) PP_2887catalase‑like peroxidase
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base T (2/2);  total (33/28)
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.63e-01

TAACCCGCCGCCCATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGCAAAGCGGCCAACCACCGGCA  >  NC_002947/3285293‑3285581
                                                                                                                                             |                                                                                                                                                    
tAACCCGCCGCCCATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCa                                                                                                                                                <  7:314000/148‑1 (MQ=255)
  aCCCGCCGCCCATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCa                                                                                                                                                >  5:93830/1‑146 (MQ=255)
  aCCCGCCGCCCATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCa                                                                                                                                                <  6:93830/146‑1 (MQ=255)
    ccgccgCCCATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTccc                                                                                                                                           >  2:344703/1‑149 (MQ=255)
        cgccCATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACcggcg                                                                                                                                       >  6:114470/1‑149 (MQ=255)
            cATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGgtt                                                                                                                                  >  2:395069/1‑149 (MQ=255)
             aTTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGgttg                                                                                                                                 >  8:294366/1‑149 (MQ=255)
              ttGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGgttgt                                                                                                                                <  5:114470/149‑1 (MQ=255)
              ttGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGgttgt                                                                                                                                <  3:299386/149‑1 (MQ=255)
                     aaCGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGcc                                                                                                                         >  8:172215/1‑149 (MQ=255)
                     aaCGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGcc                                                                                                                         >  8:177464/1‑149 (MQ=255)
                         ccgcGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTc                                                                                                                     >  4:383947/1‑149 (MQ=255)
                          cgcgCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTcc                                                                                                                    <  6:46972/149‑1 (MQ=255)
                          cgcgCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTcc                                                                                                                    >  4:289268/1‑149 (MQ=255)
                              cTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGTCTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGc                                                                                                                  >  6:84752/1‑147 (MQ=255)
                                     ggTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTgg                                                                                                         <  7:177464/149‑1 (MQ=255)
                                     ggTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTgg                                                                                                         <  7:172215/149‑1 (MQ=255)
                                          gcggcgAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCAtc                                                                                                    >  3:358953/1‑149 (MQ=255)
                                              cgAAGACCGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGGCGGGTCTGGTGCACTCGCCAGGACTTGGTGGTAGACCGGGTGGGGCGTGCTGACTGCCAGTACCGGCGGGGTGTTCATGCCGGTCAGCCAGACCTGGCCAtcggcg                                                                                                >  2:305568/1‑149 (MQ=255)
                                               gaagaaCGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCgg                                                                                               >  2:57311/1‑149 (MQ=255)
                                               gaagaaCGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCgg                                                                                               <  1:228380/149‑1 (MQ=255)
                                                aagaaCGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGt                                                                                              >  5:301964/1‑149 (MQ=255)
                                                aagaaCGCCTGCAGCTTGCCAGGGTCCGGTCTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGt                                                                                              >  7:303378/1‑149 (MQ=255)
                                                   aaCGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGTCTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTgct                                                                                           >  7:219179/1‑149 (MQ=255)
                                                     cGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCGGTACTTGCTGGTAGAATGCTTGGGGCGGGCTGCCTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCGATCGGCGGGgctcc                                                                                         >  4:121107/1‑148 (MQ=255)
                                                     cGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTgct                                                                                           >  1:383921/1‑147 (MQ=255)
                                                                tGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGc                                                                               >  5:116351/1‑148 (MQ=255)
                                                                 gCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCGTGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGc                                                                               >  1:219428/1‑147 (MQ=255)
                                                                 gCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGcc                                                                              <  4:358953/148‑1 (MQ=255)
                                                                   cAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAgg                                                                           <  1:57311/149‑1 (MQ=255)
                                                                       gTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGc                                                                       <  4:274562/149‑1 (MQ=255)
                                                                           ggTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTAc                                                                   >  7:50517/1‑149 (MQ=255)
                                                                               tGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCgg                                                                <  1:395069/148‑1 (MQ=255)
                                                                                   ggTCGCCGGGTCTGGTTCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATAc                                                           <  3:383947/149‑1 (MQ=255)
                                                                                     tCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATAccc                                                         <  6:157591/149‑1 (MQ=255)
                                                                                     tCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATAccc                                                         <  6:355951/149‑1 (MQ=255)
                                                                                        ccGGGTCTAGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTa                                                      >  2:361407/1‑149 (MQ=255)
                                                                                            gTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCg                                                  <  7:110686/149‑1 (MQ=255)
                                                                                                 gTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTgctcagctc                                                                                     <  5:26261/109‑1 (MQ=255)
                                                                                                 gTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTgctcagctc                                                                                     >  6:26261/1‑109 (MQ=255)
                                                                                                 gTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCgggg                                               <  6:116351/147‑1 (MQ=255)
                                                                                                 gTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGc                                              >  3:311000/1‑148 (MQ=255)
                                                                                                 gTGCACTCGCCAGTACTTGCTGGTAGAACGCGTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCggg                                                <  2:391045/146‑1 (MQ=255)
                                                                                                    cACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAAc                                          <  3:285101/149‑1 (MQ=255)
                                                                                                    cACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAAc                                          >  4:285101/1‑149 (MQ=255)
                                                                                                       tCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACggg                                       >  8:34295/1‑149 (MQ=255)
                                                                                                       tCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACggg                                       >  3:197232/1‑149 (MQ=255)
                                                                                                         gCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGAATTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGt                                      >  1:101579/1‑149 (MQ=255)
                                                                                                         gCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGGGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACggg                                       <  3:289268/147‑1 (MQ=255)
                                                                                                          ccAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTg                                    <  5:64553/149‑1 (MQ=255)
                                                                                                           cAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGTCTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTgg                                   <  5:84752/149‑1 (MQ=255)
                                                                                                              tACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCAc                                <  6:301964/149‑1 (MQ=255)
                                                                                                                 ttGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTgg                                   <  7:45161/143‑1 (MQ=255)
                                                                                                                 ttGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTgg                                   >  8:45161/1‑143 (MQ=255)
                                                                                                                        tAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATgg                       >  1:296022/1‑148 (MQ=255)
                                                                                                                         aGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGc                      >  7:327276/1‑148 (MQ=255)
                                                                                                                           aaCGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGCaaa                   <  7:294366/149‑1 (MQ=255)
                                                                                                                              gCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGCAAAGCg                >  6:327516/1‑149 (MQ=255)
                                                                                                                               cTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGCAAAGCgg               <  1:361407/149‑1 (MQ=255)
                                                                                                                                ttGGGGCGTGCTGTCTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGCAAAGCGGc              <  8:219179/149‑1 (MQ=255)
                                                                                                                                            gACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGCAAAGCGGCCAACCCCCGGcc  >  8:263602/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TAACCCGCCGCCCATTGGCGAAACGCCGCGCTTTCCGGGTGGGCGGCGAAGAACGCCTGCAGCTTGCCAGGGTCCGGTTTGCCGGTCGCCGGGTCTGGTGCACTCGCCAGTACTTGCTGGTAGAACGCTTGGGGCGTGCTGACTGCCAGTACCGGCGGGTTGTTCATGCCGGTCCGCCAGACCTGGCCATCGTCGGTGCTCAGCTCAATCGCCAGGCTGCGTACCGGTATACCCGTATCCGGGGCGAACGGGTTGGCACCGCCAATGGCAAAGCGGCCAACCACCGGCA  >  NC_002947/3285293‑3285581

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: