Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,911,782:1 +T 6.6% coding (399/477 nt) ccmH ← cytochrome c heme lyase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,911,7821.T 6.6% 306.8 / 14.8 76coding (399/477 nt)ccmHcytochrome c heme lyase subunit
Reads supporting (aligned to +/- strand):  ref base . (37/34);  new base T (3/2);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

AGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAC  >  NC_002947/4911650‑4911925
                                                                                                                                     |                                                                                                                                               
agcagcAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGgc                                                                                                                                 <  1:354477/149‑1 (MQ=255)
     caAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAANCNTGGGTCTGTTCNTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCgg                                                                                                                            >  5:74455/1‑149 (MQ=255)
       aaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGcc                                                                                                                                           >  7:312054/1‑132 (MQ=255)
       aaaaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGcc                                                                                                                                           <  8:312054/132‑1 (MQ=255)
         aaGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTa                                                                                                                        <  5:333305/149‑1 (MQ=255)
           gCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCTCCCTGCACGGCCGTAGCACGGAGCCGGCGTAc                                                                                                                       >  2:30265/1‑148 (MQ=255)
            cTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGCCga                                                                                                                            <  6:213323/142‑2 (MQ=255)
            cTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCGTGCTACGGCCGTGCAGGGCGCCga                                                                                                                            >  5:213323/1‑141 (MQ=255)
            cTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGa                                                                                                                     >  6:273672/1‑149 (MQ=255)
                 ggCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCAcc                                                                                                                >  3:356913/1‑149 (MQ=255)
                       gAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGc                                                                                                          >  7:140410/1‑149 (MQ=255)
                            gcagAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGcaccac                                                                                                     >  4:351476/1‑149 (MQ=255)
                               gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACg                                                                                                    >  4:173003/1‑147 (MQ=255)
                               gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  8:140410/149‑1 (MQ=255)
                               gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  3:351476/149‑1 (MQ=255)
                               gAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  5:273672/149‑1 (MQ=255)
                                    cccGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCCCCGCCAGCACCACGAAGCCaa                                                                                             >  8:92091/1‑148 (MQ=255)
                                    cccGCGCTAAGCCAGAATTCAATCATAGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  <  2:52569/144‑1 (MQ=255)
                                    cccGCGCTAAGCCAGAATTCAATCATAGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGaa                                                                                                  >  1:52569/1‑144 (MQ=255)
                                     ccGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCAcc                                                                                            >  6:388291/1‑149 (MQ=255)
                                       gcgcTAAGCCTGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCAGGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCgg                                                                                          <  1:30265/149‑1 (MQ=255)
                                              gCCACAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCaa                                                                                   >  6:77857/1‑149 (MQ=255)
                                                 aGAATTCAATCATGGGGCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGAt                                                                                <  5:387797/149‑1 (MQ=255)
                                                    aTTCAATCATGGGTATGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATccc                                                                             <  5:309364/149‑1 (MQ=255)
                                                     ttCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGTACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACAAAGCCACCGGCCAGCAAGATcccc                                                                            >  3:220526/1‑149 (MQ=255)
                                                     ttCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCa                                                                                       >  2:22632/1‑138 (MQ=255)
                                                     ttCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCa                                                                                       <  1:22632/138‑1 (MQ=255)
                                                        aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCgg                                                                          <  5:388291/148‑1 (MQ=255)
                                                        aaTCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTAAGATCACCGCCAGCACCACCAAGCCACCGGCCAGCAAGATCCCCggc                                                                         >  5:137488/1‑148 (MQ=255)
                                                          tCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                        >  5:403832/1‑148 (MQ=255)
                                                          tCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGc                                                                        >  1:43339/1‑148 (MQ=255)
                                                                     ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACcaca                                                               <  6:403832/146‑1 (MQ=255)
                                                                     ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagca                                                            >  1:203464/1‑149 (MQ=255)
                                                                     ttCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcag                                                              >  1:138235/1‑147 (MQ=255)
                                                                       cTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACagcagc                                                          >  7:312664/1‑149 (MQ=255)
                                                                        tttttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACagcagc                                                          >  1:221580/1‑148 (MQ=255)
                                                                          tttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagc                                                             >  2:199008/1‑143 (MQ=255)
                                                                          tttCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCAcagc                                                             <  1:199008/143‑1 (MQ=255)
                                                                           ttCCATCAGTTTGGAGAGACGNTCGCGTTATTCGTCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAATCCACCGGCCATCAAGATCCNCGGGCCAAACCACAGCAGCCAg                                                       <  6:74455/148‑1 (MQ=255)
                                                                            tCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAgg                                                      >  4:163674/1‑148 (MQ=255)
                                                                             ccAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                           <  2:55803/126‑1 (MQ=255)
                                                                             ccAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                           >  1:55803/1‑126 (MQ=255)
                                                                                   gTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATcccgg                                                                           >  2:164420/1‑118 (MQ=255)
                                                                                   gTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCg                                                                           <  1:164420/120‑1 (MQ=255)
                                                                                    tttGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGc                                             >  6:239579/1‑149 (MQ=255)
                                                                                      tGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                            <  5:284198/148‑1 (MQ=255)
                                                                                         cGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAgc                                        <  5:77857/149‑1 (MQ=255)
                                                                                            gaCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGTACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACAAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGcc                                     <  4:220526/149‑1 (MQ=255)
                                                                                             aCGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCg                                    <  8:366769/149‑1 (MQ=255)
                                                                                              cGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCTAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCgg                                   <  1:268711/149‑1 (MQ=255)
                                                                                              cGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCTAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCgg                                   <  1:274632/149‑1 (MQ=255)
                                                                                                 tCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCtt                                >  6:364316/1‑149 (MQ=255)
                                                                                                   gcgTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTg                               <  2:43339/148‑1 (MQ=255)
                                                                                                   gcgTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACACCAGCCAGGTACGCCCGCTGAGCGCCGGCTTg                               >  6:161696/1‑148 (MQ=255)
                                                                                                    cgTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTa                             <  2:221580/149‑1 (MQ=255)
                                                                                                     gTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAg                            >  7:345524/1‑149 (MQ=255)
                                                                                                     gTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAg                            <  5:364316/149‑1 (MQ=255)
                                                                                                      ttcttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGAACACCGCCAGCACAACGCAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCCGCTTGTAgc                           >  6:356548/1‑149 (MQ=255)
                                                                                                        cttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAgcg                          >  5:377002/1‑148 (MQ=255)
                                                                                                         ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                            >  8:304908/1‑129 (MQ=255)
                                                                                                         ttcGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCt                                            <  7:304908/129‑1 (MQ=255)
                                                                                                           cGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCAc                       <  5:402550/148‑1 (MQ=255)
                                                                                                            ggCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCAc                       <  5:381876/147‑1 (MQ=255)
                                                                                                            ggCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCAc                       >  6:381876/1‑147 (MQ=255)
                                                                                                                gACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTc                 <  3:173003/149‑1 (MQ=255)
                                                                                                                 aCAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCg                >  6:29708/1‑149 (MQ=255)
                                                                                                                    aCTCATCGGCGCCCTCCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGACCTCGCCa             <  5:127128/149‑1 (MQ=255)
                                                                                                                       cATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCtt                                >  2:76441/1‑127 (MQ=255)
                                                                                                                       cATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCtt                                <  1:76441/127‑1 (MQ=255)
                                                                                                                        gccGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  <  1:309593/122‑1 (MQ=38)
                                                                                                                        gccGGCGCCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGc                                  >  2:309593/3‑124 (MQ=38)
                                                                                                                          cGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCgg       >  2:217531/1‑149 (MQ=255)
                                                                                                                           ggCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGt      <  5:29708/149‑1 (MQ=255)
                                                                                                                            gcgcCGTGCTACGGCCGTGCAGGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTc     >  3:378683/1‑149 (MQ=255)
                                                                                                                            gcgcCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTc     >  6:177087/1‑149 (MQ=255)
                                                                                                                               cccTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAc  <  5:239579/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                               
AGCAGCAAAAAGCTGAGGGCAGCGAGCAGCAGAAGGCCCGCGCTAAGCCAGAATTCAATCATGGGTCTGTTCTTTTTCCAGCAGTTTGGCGAGACGCTCGCGTTCTTCGGCAGACAACTCATCGGCGCCCTGCACGGCCGTAGCACGGCGCCGGCGTACGATCACCGCCAGCACCACGAAGCCACCGGCCAGCAAGATCCCCGGGCCAAACCACAGCAGCCAGGTACGCCCGCTGAGCGCCGGCTTGTAGCGCACGAAGTCGCCATAGCGGTCGAC  >  NC_002947/4911650‑4911925

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: