Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,806,160 A→T 11.9% S593T (TCG→ACG)  PP_4234 ← bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,806,1600AT11.9% 154.9 / 12.8 60S593T (TCG→ACG) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base T (3/4);  total (29/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 = 4.32e-01

CCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTC  >  NC_002947/4806013‑4806304
                                                                                                                                                   |                                                                                                                                                
ccGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAAc                                                                                                                                                >  2:912755/1‑150 (MQ=255)
ccGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAAc                                                                                                                                                <  2:81578/150‑1 (MQ=255)
      aTGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGa                                                                                                                                          <  2:793355/150‑1 (MQ=255)
          cGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGGGTCCAGGCGTGTCGCGGTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgag                                                                                                                                      <  2:578162/150‑2 (MQ=255)
          cGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgag                                                                                                                                      <  2:300338/150‑2 (MQ=255)
           gACGCCAGCTTGCCGACCCAGAGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGaaa                                                                                                                                     >  2:387456/1‑150 (MQ=255)
                              aGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCt                                                                                                                  >  1:191424/1‑150 (MQ=255)
                                cgcgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                <  1:934921/150‑1 (MQ=255)
                                  cgGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGa                                                                                                              <  2:948761/150‑1 (MQ=255)
                                     gTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGTCTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGt                                                                                                           <  2:512119/150‑1 (MQ=255)
                                     gTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCcgacgag                                                                                                                                      <  1:751057/123‑2 (MQ=255)
                                     gTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGacgag                                                                                                                                      >  2:751057/1‑119 (MQ=255)
                                      ttGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTGTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTg                                                                                                          <  2:247739/150‑1 (MQ=255)
                                             gcgcGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGc                                                                                                   <  1:621284/150‑1 (MQ=255)
                                                 ggcgTCCTGCTCTGGCTAGCTGATCAGGCACGTGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATGCAGGGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGt                                                                                               <  1:347210/150‑1 (MQ=255)
                                                 ggcgTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGt                                                                                               >  1:298646/1‑150 (MQ=255)
                                                 ggcgTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGt                                                                                               <  1:1045477/150‑1 (MQ=255)
                                                      ccTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCa                                                                                          <  1:327178/150‑1 (MQ=255)
                                                       cTGCTCCGGCTGGCTGATCAGGTACGGGTTTGCGCTGGTGCTCAATTGCAGCTCTGGACATTCGGTGACTGCAGTCTTCACATCGTGGGCCGAACGGCCGAGGAAATCGCTGACAACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAg                                                                                         >  1:497064/1‑150 (MQ=255)
                                                        tGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGacgagg                                                                                                                                     >  1:240045/1‑100 (MQ=255)
                                                        tGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGTTCGGCCGacgagg                                                                                                                                     <  2:240045/105‑6 (MQ=255)
                                                             cGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGTCCGGACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGt                                                                                   >  1:619282/1‑150 (MQ=255)
                                                               gctggctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGc                                                                                 <  2:234607/150‑1 (MQ=255)
                                                                   gctgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGat                                                                             >  2:720336/1‑150 (MQ=255)
                                                                     tgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCGGAACGGCCGAGGAAATCGCTGACCACAAAAGAATGGGAGTTGGACATGCCGGTGCCCAGCTCGGTGCCGatat                                                                           >  1:518830/1‑150 (MQ=255)
                                                                     tgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGatat                                                                           >  1:593862/1‑150 (MQ=255)
                                                                     tgATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGatat                                                                           >  1:1001939/1‑150 (MQ=255)
                                                                        tCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATggg                                                                        >  1:836908/1‑149 (MQ=255)
                                                                        tCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATgaa                                                                        <  1:1062523/150‑3 (MQ=255)
                                                                        tCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGgc                                                                        <  1:387456/150‑1 (MQ=255)
                                                                            gTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAg                                                                    >  1:1006103/1‑150 (MQ=255)
                                                                            gTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAg                                                                    >  2:976219/1‑150 (MQ=255)
                                                                            gTACGGGTTTCCGCTGGTGCTCAATTGCAGCCCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAg                                                                    >  2:252253/1‑150 (MQ=255)
                                                                                 ggTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACTAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGa                                                               <  2:214634/150‑1 (MQ=255)
                                                                                    ttCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGc                                                            <  2:623985/150‑1 (MQ=255)
                                                                                       cGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGc                                                         <  1:865537/150‑1 (MQ=255)
                                                                                          tGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAATGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCgg                                                          <  2:617144/146‑1 (MQ=255)
                                                                                          tGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAATGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCgg                                                          >  1:617144/1‑146 (MQ=255)
                                                                                            gTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                <  1:10427/90‑1 (MQ=255)
                                                                                            gTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTgg                                                                                                                >  2:10427/1‑90 (MQ=255)
                                                                                             tGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCgg                                                   >  2:627453/1‑150 (MQ=255)
                                                                                                               gCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTgg                                 >  1:71702/1‑150 (MQ=255)
                                                                                                                       gTGACTGCAGTCTTCCCCTCGGCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGgc                         >  2:831271/1‑150 (MQ=255)
                                                                                                                       gTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGgc                         <  1:912755/150‑1 (MQ=255)
                                                                                                                            tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAg                                                                    >  2:318321/1‑102 (MQ=255)
                                                                                                                            tGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAg                                                                    <  1:318321/102‑1 (MQ=255)
                                                                                                                                 tCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAc               >  1:287224/1‑150 (MQ=255)
                                                                                                                                 tCTTCACCTCGTCGGACGAACGGCCGAGGAAATCGCGGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCTAGCTCGGGGCCGATATGGCGCAGGCTGATGCGGCAATCGGCGCTGAACTGGACACTGGCCATCGGCGCCTCGGAAc               <  1:309655/150‑1 (MQ=255)
                                                                                                                                  cTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAAcc              >  2:782034/1‑150 (MQ=255)
                                                                                                                                   ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt             >  1:305105/1‑150 (MQ=255)
                                                                                                                                   ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt             >  1:655355/1‑150 (MQ=255)
                                                                                                                                   ttCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCATGGGAGGCGGAAATCAGCATGCCCGGCTCGGTGCCGATAGGGCGCAGGCTGATGCGGCAATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCt             <  1:252253/150‑1 (MQ=255)
                                                                                                                                       ccTCGTCGGCCGAACGGCCGAGGAAATCGCTGTCCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCACAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTAc         >  2:626644/1‑150 (MQ=255)
                                                                                                                                       ccTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTAc         <  2:1006103/150‑1 (MQ=255)
                                                                                                                                         tcgtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGgc                         <  1:49702/132‑1 (MQ=255)
                                                                                                                                         tcgtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGgc                         >  2:49702/1‑132 (MQ=255)
                                                                                                                                           gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACCCTGGCCATCGGCGCCTCGGAACCTGTACCGaa     <  2:1031566/150‑1 (MQ=255)
                                                                                                                                           gtcgGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa     <  2:71702/150‑1 (MQ=255)
                                                                                                                                           ctcgGCCGTTCGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGaa     >  1:802562/2‑150 (MQ=255)
                                                                                                                                              ggCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCACAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTc  <  1:626644/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGACGATGGCGACGCCAGCTTGCCGACCCAGCGCGGGTTGCGCAGCGCGGCGTCCTGCTCCGGCTGGCTGATCAGGTACGGGTTTCCGCTGGTGCTCAATTGCAGCTCTGGCCATTCGGTGACTGCAGTCTTCACCTCGTCGGCCGAACGGCCGAGGAAATCGCTGACCACCAAAGCCTGGGAGGTGGACATGCCGGTGCCCAGCTCGGTGCCGATATGGCGCAGGCTGATGCGGCCATCGGCGCTGAACTCGACACTGGCCATCGGCGCCTCGGAACCTGTACCGAAGTC  >  NC_002947/4806013‑4806304

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: ‑