Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,032,390 C→T A65V (GCG→GTG)  hemG → protoporphyrin oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,032,3900CT100.0% 949.7 / NA 253A65V (GCG→GTG) hemGprotoporphyrin oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (148/105);  total (148/105)

Only 100 of 258 total aligned reads displayed.

CCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAACCGG  >  NC_000913/4032300‑4032479
                                                                                          |                                                                                         
cgaaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAgtgt                                                                                         <  1:844957/91‑1 (MQ=255)
  aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                                                                                       <  1:5936686/93‑1 (MQ=255)
  aaTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTc                                                                                       >  2:1294246/1‑93 (MQ=255)
   aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                                                                                      <  1:6315802/93‑1 (MQ=255)
   aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                                                                                      <  1:705319/93‑1 (MQ=255)
   aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                                                                                      >  2:5400726/1‑93 (MQ=255)
   aTGTGCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTcc                                                                                      <  2:5534529/93‑1 (MQ=255)
     gtgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAg                                                                                    <  2:1383870/93‑1 (MQ=255)
      tgCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAgg                                                                                   <  2:1846294/93‑1 (MQ=255)
       gCACCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGa                                                                                  <  1:1644383/93‑1 (MQ=255)
         aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                                                                                >  2:172159/1‑93 (MQ=255)
         aCCGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAt                                                                                <  2:5019670/93‑1 (MQ=255)
          ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAtt                                                                               >  1:4175942/1‑93 (MQ=255)
          ccGCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAAtt                                                                               >  2:578285/1‑93 (MQ=255)
            gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg                                                                             <  2:2531137/93‑1 (MQ=255)
            gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg                                                                             >  1:3444977/1‑93 (MQ=255)
            gCATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTg                                                                             >  2:4344053/1‑93 (MQ=255)
             cATTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGt                                                                            >  1:656387/1‑93 (MQ=255)
              aTTGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTc                                                                           <  2:105221/93‑1 (MQ=255)
               ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa                                                                          <  1:4376579/93‑1 (MQ=255)
               ttGAAGAACCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCa                                                                          <  1:2148453/93‑1 (MQ=255)
                 gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa                                                                        >  2:412432/1‑93 (MQ=255)
                 gaagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaa                                                                        <  2:1157606/93‑1 (MQ=255)
                  aagaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaa                                                                       >  1:5690364/1‑93 (MQ=255)
                   agaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaa                                                                      >  2:2358501/1‑93 (MQ=255)
                   agaaCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCaaaaa                                                                      >  2:5303163/1‑93 (MQ=255)
                      aacacgGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa                                                                   >  1:4205463/3‑93 (MQ=255)
                      aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa                                                                   <  1:4650334/93‑1 (MQ=255)
                      aCCACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACa                                                                   >  1:296893/1‑93 (MQ=255)
                       ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt                                                                  <  1:2230802/93‑1 (MQ=255)
                       ccACGGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACAt                                                                  <  1:4536429/93‑1 (MQ=255)
                         acggTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGc                                                                <  2:5354641/90‑1 (MQ=255)
                         acggTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGc                                                                >  2:2698811/4‑93 (MQ=255)
                         acggTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGc                                                                <  2:4023278/90‑1 (MQ=255)
                         acggTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGc                                                                >  1:1381782/4‑93 (MQ=255)
                         acggTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGc                                                                >  2:3777684/4‑93 (MQ=255)
                           ggTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGa                                                              >  1:5975318/2‑93 (MQ=255)
                            gTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGAc                                                             <  1:6453734/93‑1 (MQ=255)
                                gAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCgg                                                         >  2:5160058/1‑93 (MQ=255)
                                 aaaaCTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGc                                                        <  1:3539382/93‑1 (MQ=255)
                                 aaaaCTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGc                                                        >  2:4562501/1‑93 (MQ=255)
                                   aaCTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTg                                                      >  1:4422833/1‑93 (MQ=255)
                                    aCTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGa                                                     >  2:263503/1‑93 (MQ=255)
                                    aCTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGa                                                     >  1:1567199/1‑93 (MQ=255)
                                    aCTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGa                                                     >  1:5089940/1‑93 (MQ=255)
                                      tATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAAt                                                   >  1:5386311/1‑93 (MQ=255)
                                      tATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAAt                                                   >  2:5572438/1‑93 (MQ=255)
                                      caTGACCGTGTGGTCATTGCTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCAAGGAATTTGTCAAAAAACATGCGACGCGGCTGAAt                                                   <  2:3377179/92‑1 (MQ=255)
                                         gACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCg                                                <  1:6141359/93‑1 (MQ=255)
                                          aCCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGa                                               >  2:4797436/1‑93 (MQ=255)
                                          aCCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGa                                               <  2:226145/93‑1 (MQ=255)
                                          aCCGAGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGa                                               <  2:3662424/93‑1 (MQ=255)
                                           ccGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGAt                                              >  2:9134/1‑93 (MQ=255)
                                           ccGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGAt                                              <  2:682556/93‑1 (MQ=255)
                                              tgtgGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGcc                                           >  2:6265109/1‑93 (MQ=255)
                                              tgtgGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGcc                                           >  1:3813817/1‑93 (MQ=255)
                                                tgGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGa                                         >  1:2608274/1‑93 (MQ=255)
                                                  gTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAgc                                       <  2:5620482/93‑1 (MQ=255)
                                                  gTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAgc                                       <  2:5200203/93‑1 (MQ=255)
                                                     aTTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGcc                                    >  2:6590462/1‑93 (MQ=255)
                                                     aTTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGcc                                    <  1:3179193/93‑1 (MQ=255)
                                                      ttGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCt                                   >  2:1251504/1‑93 (MQ=255)
                                                      ttGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCt                                   <  2:6059590/93‑1 (MQ=255)
                                                      ttGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCt                                   >  2:169373/1‑93 (MQ=255)
                                                      ttGGTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCt                                   >  1:4993738/1‑93 (MQ=255)
                                                        ggTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCttt                                 >  1:102343/1‑93 (MQ=255)
                                                        ggTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCttt                                 <  1:1115080/93‑1 (MQ=255)
                                                         gTGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCtttt                                >  2:4527568/1‑93 (MQ=255)
                                                          tGCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTa                               <  2:3457007/93‑1 (MQ=255)
                                                           gCTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTAc                              <  1:1822548/93‑1 (MQ=255)
                                                            cTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACt                             <  2:3071964/93‑1 (MQ=255)
                                                            cTTCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACt                             >  2:6029825/1‑93 (MQ=255)
                                                             ttCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTc                            <  2:3525117/93‑1 (MQ=255)
                                                             ttCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTc                            >  1:6419032/1‑93 (MQ=255)
                                                              tCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTcc                           <  2:848172/93‑1 (MQ=255)
                                                              tCTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTcc                           >  1:3923273/1‑93 (MQ=255)
                                                               cTATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCg                          >  1:6077462/1‑93 (MQ=255)
                                                                tATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGt                         <  2:5828372/93‑1 (MQ=255)
                                                                tATTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGt                         >  1:2753843/1‑93 (MQ=255)
                                                                 aTTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTg                        >  1:2638209/1‑93 (MQ=255)
                                                                 aTTCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTg                        >  1:2002016/1‑93 (MQ=255)
                                                                  ttCGCTATGGTCACTACCATTCATTGCTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCCTTTACTCCGTGa                       >  1:2739542/1‑93 (MQ=255)
                                                                  ttCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGa                       >  2:4299724/1‑93 (MQ=255)
                                                                  ttCGCTATGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGa                       <  2:3905188/93‑1 (MQ=255)
                                                                        aTGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCtgg                 >  1:6688000/1‑93 (MQ=255)
                                                                        aTGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCtgg                 >  2:2208483/1‑93 (MQ=255)
                                                                         tGGTCACTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCtggt                >  1:1236389/1‑93 (MQ=255)
                                                                               cTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGgcgcg          >  1:2500174/1‑93 (MQ=255)
                                                                               cTACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGgcgcg          >  2:2948339/1‑93 (MQ=255)
                                                                                tACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGgcgcgc         >  2:3275616/1‑93 (MQ=255)
                                                                                tACCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGgcgcgc         >  1:2314050/1‑93 (MQ=255)
                                                                                 aCCATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGCCGCGCa        >  2:3137245/1‑93 (MQ=255)
                                                                                   cATTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCaaa      >  1:5722101/1‑93 (MQ=255)
                                                                                    aTTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAAc     >  2:6514574/1‑93 (MQ=255)
                                                                                    aTTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAAc     >  1:3250329/1‑93 (MQ=255)
                                                                                    aTTCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGACGAGAGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAAc     >  1:4024470/1‑93 (MQ=255)
                                                                                     atCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGGAAAcc    >  1:6436351/2‑93 (MQ=255)
                                                                                      tCAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAACCg   <  2:3970303/93‑1 (MQ=255)
                                                                                       cAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAACCgg  <  1:2833985/93‑1 (MQ=255)
                                                                                       cAGTGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAACCgg  <  1:3350828/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCAATGTGCACCGCATTGAAGAACCACAGTGGGAAAACTATGACCGTGTGGTCATTGGTGCTTCTATTCGCTATGGTCACTACCATTCAGCGTTCCAGGAATTTGTCAAAAAACATGCGACGCGGCTGAATTCGATGCCGAGCGCCTTTTACTCCGTGAATCTGGTGGCGCGCAAACCGG  >  NC_000913/4032300‑4032479

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