Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,958 G→T 15.4% intergenic (‑45/+28) ridA ← / ← pyrI enamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,470,9580GT15.4% 473.8 / 63.7 183intergenic (‑45/+28)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (80/74);  new base T (10/19);  total (90/93)
Fisher's exact test for biased strand distribution p-value = 1.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 196 total aligned reads displayed.

TACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGC  >  NC_000913/4470860‑4471057
                                                                                                  |                                                                                                   
tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggg                                                                                                    >  1:2633077/1‑100 (MQ=255)
 aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                                   >  2:2670240/1‑100 (MQ=255)
  cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                                  >  1:4039800/1‑100 (MQ=255)
   gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTACTTTATTACCGCGGTAACAAGTCTGGCTCATTATAGGGAGc                                                                                                 >  2:1136154/1‑100 (MQ=255)
      aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAg                                                                                              <  2:240326/100‑1 (MQ=38)
       ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                             <  1:1442708/100‑1 (MQ=38)
       ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                             <  2:3510728/100‑1 (MQ=255)
        gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                            <  2:2151547/100‑1 (MQ=38)
        gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                            <  2:1025628/100‑1 (MQ=255)
        gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                            >  1:2175389/1‑100 (MQ=255)
        gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                            >  1:195259/1‑100 (MQ=255)
         aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                           <  2:1257649/100‑1 (MQ=38)
         aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                           <  1:1045233/100‑1 (MQ=255)
         aCCGATAGCTGCCGGTACATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                           >  1:3682908/1‑100 (MQ=255)
         aCCGATAGCCGCCGGTGCATTTTCCGTCGCTATAGTTTAGCTCATGATTCCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                            <  2:2907445/99‑1 (MQ=255)
          ccGATAGCTGCCGGTGCTTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                          <  2:3811799/100‑1 (MQ=38)
          ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                          <  1:3316202/100‑1 (MQ=38)
           cGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACttt                                                                                         <  2:3216330/100‑1 (MQ=255)
            gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGTAGCCAGACtttt                                                                                        >  2:3328665/1‑100 (MQ=255)
            gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                        <  1:2747102/100‑1 (MQ=255)
              tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                      <  1:2356947/100‑1 (MQ=255)
              tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                      <  1:1052818/100‑1 (MQ=255)
               aGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATGGGGAGCCAGACTTttatt                                                                                     <  2:3961003/100‑1 (MQ=255)
               aGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatt                                                                                     <  1:1015889/100‑1 (MQ=255)
                gCTGCCTGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTTGCTCATTATAGGGAGCCAGACTTttatta                                                                                    <  2:1738377/100‑1 (MQ=255)
                gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtttttacc                                                                                    <  1:1699475/100‑5 (MQ=25)
                gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                                    >  2:2731298/1‑100 (MQ=255)
                   gCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                                 <  1:3212075/100‑1 (MQ=255)
                      ggTGCATTTTCCGTCGCGATAGTTTTGCTAATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                              <  2:1427823/100‑1 (MQ=255)
                       gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                             <  2:695417/100‑1 (MQ=255)
                         gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                           <  1:3913398/100‑1 (MQ=255)
                         gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCATTTAAAACGGCGGTAAAAAGTCTCGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                                 <  2:3384591/94‑1 (MQ=255)
                          cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaa                                                                          >  2:809519/1‑100 (MQ=255)
                           aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGTGATCCAGACTTTTATTACCAACCGCaat                                                                         >  2:178207/1‑100 (MQ=255)
                           aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                         <  2:2011529/100‑1 (MQ=255)
                                 cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCTAGACTTTTATTACCAACCGCAATTAATTg                                                                   >  2:2481980/1‑100 (MQ=255)
                                 cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                    >  2:2364646/1‑99 (MQ=255)
                                 cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTg                                                                   >  2:2256074/1‑100 (MQ=255)
                                 cGTCGCGATAGTTTTGCTCATGATTTCTCCTTGATTACATCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                                    <  2:974682/99‑1 (MQ=255)
                                  gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                  <  2:1594207/100‑1 (MQ=255)
                                  gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                                  <  2:783201/100‑1 (MQ=255)
                                   tagcgATAGTTTTGCTCATGTTTTATCCTTTAGTACAGCGGTAATAAGTCTGGCTCATTATAGGGAGACAGACTTTTATTACCAACCGCAATTAATTgg                                                                  <  2:2698442/97‑1 (MQ=255)
                                    cgcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                                <  1:2743975/100‑1 (MQ=255)
                                    cgcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTGTTATTACCAACCGCAATTAATTGGcc                                                                <  2:1718015/100‑1 (MQ=255)
                                         tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                           <  1:2275676/100‑1 (MQ=255)
                                            ttttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCACCAGCAATTAATTGGCCAGcaac                                                          >  1:4120280/1‑96 (MQ=255)
                                              ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                      <  1:467762/100‑1 (MQ=255)
                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                 <  1:766305/100‑1 (MQ=255)
                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                 <  2:4008684/100‑1 (MQ=255)
                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                 >  2:3119915/1‑100 (MQ=255)
                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                 >  1:1657747/1‑100 (MQ=255)
                                                     tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                               >  2:2560326/1‑100 (MQ=255)
                                                     tGATTTCTCCTTTATTACAACGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                               >  1:1427843/1‑100 (MQ=255)
                                                     tGATTTCTCCTTTATTAAAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTATATTACCAACCGCAATTAATTGGCAAGCAGCACATTATGGGa                                               >  2:371378/1‑100 (MQ=255)
                                                        tttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                            >  1:1321590/1‑100 (MQ=255)
                                                           ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCACCCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                         >  1:3632441/1‑100 (MQ=255)
                                                             ccTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                       >  1:647303/1‑100 (MQ=255)
                                                              cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCttt                                      >  1:3889710/1‑100 (MQ=255)
                                                              cTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATGAATTGGCCAGCACCACATTATGGGAAAACTCttt                                      >  2:2168419/1‑100 (MQ=255)
                                                                 tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                   <  2:3490998/100‑1 (MQ=255)
                                                                 tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                   <  1:1728976/100‑1 (MQ=255)
                                                                 tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                   <  1:1592313/100‑1 (MQ=255)
                                                                  cgtaCAGCGGTAAATAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACGAACCGCAATTAATTGGCCAGCACCACATTATGGGAGAACTCTTTTTc                                   <  2:3583792/97‑1 (MQ=255)
                                                                  attaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                                  >  2:1588130/1‑100 (MQ=255)
                                                                   ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                 <  2:378493/100‑1 (MQ=255)
                                                                   ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                                 <  1:904008/100‑1 (MQ=255)
                                                                    ggttggtaataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                                >  1:198679/11‑100 (MQ=21)
                                                                     aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                               <  1:3436503/100‑1 (MQ=255)
                                                                      cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcacatt                              >  2:553452/1‑98 (MQ=255)
                                                                      cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                               >  2:211998/1‑99 (MQ=255)
                                                                      cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGt                              >  2:3005786/1‑100 (MQ=255)
                                                                        gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                            >  1:563130/1‑100 (MQ=255)
                                                                        gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                            >  2:149311/1‑100 (MQ=255)
                                                                         cGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                           <  2:1164089/100‑1 (MQ=255)
                                                                         cGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                           <  1:1573805/100‑1 (MQ=255)
                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                         <  1:602726/100‑1 (MQ=255)
                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAATa                             >  1:1295118/1‑96 (MQ=255)
                                                                           gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGCAAAACACCTTTCCACAGTAt                            >  1:1787819/1‑97 (MQ=255)
                                                                            tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                        <  1:1892821/100‑1 (MQ=255)
                                                                            tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                        >  1:3565710/1‑100 (MQ=255)
                                                                               aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt                     >  1:2143174/1‑100 (MQ=37)
                                                                                aaGTCTGGCTCAGTATAGGGAGCCAGACTTTTCTTACCAACCGCAATTAATTGGCCAGCACCGCCTTATGGGAAAACTATTTTTCACAGCATCTGCAtt                     >  2:2810701/1‑99 (MQ=255)
                                                                                  gTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTgagatcgg                  >  1:1879426/1‑93 (MQ=37)
                                                                                  gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                  <  1:495569/100‑1 (MQ=255)
                                                                                  gTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGa                  >  1:1551336/1‑100 (MQ=255)
                                                                                   tCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                 >  1:2461960/1‑100 (MQ=38)
                                                                                   tCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg                 >  1:3989003/1‑100 (MQ=38)
                                                                                       gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg             <  2:2987237/100‑1 (MQ=255)
                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa            >  1:3837472/1‑100 (MQ=255)
                                                                                        cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCACCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa            <  2:361342/100‑1 (MQ=255)
                                                                                            ttATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATc        >  1:2869684/1‑100 (MQ=255)
                                                                                              ataGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCAt      >  2:3563942/1‑100 (MQ=255)
                                                                                                aGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTg    <  1:967224/100‑1 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg   <  1:316306/100‑1 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg   >  2:1332538/1‑100 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg   >  1:3489838/1‑100 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg   >  1:723298/1‑100 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg   <  1:1383054/100‑1 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTAAGCGCGATCTCATTg    >  1:3127334/1‑99 (MQ=255)
                                                                                                  ggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGgc  <  1:3771119/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
TACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGGC  >  NC_000913/4470860‑4471057

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: