Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,268,999 T→A 15.0% intergenic (+19/+17) PP_4643 → / ← radA xanthine/uracil permease family protein/DNA repair enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,268,9990TA15.0% 155.5 / 18.8 60intergenic (+19/+17)PP_4643/radAxanthine/uracil permease family protein/DNA repair enzyme
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base A (5/4);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

CGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCTTC  >  NC_002947/5268851‑5269148
                                                                                                                                                     |                                                                                                                                                     
cgcgCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCtt                                                                                                                                                      >  2:1039414/1‑150 (MQ=255)
   gcTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgct                                                                                                                                                   <  1:312591/150‑1 (MQ=255)
   gcTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgct                                                                                                                                                   <  1:1181114/150‑1 (MQ=255)
    cTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctc                                                                                                                                                  <  2:999522/150‑1 (MQ=255)
      gCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTgctcgc                                                                                                                                                <  1:645405/150‑1 (MQ=255)
        gTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAgctcgctc                                                                                                                                              <  1:942634/150‑1 (MQ=255)
          gATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCt                                                                                                                                            <  2:359450/150‑1 (MQ=255)
          gATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCAAGCTCGCTCCt                                                                                                                                            <  2:1071087/150‑1 (MQ=255)
                 actactTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCt                                                                                                                                     >  2:612101/1‑150 (MQ=255)
                     ctTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTc                                                                                                                                 >  2:541365/1‑150 (MQ=255)
                      ttCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCg                                                                                                                                >  2:104016/1‑150 (MQ=255)
                         aaGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAAc                                                                                                                             <  1:612101/150‑1 (MQ=255)
                          aGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACa                                                                                                                            <  1:1009253/150‑1 (MQ=255)
                            gATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAAt                                                                                                                          <  1:545652/150‑1 (MQ=255)
                                       gTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAgcgc                                                                                                               <  1:697171/150‑1 (MQ=255)
                                       gTTCGATGCCGAGGGTAACGAAGTGACGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGCCGGGTTGCAGGAGCGAGGTTGCGTGCGCCTTCACGCTATTCGAACAATGCATCCAgcgc                                                                                                               >  1:499581/1‑150 (MQ=255)
                                        ttCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGcc                                                                                                              >  2:908396/1‑150 (MQ=255)
                                                   gggTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCg                                                                                                   >  2:227646/1‑150 (MQ=255)
                                                    ggTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACTGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCgg                                                                                                  <  1:559742/150‑1 (MQ=255)
                                                       aaCGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGAACGGCCGTTGCCAAGCACTGTCCGATGCGGTGGGCTAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTa                                                                                               >  1:1051358/1‑150 (MQ=255)
                                                        aCGAAGCGCCGATGATGGAACTGGACGAGCACGGCAATGAACCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGCGTTGCAGGCGCGCGCTTGTACGATCGTGCCCCGTATTCGAATAATGCGTCCATCCGTTGTGCCAGCCGGGTaa                                                                                              <  2:499581/150‑1 (MQ=255)
                                                                    tGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCt                                                                                  <  1:908396/150‑1 (MQ=255)
                                                                           aCTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAacc                                                                           >  1:944924/1‑150 (MQ=255)
                                                                            cTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCggg                                                                                                 >  1:156411/1‑127 (MQ=255)
                                                                            cTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCggg                                                                                                 <  2:156411/127‑1 (MQ=255)
                                                                                      aCGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCgg                                                                                                  >  2:865070/1‑116 (MQ=255)
                                                                                      aCGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCgg                                                                                                  <  1:865070/116‑1 (MQ=255)
                                                                                        ggCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGcttct                                                              >  1:1184631/1‑150 (MQ=255)
                                                                                          cAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCttt                                                            >  1:383565/1‑150 (MQ=255)
                                                                                          cAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCttt                                                            <  2:1051358/150‑1 (MQ=255)
                                                                                           aaTGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAACCCGCAGGCGCTTCTTTc                                                           <  2:944924/150‑1 (MQ=255)
                                                                                              gAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTc                                                           >  1:682722/1‑147 (MQ=255)
                                                                                              gAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTc                                                           <  2:682722/147‑1 (MQ=255)
                                                                                              gAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTCCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGgc                                                        >  2:1171482/1‑150 (MQ=255)
                                                                                                 aCGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCa                                                     >  1:104056/1‑150 (MQ=255)
                                                                                                         gcgcACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCtt                                             >  1:639957/1‑150 (MQ=255)
                                                                                                                   gTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCgg                                                                                                  >  1:995332/1‑87 (MQ=255)
                                                                                                                   gTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCgg                                                                                                  <  2:995332/87‑1 (MQ=255)
                                                                                                                   gTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATg                                   <  2:104056/150‑1 (MQ=255)
                                                                                                                     tGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGc                                 <  1:104016/150‑1 (MQ=255)
                                                                                                                           gCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCtt                           <  2:639957/150‑1 (MQ=255)
                                                                                                                               tGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCCGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGAt                                    >  2:644571/1‑137 (MQ=255)
                                                                                                                               tGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCt                                                                                  >  1:362548/1‑91 (MQ=255)
                                                                                                                               tGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCt                                                                                  <  2:362548/91‑1 (MQ=255)
                                                                                                                               tGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGAt                                    <  1:644571/137‑1 (MQ=255)
                                                                                                                               tGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAg                       >  2:290137/1‑150 (MQ=255)
                                                                                                                                   gggTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGt                   <  1:1039414/150‑1 (MQ=255)
                                                                                                                                     gTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAc                                                                             <  1:1100298/90‑1 (MQ=255)
                                                                                                                                     gTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAAc                                                                             >  2:1100298/1‑90 (MQ=255)
                                                                                                                                         cAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTgg             >  1:50235/1‑150 (MQ=255)
                                                                                                                                         cAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTgg             <  1:227646/150‑1 (MQ=255)
                                                                                                                                          aGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCa                                        <  2:13664/122‑1 (MQ=255)
                                                                                                                                          aGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCa                                        >  1:13664/1‑122 (MQ=255)
                                                                                                                                          aGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGc            >  2:568481/1‑150 (MQ=255)
                                                                                                                                          aGGAGCGAGCAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTCCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGc            >  1:353297/1‑150 (MQ=255)
                                                                                                                                             gagcgagcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCg          >  1:200660/1‑150 (MQ=255)
                                                                                                                                             gagcgagcAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTg                                            >  2:351929/1‑116 (MQ=255)
                                                                                                                                             gagcgagcAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTg                                            <  1:351929/116‑1 (MQ=255)
                                                                                                                                              agcgagcAAGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGc         >  2:1118540/1‑150 (MQ=255)
                                                                                                                                                   gcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttcct    >  2:410400/1‑150 (MQ=255)
                                                                                                                                                   gcTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttcct    >  1:150920/1‑150 (MQ=255)
                                                                                                                                                     ttGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGcttccttc  >  1:608639/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCTGCGTGATGGGGACTACTTCAAGGATGGCAAGCTGTTCGATGCCGAGGGTAACGAAGTGCCGATGCTGGAACTGGACGAGCACGGCAATGAAACGGTCAGGCGCACGGCCGTTGCCGAGCACTGACGGGTTGCAGGAGCGAGCTTGCTCGCTCCTGCACGCTATTCGAACAATGCATCCAGCGCCTGTTCCAGCCGGGTAACGGCAATCACCTGCAACCCCGCAGGCGCTTCTTTCGGCGCATTGCCCTTGGGCACGATGGCACGCTTGAAGCCGTGCTTGGCCGCTTCCTTC  >  NC_002947/5268851‑5269148

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