Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,923,308 T→C G159G (GGA→GGG)  truC ← tRNA(Ile1,Asp) pseudouridine(65) synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,923,3080TC98.3% 192.5 / ‑3.8 59G159G (GGA→GGG) truCtRNA(Ile1,Asp) pseudouridine(65) synthase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (33/25);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.73e-01

AATGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGC  >  NC_000913/2923162‑2923451
                                                                                                                                                  |                                                                                                                                               
aaTGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCgg                                                                                                                                             <  1:201943/151‑1 (MQ=255)
      aaCTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCg                                                                                                                                       >  2:37344/1‑151 (MQ=255)
          gactgGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCt                                                                                                                                    <  2:709633/150‑1 (MQ=255)
           actgGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTcc                                                                                                                                  >  2:711011/1‑151 (MQ=255)
           actgGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTcc                                                                                                                                  <  2:658527/151‑1 (MQ=255)
                  tGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTcc                                                                                                                            >  1:117698/1‑150 (MQ=255)
                   gTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTcca                                                                                                                           <  2:404875/150‑1 (MQ=255)
                      gCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAg                                                                                                                       <  1:756505/151‑1 (MQ=255)
                            aCCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGt                                                                                                                  >  1:319516/1‑150 (MQ=255)
                               gCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAgcg                                                                                                              >  1:535562/1‑151 (MQ=255)
                                   gAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTg                                                                                                          >  2:276860/1‑151 (MQ=255)
                                      gCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTc                                                                                                       >  2:455684/1‑151 (MQ=255)
                                         aaaaTGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCggg                                                                                                    >  1:810771/1‑151 (MQ=255)
                                            aTGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTaa                                                                                                 <  1:349320/151‑1 (MQ=255)
                                            aTGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTa                                                                                                  >  1:48294/1‑150 (MQ=255)
                                               cTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGt                                                                                              >  2:769117/1‑151 (MQ=255)
                                                 cAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTcc                                                                                            >  1:76942/1‑151 (MQ=255)
                                                  aGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCg                                                                                           >  2:395336/1‑151 (MQ=255)
                                                             gATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCgg                                                                                >  2:462248/1‑151 (MQ=255)
                                                              aTTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGc                                                                               <  2:509946/151‑1 (MQ=255)
                                                                  tGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCAttt                                                                           >  1:467592/1‑151 (MQ=255)
                                                                   ggCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTc                                                                          <  2:41651/151‑1 (MQ=255)
                                                                   ggCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCAttt                                                                           <  1:680064/150‑1 (MQ=255)
                                                                   ggCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCAtt                                                                            <  2:337339/149‑1 (MQ=255)
                                                                    gCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCt                                                                         >  1:421916/1‑151 (MQ=255)
                                                                               cATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGcc                                                                >  2:839847/1‑149 (MQ=255)
                                                                                          tCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATaa                                                    >  2:184842/1‑150 (MQ=255)
                                                                                           cGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATg                                                  <  2:859356/151‑1 (MQ=255)
                                                                                            gCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGc                                                 >  1:300973/1‑151 (MQ=255)
                                                                                                ataatCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTc                                              >  2:598620/1‑150 (MQ=255)
                                                                                                   atCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCActg                                          >  1:413626/1‑151 (MQ=255)
                                                                                                      ggATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCActgctg                                       <  2:252869/151‑1 (MQ=255)
                                                                                                       gATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTgctgg                                      <  1:33033/151‑1 (MQ=255)
                                                                                                       gATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTgctgg                                      >  1:778131/1‑151 (MQ=255)
                                                                                                       gATGACGCAAATGCGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTgctgg                                      <  2:212222/151‑1 (MQ=255)
                                                                                                            cGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCg                                 >  1:526397/1‑151 (MQ=255)
                                                                                                             gCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCgg                                >  2:543740/1‑151 (MQ=255)
                                                                                                                aaTGGGCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGcc                             >  2:192404/1‑151 (MQ=255)
                                                                                                                   gggCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCttt                          <  1:57269/151‑1 (MQ=255)
                                                                                                                     gCAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTAt                        <  1:106600/151‑1 (MQ=255)
                                                                                                                      cAAGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATc                       >  1:64094/1‑151 (MQ=255)
                                                                                                                       aaGATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCt                      >  2:228869/1‑151 (MQ=255)
                                                                                                                         gATGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTc                    <  2:319539/151‑1 (MQ=255)
                                                                                                                          aTGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTcg                   <  2:190281/151‑1 (MQ=255)
                                                                                                                           tGTCGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTcgc                  >  2:752544/1‑151 (MQ=255)
                                                                                                                              cGGCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCggg               >  2:96740/1‑151 (MQ=255)
                                                                                                                               ggCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGc              <  1:701378/151‑1 (MQ=255)
                                                                                                                               ggCGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGACGTAGCGCGTGGTCGGGCAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCAAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTtgc                  >  1:458257/1‑147 (MQ=255)
                                                                                                                                 cGAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCaa            <  2:213077/151‑1 (MQ=255)
                                                                                                                                  gAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCGGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCaaa           <  1:601630/151‑1 (MQ=255)
                                                                                                                                  gAAGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCaaa           <  1:160824/151‑1 (MQ=255)
                                                                                                                                   aaGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAt          >  1:199215/1‑151 (MQ=255)
                                                                                                                                    aGCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAtt         >  2:232367/1‑151 (MQ=255)
                                                                                                                                     gCTGATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAAtt         <  1:652559/150‑1 (MQ=255)
                                                                                                                                        gATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTa       <  1:625517/149‑1 (MQ=255)
                                                                                                                                        gATGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTAt      <  1:196633/150‑1 (MQ=255)
                                                                                                                                         aTGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATc     >  2:677978/1‑150 (MQ=255)
                                                                                                                                         aTGTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCa    <  2:49566/151‑1 (MQ=255)
                                                                                                                                           gTTTGCGCCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGc  >  1:94312/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
AATGACAACTGACTGGCATGTAGCATTAACCGCTGGAGGCCAAAATGCTCAGCACCGCTGCGATTCTGGCGTAAATCGCCATGTTTGCTATCGCCAATAATCGGATGACGCAAATGGGCAAGATGTCGGCGAAGCTGATGTTTGCGTCCGGTTTTCGGCTCCAGTTCCACCAGGCCGTAGCGCGTGGTCGGGTAACGTCCGGTCGCTACCGGCATTTCTACGGTCGCCAGACCGCGATAATGCGTCACTGCTGGCTGCGGGCCTTTATCTTCGCGGGCAAATTTATCAGC  >  NC_000913/2923162‑2923451

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑