Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,699,638 T→A 6.0% Q12L (CAG→CTG)  clsB ← cardiolipin synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,699,6380TA 6.0% 205.8 / 2.8 67Q12L (CAG→CTG) clsBcardiolipin synthase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base A (2/2);  total (32/35)
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.34e-01

GCAGCCTCGATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGATGCGGTCCAGGC  >  NC_002947/3699497‑3699784
                                                                                                                                             |                                                                                                                                                  
gCAGCCTCGATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCt                                                                                                                                             <  5:125249/149‑1 (MQ=255)
    ccTCGATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTg                                                                                                                                         >  4:72912/1‑149 (MQ=255)
    ccTCGATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTg                                                                                                                                         >  1:96602/1‑149 (MQ=255)
       cGATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCg                                                                                                                                      >  5:243297/1‑149 (MQ=255)
        gATCAACACCTGCTGCTATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGGAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACAGCGTGCAGTTGctga                                                                                                                                     <  6:78516/149‑4 (MQ=255)
        gATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGcc                                                                                                                                       <  2:157894/147‑1 (MQ=255)
        gATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGcc                                                                                                                                       <  2:162963/147‑1 (MQ=255)
           cAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCTTCTa                                                                                                                                  >  8:152301/1‑149 (MQ=255)
            aaCACCTGCTTCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACAGCGTGCAGTTGctgatcaa                                                                                                                                 <  5:111621/149‑8 (MQ=255)
             acacCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTAcc                                                                                                                                >  4:128930/1‑149 (MQ=255)
              cacCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTAcc                                                                                                                                >  8:121281/1‑148 (MQ=255)
                 ctgctgCAATTGCTTGCCCACTTTGTCGTCGAAGATTATGTAGGTTTCCAGCATTATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAg                                                                                                                            <  7:69030/149‑1 (MQ=255)
                  tgctgcAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAgg                                                                                                                           >  8:312667/1‑149 (MQ=255)
                   gctgcAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCACGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGt                                                                                                                          >  7:114764/1‑149 (MQ=255)
                       cAATTGCTTTCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCgg                                                                                                                       <  8:133582/148‑1 (MQ=255)
                        aaTTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGtt                                                                                                                     <  8:114764/149‑1 (MQ=255)
                        aaTTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTACACGCTGTTGCCGTCTACCCAAGGTCGGtt                                                                                                                     <  7:91727/149‑1 (MQ=255)
                             cttgcCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTcaca                                                                                                                 <  7:312667/148‑1 (MQ=255)
                             cttgcCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAg                                                                                                                <  7:315224/149‑1 (MQ=255)
                             cttgcCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAg                                                                                                                >  7:167745/1‑149 (MQ=255)
                             cttgcCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCAGTTCACAg                                                                                                                <  6:243297/149‑1 (MQ=255)
                                gcCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGTACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGt                                                                                                              <  1:226171/148‑1 (MQ=255)
                                gcCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGt                                                                                                              <  5:31213/148‑1 (MQ=255)
                                  ccACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGa                                                                                                           >  4:126727/1‑149 (MQ=255)
                                  ccACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGa                                                                                                           >  3:63686/1‑149 (MQ=255)
                                  ccACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTCGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTg                                                                                                            >  5:97934/1‑148 (MQ=255)
                                    aCTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGAtc                                                                                                         <  3:126727/149‑1 (MQ=255)
                                              cgAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGc                                                                                                 >  8:32313/1‑147 (MQ=255)
                                               gAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAg                                                                                              <  4:57574/149‑1 (MQ=255)
                                                     gatgaAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACGGGTTGATCTCCGCAGCCAGTGGTAc                                                                                        >  3:157210/1‑149 (MQ=255)
                                                       tgaAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGc                                                                                                    <  4:240036/135‑1 (MQ=255)
                                                       tgaAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGc                                                                                                    >  3:240036/1‑135 (MQ=255)
                                                       tgaAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACg                                                                                       <  3:223630/148‑1 (MQ=255)
                                                        gaAGGTTTCCAGCAATATTTCTTCCCGCGTCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTg                                                                                     >  2:188686/1‑149 (MQ=255)
                                                          aGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGt                                                                                   >  1:60924/1‑149 (MQ=255)
                                                             tttCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCg                                                                                                                        <  1:92174/109‑1 (MQ=255)
                                                             tttCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCg                                                                                                                        >  2:92174/1‑109 (MQ=255)
                                                                   gCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGAt                                                                          <  5:304626/149‑1 (MQ=255)
                                                                   gCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGAt                                                                          <  7:121281/149‑1 (MQ=255)
                                                                     aaTATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGTTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGc                                                                        <  6:97934/149‑1 (MQ=255)
                                                                     aaTATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATg                                                                         <  7:152301/148‑1 (MQ=255)
                                                                         tttCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTg                                                                                            >  5:181215/1‑125 (MQ=255)
                                                                         tttCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTg                                                                                            <  6:181215/125‑1 (MQ=255)
                                                                          ttcttcCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATg                                                                         <  3:130674/143‑1 (MQ=255)
                                                                          ttcttcCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATg                                                                         >  4:130674/1‑143 (MQ=255)
                                                                             ttcCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACgg                                                                >  5:97110/1‑149 (MQ=255)
                                                                             ttcCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGCTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACgg                                                                >  5:5891/1‑149 (MQ=255)
                                                                               cccGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACgg                                                                <  2:60924/147‑1 (MQ=255)
                                                                                      cTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGa                                                       <  7:32313/149‑1 (MQ=255)
                                                                                      cTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGAGCAGCAACTGCACGCGGTTGCCGTCTACCCAAGGTAGGTTCACAGGGTGATGTCCGCAGCCAGTGGTACGTGGTCGGACAGCTGCCACCCCGGGTATTTggg                                                       >  1:246472/1‑148 (MQ=255)
                                                                                       tGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACgg                                                                >  8:24080/1‑139 (MQ=255)
                                                                                       tGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACgg                                                                <  7:24080/139‑1 (MQ=255)
                                                                                          gCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTAGGACAGATGCGACCACGGGTATTTGGACAg                                                    >  5:244063/1‑148 (MQ=255)
                                                                                             aTGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCAcc                                                <  6:97110/149‑1 (MQ=255)
                                                                                                        cacGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGAGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGcc                                     >  5:249852/1‑149 (MQ=255)
                                                                                                           gcgGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCgg                                  >  8:118766/1‑149 (MQ=255)
                                                                                                                    aTACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGt                          <  2:96602/148‑1 (MQ=255)
                                                                                                                          cTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAg                   <  4:63686/149‑1 (MQ=255)
                                                                                                                             gCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTa                <  6:249852/149‑1 (MQ=255)
                                                                                                                                aTTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGAt             <  8:167745/149‑1 (MQ=255)
                                                                                                                                  tGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGATGc           >  5:224272/1‑149 (MQ=255)
                                                                                                                                   gATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACggg                                                               >  4:243917/1‑96 (MQ=255)
                                                                                                                                   gATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACggg                                                               <  3:243917/96‑1 (MQ=255)
                                                                                                                                    aTCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGATGCgg         >  1:253298/1‑149 (MQ=255)
                                                                                                                                      cggCAACAGCGTGCAGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGATGCGGTc       >  1:234402/3‑149 (MQ=255)
                                                                                                                                       ggCAACAGCGTGCAGTTGCTGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGATGCGGTc       >  7:171178/2‑148 (MQ=255)
                                                                                                                                           aCTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGATGCGGTCCAGgc  <  2:253298/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCAGCCTCGATCAACACCTGCTGCAATTGCTTGCCCACTTTGTCGTCGAAGATGATGAAGGTTTCCAGCAATATTTCTTCCCGCGCCTGGGCCATGGCCTCGAACACGCGGGGGTAATACTGCTCGCCATTGATCAGCAACTGCACGCTGTTGCCGTCTACCCAAGGTCGGTTCACAGGTTGATCTCCGCAGCCAGTGGTACGTGGTCGGACAGATGCGACCACGGGTATTTGGACAGCACCTGGGCTGCGCCCGGCACCGCGTTTCGCAGGTAGATGCGGTCCAGGC  >  NC_002947/3699497‑3699784

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑