Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,988,945 C→A 34.6% intergenic (+183/‑363) kefB‑III → / → PP_mr64 glutathione‑regulated potassium/H+ antiporter/C4 ncRNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,988,9450CA34.6% 106.6 / 47.0 61intergenic (+183/‑363)kefB‑III/PP_mr64glutathione‑regulated potassium/H+ antiporter/C4 ncRNA
Reads supporting (aligned to +/- strand):  ref base C (21/18);  new base A (9/12);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCA  >  NC_002947/5988800‑5989086
                                                                                                                                                    |                                                                                                                                             
tgAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCt                                                                                                                                             >  1:411328/1‑148 (MQ=255)
tgAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCt                                                                                                                                             <  2:411328/148‑1 (MQ=255)
       gagcgcCGCCCGCGCGGCGCTTCGCGGGTAACCCCGCTCCTACGTTTGTTTCGGCCCCATGTTTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCagga                                                                                                                                    <  1:529771/148‑1 (MQ=2)
        ggCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggag                                                                                                                                   <  2:363920/150‑1 (MQ=255)
        ggCGCCGCCCGCGCGGCGCTTCGCGCGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGGCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCaggag                                                                                                                                   <  1:57143/150‑1 (MQ=14)
         gcgcCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGTCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggagg                                                                                                                                  <  1:504982/150‑1 (MQ=255)
         gcgcCGCCCCCGCGGCACTTCGCGGGAAAACCCGCCCCTACGTTTGTTTCGTGCCAATTATCCCTGCGGCATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTAGATATCGCGGCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCaggagg                                                                                                                                  <  2:570001/150‑1 (MQ=2)
               cgccggcggcgCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTgg                                                                                                                            <  1:93247/145‑1 (MQ=255)
               cccgcGCGGCGCTTCGCGGGTAAACCCGCTCCTACGGTTGTGTGAGGCCAGTTATTCCTGTGGGACGCGCAGGCGAAATTTACCCCAGCGAGCG‑‑‑CGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACaggagatcggaa                                                                                                                            >  1:779328/1‑143 (MQ=16)
                     cggcgCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGccgaaa                                                                                                                      <  2:93414/150‑5 (MQ=9)
                           ttCGCGGGTAACCCCGCTCCTACGTTCGTTACCGGCAACTTAGGGCTGGGGGACGCGCACGCGA‑‑‑GCCCCATGTCCCGCGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                                                                                <  1:775588/150‑1 (MQ=21)
                              gcgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGAGTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTAGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAg                                                                                                            >  2:47604/1‑150 (MQ=9)
                               cgGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTACTGTGGGACCCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGAGGCTCGCGGCTGGCACGGGAGGAATTGACCAGaaaaaaaccgag                                                                                                            >  2:324207/1‑146 (MQ=25)
                                gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGAtt                                                                                                                              >  1:415183/1‑131 (MQ=255)
                                gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGAtt                                                                                                                              <  2:415183/131‑1 (MQ=255)
                                gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaa                                                                                                                  >  2:654768/1‑143 (MQ=255)
                                gggTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaa                                                                                                                  <  1:654768/143‑1 (MQ=255)
                                 gttAAACCCGCTCCTACGTTTGTTTCGGGCCACATGATTAAGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAg                                                                                                         <  1:439329/148‑1 (MQ=2)
                                  gTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAg                                                                                                         >  2:424708/1‑150 (MQ=255)
                                   tAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                                                                                <  1:522436/142‑1 (MQ=255)
                                   tAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                                                                                >  2:522436/1‑142 (MQ=255)
                                   tAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGc                                                                                                        >  2:455437/1‑150 (MQ=255)
                                    aaaCCCGCTCCTACGTTTGTTTCGGGCCAATTATTCCTCCGGGCTTTGCGCGCCC‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCgg                                                                                                      <  1:398404/150‑1 (MQ=2)
                                       cccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGt                                                                                                    <  1:687689/150‑1 (MQ=255)
                                        ccGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGtt                                                                                                   >  1:118842/1‑150 (MQ=255)
                                          gCTCCTACGTTTTTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAc                                                                                                <  2:368434/150‑1 (MQ=2)
                                          gCTCCTACGTTTGTTTCGGGCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAc                                                                                                <  2:77554/150‑1 (MQ=2)
                                            tCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTgg                                                                                                                                           <  1:67752/106‑1 (MQ=255)
                                            tCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTgg                                                                                                                                           >  2:67752/1‑106 (MQ=255)
                                              cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                                                                                <  2:440426/131‑1 (MQ=255)
                                              cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGaaacaaac                                                                                                                >  1:440426/1‑131 (MQ=255)
                                              cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg                                                                                                            <  1:400833/135‑1 (MQ=255)
                                              cTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAg                                                                                                            >  2:400833/1‑135 (MQ=255)
                                               tACGTTTGTTTCGGGCAAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCcgcg                                                                                           >  2:469513/1‑150 (MQ=2)
                                                      gtttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGcgcc                                                                                    >  1:609278/1‑150 (MQ=9)
                                                       tttCGGGCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGcgccg                                                                                   >  2:436596/1‑150 (MQ=9)
                                                           gggCCAATTATTCCTGCGGGATTTGCGCGCGA‑‑‑CCGCCTTGGGGCATGTCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGGCTGGCCGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCcgcgcg                                                                               >  1:155630/1‑150 (MQ=2)
                                                                               aCGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATg                                                            >  1:674720/1‑150 (MQ=255)
                                                                                cgcgCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATgg                                                           <  1:424708/150‑1 (MQ=255)
                                                                                              gCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCaaaaa                                                >  2:784424/1‑150 (MQ=255)
                                                                                                   cTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATc                                           <  2:118842/150‑1 (MQ=255)
                                                                                                        gcatgtctcgataTCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgat                                                                  >  1:700693‑M2/7‑90 (MQ=255)
                                                                                                        gcatgtctcgataTCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgat                                                                  <  2:700693‑M2/115‑32 (MQ=255)
                                                                                                        gcACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCAGCGCGGGCGGCGGTCGATCTCATTGGCGCCCAAAAACTATCGCCAg                                      >  2:728886/1‑150 (MQ=255)
                                                                                                               tagatATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagc                                                                                       >  2:290721‑M2/3‑83 (MQ=255)
                                                                                                               tCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagc                                                                                       <  1:290721‑M2/93‑11 (MQ=255)
                                                                                                                 gATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGa                             >  2:345793/1‑150 (MQ=255)
                                                                                                                 gATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGACCAAAAACTTTTGGCAGCCACCTGGa                             >  2:382074/1‑150 (MQ=21)
                                                                                                                       gcgTCGTACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattcca                      >  1:47604‑M2/1‑75 (MQ=255)
                                                                                                                         gtcgtACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATcc                     >  2:559292/1‑150 (MQ=255)
                                                                                                                             tACAAACAAGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggagctagatctcacattcgctgaaaatcctccggcaaacacctcaaaattccaacgccg                <  1:651744‑M2/150‑82 (MQ=255)
                                                                                                                                aaacaaGGCGGTCGCGCGCGCCTGACATGAGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAAtt              >  1:543262/1‑150 (MQ=255)
                                                                                                                                aaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaa                           <  2:640484‑M2/136‑71 (MQ=255)
                                                                                                                                aaacaaGGCGGTCGCGCGCGACTGGTGCAGGAGGGACTGGCCGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaa                           >  1:640484‑M2/1‑66 (MQ=255)
                                                                                                                                 aacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcccgcgcggcgcttcgc                                                                  >  2:564724‑M2/1‑66 (MQ=255)
                                                                                                                                 aacaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTAcccgcgaagcgccgcccgcgcggcgcttcgc                                                                  <  1:564724‑M2/97‑32 (MQ=255)
                                                                                                                                  acaaGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTc            >  1:390109/1‑150 (MQ=255)
                                                                                                                                        cGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACg      <  2:390109/150‑1 (MQ=255)
                                                                                                                                          gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGc                                         >  2:94895/1‑113 (MQ=255)
                                                                                                                                          gTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGc                                         <  1:94895/113‑1 (MQ=255)
                                                                                                                                            cgcgcgcgCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCa  <  1:345793/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
TGAGATTGGGCGCCGCCCGCGCGGCGCTTCGCGGGTAAACCCGCTCCTACGTTTGTTTCGGGCCAGTTATTCCTGTGGGACGCGCACGCGA‑‑‑GCCCCCTGGCGCACGCCTCGATATCGCGTCGTACAAACAAGGCGGTCGCGCGCGCCTGGCACAGGAGGGATTGGCCAGAAACAAACGTAGGAGCGGGTTTACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCATGGGCGCCCAAAAACTATCGCCAGGCACCTGGAAGCCATCCTGCAATTTCCTGACGATCA  >  NC_002947/5988800‑5989086

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