Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,429,581 G→A 7.5% intergenic (+91/‑389) mqo‑II → / → PP_1252 malate:quinone oxidoreductase/group II intron‑encoding maturase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,429,5810GA 7.5% 173.4 / 2.8 62intergenic (+91/‑389)mqo‑II/PP_1252malate:quinone oxidoreductase/group II intron‑encoding maturase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (3/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACCGCGA  >  NC_002947/1429432‑1429684
                                                                                                                                                     |                                                                                                       
tCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCATACGCGCAAcgcg                                                                                                         <  2:719862/150‑1 (MQ=255)
    tggtgCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTg                                                                                                     <  1:468707/150‑1 (MQ=255)
       tgCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGtt                                                                                                  <  2:56251/150‑1 (MQ=255)
          ggCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTtgtg                                                                                               <  1:24575/150‑1 (MQ=255)
           gcaagcaaGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtgg                                                                                              >  2:51754/1‑150 (MQ=255)
            caagcaagCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                                             >  1:1088225/1‑150 (MQ=255)
            caagcaagCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                                             <  1:423625/150‑1 (MQ=255)
             aagcaagcCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGGCCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggtg                                                                                            >  2:464719/1‑150 (MQ=255)
               gcaagcCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGc                                                                                          <  1:258206/150‑1 (MQ=255)
               gcaagcCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGc                                                                                          <  2:934884/150‑1 (MQ=255)
                   gcCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGATACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAg                                                                                      >  2:492179/1‑150 (MQ=255)
                          gagaGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                                   <  2:610215/146‑1 (MQ=255)
                          gagaGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                                   >  1:610215/1‑146 (MQ=255)
                               cAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTgg                                                                          >  2:847858/1‑150 (MQ=255)
                                aaGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                                         >  1:325524/1‑150 (MQ=255)
                                aaGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGa                                                                         <  2:325524/150‑1 (MQ=255)
                                    cTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                                     <  1:847858/150‑1 (MQ=255)
                                    cTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCACTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                                     >  2:85717/1‑150 (MQ=255)
                                    cTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTGTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTg                                                                     >  1:1390402/1‑150 (MQ=255)
                                     tGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGa                                                                    >  2:444366/1‑150 (MQ=255)
                                            gACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCTACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                                             <  1:812418/150‑1 (MQ=255)
                                            gACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCct                                                             >  1:75587/1‑150 (MQ=255)
                                                tGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTac                                                         <  1:1333123/150‑1 (MQ=255)
                                                  gcgcTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCGAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacac                                                       >  2:957394/1‑150 (MQ=255)
                                                    gcTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                                     >  2:419417/1‑150 (MQ=255)
                                                    gcTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCACTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTacacac                                                     >  2:292583/1‑150 (MQ=255)
                                                     cTGTAATCGGCCATAGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCACTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACAcc                                                    <  1:292583/150‑1 (MQ=255)
                                                      tGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                                        <  1:1286604/129‑1 (MQ=255)
                                                      tGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                                        >  2:1286604/1‑129 (MQ=255)
                                                           tCGGCCATTGTTGTGAAAGAAAGCCACGGGGGGGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                                       <  2:52648/109‑1 (MQ=255)
                                                           tCGGCCATTGTTGTGAAAGAAAGCCACGGGGGGGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCaa                                                                                       >  1:52648/1‑109 (MQ=255)
                                                           tCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                                        >  1:772560/1‑124 (MQ=255)
                                                           tCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAg                                                                        <  2:772560/124‑1 (MQ=255)
                                                              gCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAggg                                           <  1:51754/150‑1 (MQ=255)
                                                                cATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                                            >  2:782688/1‑147 (MQ=255)
                                                                cATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAgg                                            <  1:782688/147‑1 (MQ=255)
                                                                cATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGGTGTGGTGGCCAAGTGGTTGACTTGTAGGTGAAAGGCCTCTACACACCCGGCAAGGGGa                                         >  2:1366151/1‑150 (MQ=255)
                                                                       gtgAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTa                                  <  2:1088225/150‑1 (MQ=255)
                                                                        tgAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg                                 >  2:1282004/1‑150 (MQ=255)
                                                                           aagaaagCCACGGGGGATACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg                              <  1:492179/150‑1 (MQ=255)
                                                                           aagaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCg                              >  2:1355997/1‑150 (MQ=255)
                                                                            agaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                                   <  2:1414163/96‑1 (MQ=255)
                                                                            agaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTgg                                                                                   >  1:1414163/1‑96 (MQ=255)
                                                                             gaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg                                 >  2:38777/1‑145 (MQ=255)
                                                                             gaaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAg                                 <  1:38777/145‑1 (MQ=255)
                                                                              aaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                           >  2:795426/1‑150 (MQ=255)
                                                                              aaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                           <  2:75587/150‑1 (MQ=255)
                                                                              aaagCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAg                           <  1:1239578/150‑1 (MQ=255)
                                                                                  ccACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                                     >  2:637428/1‑88 (MQ=255)
                                                                                  ccACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGt                                                                                     <  1:637428/88‑1 (MQ=255)
                                                                                  ccACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCaa                       <  1:946288/150‑1 (MQ=255)
                                                                                      gggggAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTCGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGgt                   <  2:1023349/150‑1 (MQ=255)
                                                                                           aggccTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTcgc              <  1:464719/147‑1 (MQ=37)
                                                                                             aCCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGTTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgg            <  1:554790/150‑1 (MQ=38)
                                                                                               cTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                                             <  2:1424281/67‑1 (MQ=255)
                                                                                               cTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTgtggt                                                                                             >  1:1424281/1‑67 (MQ=255)
                                                                                               cTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCgggg          >  1:596819/1‑149 (MQ=21)
                                                                                                tCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTg         >  2:907009/1‑150 (MQ=35)
                                                                                                  gTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGAc       >  1:994006/1‑150 (MQ=34)
                                                                                                  gTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCACTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGAc       <  1:85717/150‑1 (MQ=18)
                                                                                                       ttttttGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACCGCGa  <  1:2722/150‑1 (MQ=11)
                                                                                                       ttttttGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCAATTGAGGGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACCGCGa  >  1:260354/1‑150 (MQ=11)
                                                                                                       ttttttGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACACGCTAGTTGGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTAAACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACCGCGa  >  1:260491/1‑150 (MQ=11)
                                                                                                                                                     |                                                                                                       
TCGGTGGTGCGGCAAGCAAGCCGGTCGAGAGCAAGCCTGAGAACGACATGGCGCTGTAATCGGCCATTGTTGTGAAAGAAAGCCACGGGGGAGACCTCGTGGCTTTTTTGTGTGCGCCAGGCATGGCGCGTTGCGCGTATGCGCAACCCGCTTGGTTGTGGTGGCCAAGTGGGTGACTTGGAGGTGAAAGTCCTCTACACACCCGGCAAGGGGAAGTGTTAGCCGGAGGCAAGGGTGTCGCGGGTGACCGCGA  >  NC_002947/1429432‑1429684

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