Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,568,848 C→T A408T (GCC→ACC)  gadB ← glutamate decarboxylase isozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,568,8480CT99.2% 327.7 / NA 133A408T (GCC→ACC) gadBglutamate decarboxylase isozyme
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (55/77);  total (56/77)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

ATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAA  >  NC_000913/1568799‑1568937
                                                 |                                                                                         
aTAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCa                                                <  1:667082/93‑1 (MQ=14)
 tAATGCGCATCACCTCGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                               >  1:4067467/1‑93 (MQ=12)
 tAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                               <  1:1666630/93‑1 (MQ=17)
 tAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCAt                                               <  2:4614590/93‑1 (MQ=17)
  aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCATCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              <  2:2248769/93‑1 (MQ=255)
  aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              >  1:6641725/1‑93 (MQ=18)
  aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              >  2:1395914/1‑93 (MQ=18)
  aaTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATa                                              <  1:1250864/93‑1 (MQ=18)
   aTGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATAc                                             <  2:3673939/93‑1 (MQ=18)
    tGCTCACCACCTCGATGTCGGTGGCTTCGCCGCCGAGAGTGAAGGGCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACa                                            <  2:3097612/93‑1 (MQ=11)
    agcgcATCACCGCGATGTCGGTGGCTACACCGCCGAAAGTGAAGGTCGGACCCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACa                                            <  2:6623945/92‑1 (MQ=11)
     gcgcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAg                                           <  2:4204289/93‑1 (MQ=21)
       gcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         <  2:4419268/93‑1 (MQ=21)
       gcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         >  1:3327914/1‑93 (MQ=18)
       gcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         <  1:5810977/93‑1 (MQ=21)
       gcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         >  2:4176990/1‑93 (MQ=18)
       gcATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         >  2:115256/1‑93 (MQ=18)
       gcATCACCACGAGGTCGGTGGCTTCCCCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAggg                                         <  1:3877547/93‑1 (MQ=255)
        cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgt                                        <  1:4777536/93‑1 (MQ=21)
        cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgt                                        >  2:3944417/1‑93 (MQ=18)
        cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgt                                        >  2:4573023/1‑93 (MQ=18)
        cATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgt                                        <  1:4128902/93‑1 (MQ=21)
         aTCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtg                                       >  1:4087164/1‑93 (MQ=18)
         aTCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtg                                       >  1:5585300/1‑93 (MQ=18)
         aTCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtg                                       <  2:4888372/93‑1 (MQ=21)
         aTCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtg                                       <  2:395944/93‑1 (MQ=18)
          tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                                      >  1:1909463/1‑93 (MQ=18)
          tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                                      >  2:4653632/1‑93 (MQ=18)
          tCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGgtgt                                      <  1:5347786/93‑1 (MQ=21)
           ccccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTa                                     <  2:1141302/91‑1 (MQ=16)
           caccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTa                                     <  2:6539463/93‑1 (MQ=21)
           caccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTa                                     <  2:4112071/93‑1 (MQ=21)
            accacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTAt                                    >  2:2788776/1‑93 (MQ=18)
            accacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTAt                                    >  1:3681466/1‑93 (MQ=18)
             ccacGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATc                                   <  2:1454741/93‑1 (MQ=18)
               ccGATGTCGGTGTCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                                 <  2:4276754/92‑1 (MQ=255)
               ccGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                                 <  1:4407900/92‑1 (MQ=18)
               acGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                                 <  1:1748246/93‑1 (MQ=18)
               acGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATccc                                 >  1:520727/1‑93 (MQ=18)
                cGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCg                                >  2:4203956/1‑93 (MQ=18)
                 gATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCgg                               >  2:6297913/1‑93 (MQ=18)
                 gATGTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCgg                               >  1:1079041/1‑93 (MQ=18)
                   tgtgGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGAt                             <  1:2467134/93‑1 (MQ=255)
                    gTCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATc                            <  1:4314357/93‑1 (MQ=21)
                     tCGGTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCt                           <  1:318508/93‑1 (MQ=21)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc                         >  1:5748418/1‑93 (MQ=18)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc                         >  2:6308575/1‑93 (MQ=18)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc                         <  2:5687/93‑1 (MQ=21)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc                         <  2:4222411/93‑1 (MQ=21)
                       ggTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTc                         <  1:3680400/93‑1 (MQ=21)
                        gTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa                        <  1:659556/93‑1 (MQ=21)
                        gTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa                        >  2:5172743/1‑93 (MQ=18)
                        gTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa                        >  1:1944606/1‑93 (MQ=18)
                        gTGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCa                        >  2:736769/1‑93 (MQ=18)
                         tGGCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCAc                       <  1:1087794/93‑1 (MQ=21)
                          ggCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCAcc                      <  2:2469979/93‑1 (MQ=18)
                          ggCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCAcc                      >  2:1401848/1‑93 (MQ=18)
                           gCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCa                     >  1:2743719/1‑93 (MQ=18)
                           gCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCa                     <  2:4096031/93‑1 (MQ=21)
                           gCTTCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCa                     >  1:1009902/1‑93 (MQ=18)
                            ctccACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCAt                    <  2:2958143/90‑1 (MQ=16)
                            cTTCCCCGCCGAGAGTGAAAGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCAt                    <  1:3833110/93‑1 (MQ=255)
                             ttCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATc                   >  1:4236145/1‑93 (MQ=18)
                             ttCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATc                   >  1:5563005/1‑93 (MQ=18)
                              tCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCt                  <  1:719939/93‑1 (MQ=21)
                              tCACCGCCGAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCt                  <  1:894312/93‑1 (MQ=21)
                                 ccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTc               <  2:3380140/93‑1 (MQ=18)
                                 ccgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTc               <  2:6093366/93‑1 (MQ=18)
                                  cgctgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCa              <  1:115684/93‑1 (MQ=16)
                                  cgccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCa              <  2:5033946/93‑1 (MQ=21)
                                   gccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAg             >  1:4277308/1‑93 (MQ=18)
                                   gccgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAg             >  2:5997293/1‑93 (MQ=18)
                                     cgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGtt           >  1:6379012/1‑93 (MQ=18)
                                     cgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGtt           <  1:505911/93‑1 (MQ=21)
                                     cgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGtt           <  1:4694672/93‑1 (MQ=21)
                                     cgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGtt           >  1:5880695/1‑93 (MQ=18)
                                     cgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGtt           <  1:4260448/93‑1 (MQ=21)
                                     cgAGAGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGAGTATCCCGGATCTTCACCATCTTTCAGtt           >  2:2111553/1‑93 (MQ=16)
                                      gagaGTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGttt          >  1:2810750/1‑93 (MQ=18)
                                        gagTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGa        <  2:5755923/93‑1 (MQ=21)
                                        gagTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGa        <  2:2297987/93‑1 (MQ=21)
                                        gagTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGa        <  1:1138831/93‑1 (MQ=21)
                                         agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGaa       <  1:6351623/93‑1 (MQ=18)
                                         agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGaa       >  1:3375874/1‑93 (MQ=18)
                                         agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGaa       >  2:162618/1‑93 (MQ=18)
                                         agTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGaa       <  2:3459423/93‑1 (MQ=21)
                                          gTGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAg      >  2:48958/1‑93 (MQ=18)
                                           tGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGc     <  2:1589985/93‑1 (MQ=21)
                                           tGAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGc     >  1:6324409/1‑93 (MQ=18)
                                            taaGGTCGGACCCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCa    <  1:3437978/92‑1 (MQ=255)
                                            gAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCa    >  1:4653853/1‑93 (MQ=18)
                                            gAAGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCa    <  2:5204031/93‑1 (MQ=21)
                                             aaGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaa   >  2:1632605/1‑93 (MQ=18)
                                             aaGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaa   >  2:3430800/1‑93 (MQ=18)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  <  2:131614/93‑1 (MQ=21)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  >  2:6189192/1‑93 (MQ=18)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  >  2:602474/1‑93 (MQ=18)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  >  1:1582544/1‑93 (MQ=18)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  <  1:3812429/93‑1 (MQ=21)
                                              aGGTCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCaaa  <  2:2240436/93‑1 (MQ=18)
                                                 |                                                                                         
ATAATGCGCATCACCACGATGTCGGTGGCTTCACCGCCGAGAGTGAAGGCCGGAACCTGCCAGCCGCGCAGACGCAGACGTTCAGAGAGGTCATACAGGGTGTATCCCGGATCTTCACCATCTTTCAGTTTGAAGCAAA  >  NC_000913/1568799‑1568937

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