Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,732,541 G→C 6.2% A89G (GCC→GGC)  sucB ← 2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,732,5410GC 6.2% 203.4 / 4.2 65A89G (GCC→GGC) sucB2‑oxoglutarate dehydrogenase dihydrolipoyltranssuccinylase subunit
Reads supporting (aligned to +/- strand):  ref base G (30/31);  new base C (2/2);  total (32/33)
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.68e-01

GTGACGCGACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAACCTTGT  >  NC_002947/4732397‑4732682
                                                                                                                                                |                                                                                                                                             
gTGACGCGACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGcc                                                                                                                                          >  2:1410880/1‑150 (MQ=255)
     gcgACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCCCGGc                                                                                                                                     >  2:1660930/1‑150 (MQ=255)
        aCCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCGCCGGCAGCAGCGGCCGGAGCAGCTGCAGGCGCGCCTGCCGCGGCAgc                                                                                                                                  >  2:1707561/1‑149 (MQ=255)
          cGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGgcg                                                                                                                                >  2:7634/1‑150 (MQ=255)
            ccTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCGCCGGCGGGCGCa                                                                                                                              >  2:542259/1‑150 (MQ=255)
              tttACCGGTACGTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGCCTGCCGCGGCAGCCGCacc                                                                                                                            <  1:1707561/150‑3 (MQ=255)
                 aCCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCACGGCAGGCGCAGCGgc                                                                                                                         <  1:650889/150‑1 (MQ=255)
                  ccGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGgca                                                                                                                        <  1:1390819/150‑1 (MQ=255)
                    ggTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGgcagc                                                                                                                      >  2:1730599/1‑150 (MQ=255)
                      tACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCAc                                                                                                                    <  1:947809/150‑1 (MQ=255)
                           gCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCCGGCACAGAGGCAGCACCGccc                                                                                                               <  1:1500540/150‑1 (MQ=255)
                              aCGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCGGCACCGcc                                                                                                                >  2:177247/1‑146 (MQ=255)
                              aCGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCGGCACCGcc                                                                                                                <  1:177247/146‑1 (MQ=255)
                                    gCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATc                                                                                                      >  1:122893/1‑150 (MQ=255)
                                    gCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATc                                                                                                      <  2:1097594/150‑1 (MQ=255)
                                     ccAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCg                                                                                                     >  1:718671/1‑150 (MQ=255)
                                      cAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGa                                                                                                    >  1:627435/1‑150 (MQ=255)
                                          tCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAAccc                                                                                                >  2:1007844/1‑150 (MQ=255)
                                            gATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAg                                                                                              >  2:185028/1‑150 (MQ=255)
                                             aTGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGc                                                                                             >  2:435200/1‑150 (MQ=255)
                                                      tcttcAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATtcgtcg                                                                                    <  2:984238/150‑1 (MQ=255)
                                                      tcttcAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATtcgtcg                                                                                    <  2:1631485/150‑1 (MQ=255)
                                                      tctgcAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATtcgtcg                                                                                    <  1:767971/150‑1 (MQ=255)
                                                        ttcAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCATCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGa                                                                                  >  1:1429390/1‑150 (MQ=255)
                                                        ttcAGCCAGCTTGCGCACGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGa                                                                                  >  2:815839/1‑150 (MQ=255)
                                                          cagccagcTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACa                                                                                >  1:1503981/1‑150 (MQ=255)
                                                             ccagcTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGa                                                                             >  2:204792/1‑150 (MQ=255)
                                                               agcTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACg                                                                           >  2:1238638/1‑150 (MQ=255)
                                                                  ttGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGgtg                                                                        <  1:853281/150‑1 (MQ=255)
                                                                      gcgcgGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGtcgc                                                                    >  2:106002/1‑150 (MQ=255)
                                                                        gcgGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGtcgccc                                                                  <  2:122893/150‑1 (MQ=255)
                                                                             cGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGtcgccct                                                                 >  2:695469/1‑146 (MQ=255)
                                                                             cGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGtcgccct                                                                 <  1:695469/146‑1 (MQ=255)
                                                                              ggggCTGCGATTGGGTCATCTTCGCGGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCCCCCTCAACGCTCGACCCCAGCAACTCGTCGGACAGGTCGGTGtcgccctcgccc                                                            <  2:7724/150‑1 (MQ=255)
                                                                              ggggCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCACAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGtcgccctcgccc                                                            <  2:1242900/150‑1 (MQ=255)
                                                                                   tGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGAc                                                       >  2:236580/1‑150 (MQ=255)
                                                                                     cGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGa                                                     <  1:542259/150‑1 (MQ=255)
                                                                                      gATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGAt                                                    <  1:7634/150‑1 (MQ=255)
                                                                                           ggTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGc                                               >  2:1230137/1‑150 (MQ=255)
                                                                                             tCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCGAGCTATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGccc                                             >  1:610896/1‑150 (MQ=255)
                                                                                             tCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGccc                                             <  1:266840/150‑1 (MQ=255)
                                                                                                 cTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAAccc                                                                                                <  2:1746892/95‑1 (MQ=255)
                                                                                                 cTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAAccc                                                                                                >  1:1746892/1‑95 (MQ=255)
                                                                                                 cTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCccacca                                         >  2:386665/1‑150 (MQ=255)
                                                                                                      ccGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCa                                    >  1:1143348/1‑150 (MQ=255)
                                                                                                       cGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCAt                                   >  1:555612/1‑150 (MQ=255)
                                                                                                            tCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCa                              <  2:1202564/150‑1 (MQ=255)
                                                                                                               gCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTa                           <  1:204792/150‑1 (MQ=255)
                                                                                                               gCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTa                           <  2:555612/150‑1 (MQ=255)
                                                                                                                     gcagcGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGt                     >  2:1005216/1‑150 (MQ=255)
                                                                                                                       agcGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTAc                   <  1:1730599/150‑1 (MQ=255)
                                                                                                                         cGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACtt                 <  1:1384036/150‑1 (MQ=255)
                                                                                                                            ccGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCGCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCa              <  1:1692795/150‑1 (MQ=255)
                                                                                                                             cGGAGCATCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGCACTTCCAg             <  2:1429390/150‑1 (MQ=255)
                                                                                                                             cGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAg             <  1:1238638/150‑1 (MQ=255)
                                                                                                                             cGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAg             <  1:1230137/150‑1 (MQ=255)
                                                                                                                             cGGAGAAGCTGCAGGTGCGGATGCCGCCGGAGGCGCAGCGGCAGCAGCTCTCTCAACGATCTAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAg             <  2:610896/150‑1 (MQ=255)
                                                                                                                               gAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGAc                                                       <  2:16619/106‑1 (MQ=255)
                                                                                                                               gAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGAc                                                       >  1:16619/1‑106 (MQ=255)
                                                                                                                                    gcTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAg                                                                               >  1:195402/1‑77 (MQ=255)
                                                                                                                                    gcTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAg                                                                               <  2:195402/77‑1 (MQ=255)
                                                                                                                                    gcTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGCCAGCGCCGCCATTAACGTTCGAACCCAGCAATTCGTAGGCATTTACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAAcc      <  1:435200/150‑1 (MQ=255)
                                                                                                                                        cAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAACCTTGt  >  2:623408/1‑150 (MQ=255)
                                                                                                                                          gctGCGCCTGCCGCGGCAGCCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGtcgccctcgccc                                                            <  2:654253/88‑1 (MQ=39)
                                                                                                                                          gctGCGCCTGCCGCGGCAGCCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGtcgccctcgccc                                                            >  1:654253/3‑90 (MQ=39)
                                                                                                                                                |                                                                                                                                             
GTGACGCGACCGCCTTTACCGGTACCTGCAACGGTAGCCAGGTCGATGCCGTTTTCTTCAGCCAGCTTGCGCGCGGCCGGGGCTGCGATCGGGTCATCTTCGCCGGCGTCAGCGGCAGCAGCGGCCGGAGCAGCTGCAGGTGCGGCTGCCGCGGCAGGCGCAGCGGCAGCACCGCCCTCAACGATCGAACCCAGCAATTCGTCGGACAGGACGGTGTCGCCCTCGCCCTTGACGATGTCGCCCAGCACGCCATCGGCAGTAGCCAGTACTTCCAGGACAACCTTGT  >  NC_002947/4732397‑4732682

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑