Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,957 G→A 15.2% intergenic (‑44/+29) 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,9570GA15.2% 440.4 / 56.4 183intergenic (‑44/+29)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (81/74);  new base A (9/19);  total (90/93)
Fisher's exact test for biased strand distribution p-value = 6.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 196 total aligned reads displayed.

TACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470860‑4471056
                                                                                                 |                                                                                                   
tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggg                                                                                                   >  1:2633077/1‑100 (MQ=255)
 aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAg                                                                                                     <  1:3717694/97‑1 (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:1052818/100‑1 (MQ=255)
              tAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttat                                                                                     <  1:2356947/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:1657747/1‑100 (MQ=255)
                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                <  2:4008684/100‑1 (MQ=255)
                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                <  1:766305/100‑1 (MQ=255)
                                                   cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                                >  2:3119915/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                                  <  1:1728976/100‑1 (MQ=255)
                                                                 tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                                  <  2:3490998/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                           >  2:149311/1‑100 (MQ=255)
                                                                        gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                           >  1:563130/1‑100 (MQ=255)
                                                                         cGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                          <  1:1573805/100‑1 (MQ=255)
                                                                         cGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAtt                          <  2:1164089/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  >  2:1332538/1‑100 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  >  1:3489838/1‑100 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  <  1:316306/100‑1 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  >  1:723298/1‑100 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTgg  <  1:1383054/100‑1 (MQ=255)
                                                                                                 gggAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTAAGCGCGATCTCATTg   >  1:3127334/1‑99 (MQ=255)
                                                                                                 |                                                                                                   
TACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATCATTGG  >  NC_000913/4470860‑4471056

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: