Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,952 T→C 15.0% intergenic (‑39/+34) 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,9520TC15.0% 416.8 / 63.0 182intergenic (‑39/+34)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base T (77/76);  new base C (9/18);  total (88/94)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

Only 100 of 184 total aligned reads displayed.

CCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATC  >  NC_000913/4470855‑4471051
                                                                                                 |                                                                                                   
cccTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                    <  2:3509097/99‑1 (MQ=255)
cccTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATta                                                                                                   >  2:468676/1‑100 (MQ=255)
 ccTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATtat                                                                                                  <  1:438069/100‑1 (MQ=255)
 ccTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATtat                                                                                                  <  1:1855415/100‑1 (MQ=255)
      aCGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGa                                                                                             <  2:2997239/100‑1 (MQ=255)
        gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGc                                                                                           >  1:973954/1‑100 (MQ=255)
           aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAg                                                                                        <  2:240326/100‑1 (MQ=38)
            ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                       <  2:3045096/100‑1 (MQ=38)
            ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGa                                                                                       <  1:1442708/100‑1 (MQ=38)
             gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                      <  1:3707759/100‑1 (MQ=38)
             gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                      <  1:3657616/100‑1 (MQ=38)
             gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                      <  2:2151547/100‑1 (MQ=38)
             gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGAc                                                                                      >  1:2175389/1‑100 (MQ=255)
              aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                     <  1:1045233/100‑1 (MQ=255)
               ccGATAGCTGCCGGTGCTTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                    <  2:3811799/100‑1 (MQ=38)
               ccGATAGCTGCCGGTGCATTTTCCGTCGCTATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                    <  1:3498641/100‑1 (MQ=255)
                cGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACttn                                                                                   <  1:266804/100‑2 (MQ=37)
                 gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGTAGCCAGACtttt                                                                                  >  2:3328665/1‑100 (MQ=255)
                 gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                  <  1:2747102/100‑1 (MQ=255)
                  aTAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTtta                                                                                 <  1:1883847/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)
                     gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                              >  2:2731298/1‑100 (MQ=255)
                      cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAc                                                                             >  1:64211/1‑100 (MQ=255)
                        gCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                            >  2:2872381/1‑99 (MQ=255)
                          cGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                          <  2:113416/99‑1 (MQ=255)
                           ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                        >  1:512714/1‑100 (MQ=255)
                           ggTGCATTTTCCGTCGCGATAGTTTTGCTAATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                        <  2:1427823/100‑1 (MQ=255)
                            gTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCg                                                                       <  2:695417/100‑1 (MQ=255)
                             tGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                      <  2:1154624/100‑1 (MQ=255)
                              gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCATTTAAAACGGCGGTAAAAAGTCTCGCTCATTATAGGGAGCCAGACTTTTATTACCa                                                                           <  2:3384591/94‑1 (MQ=255)
                               cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaa                                                                    >  2:809519/1‑100 (MQ=255)
                               cATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTAAAGCGGTAAAAAGTCTGGCTCATTATAGGGTGCCAGACTTTTATTACCAACCGCaa                                                                    >  2:1718869/1‑100 (MQ=255)
                                aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGTGATCCAGACTTTTATTACCAACCGCaat                                                                   >  2:178207/1‑100 (MQ=255)
                                aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                   <  2:2011529/100‑1 (MQ=255)
                                aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                   >  2:3515997/1‑100 (MQ=255)
                                aTTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTAATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                   >  1:659685/1‑100 (MQ=255)
                                aTTTTCCGTCGCGATAGATTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaat                                                                   >  1:1893822/1‑100 (MQ=255)
                                 ttttCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaatt                                                                  >  2:3282347/1‑100 (MQ=255)
                                    tCCGCCGAGATAGTGTTGCTCTTGATTGCTCCTTTATTACAGCGTTAAAACGTTTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAACGCtatgaat                                                               <  2:3593602/100‑1 (MQ=255)
                                      cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                              >  2:2364646/1‑99 (MQ=255)
                                      cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTg                                                             >  1:4047987/1‑100 (MQ=255)
                                      cGTCGCGATAGTTTTGCTCATGATTTCTCCTTGATTACATCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCaattaatt                                                              <  2:974682/99‑1 (MQ=255)
                                       gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                            >  1:1728274/1‑100 (MQ=255)
                                       gTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTgg                                                            <  2:783201/100‑1 (MQ=255)
                                          gcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCa                                                         <  1:3582476/100‑1 (MQ=255)
                                            gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                        >  2:2941748/1‑99 (MQ=255)
                                              tAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                     <  1:2275676/100‑1 (MQ=255)
                                               aGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacc                                                    >  1:1604550/1‑100 (MQ=255)
                                               aGTTTTGCTCATGATTTCTCCTTTAGTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcac                                                     <  2:1430570/99‑1 (MQ=255)
                                                 ttttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                  <  2:2681244/100‑1 (MQ=255)
                                                   ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                <  2:3134299/100‑1 (MQ=255)
                                                   ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                <  1:467762/100‑1 (MQ=255)
                                                   ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                                <  1:1555351/100‑1 (MQ=255)
                                                        cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                           <  1:468827/100‑1 (MQ=255)
                                                        cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                           <  2:4008684/100‑1 (MQ=255)
                                                        cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                           <  1:766305/100‑1 (MQ=255)
                                                          tGATTTCTCCTTTTTTACAGCGGTCAAAAGTCTGGCTCATTATAGGGGGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                          <  2:1708828/99‑1 (MQ=255)
                                                          tGATTTCTCCTTTATTACAACGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                         >  1:1427843/1‑100 (MQ=255)
                                                          tGATTTCTCCTTTATTAAAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTATATTACCAACCGCAATTAATTGGCAAGCAGCACATTATGGGa                                         >  2:371378/1‑100 (MQ=255)
                                                             tttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaaa                                      >  1:1321590/1‑100 (MQ=255)
                                                              ttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                     <  2:2335719/100‑1 (MQ=255)
                                                               tctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAct                                    <  1:168174/100‑1 (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:935174/1‑100 (MQ=255)
                                                                      tattaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTc                             <  1:3238596/100‑1 (MQ=255)
                                                                       attaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTca                            >  2:1588840/1‑100 (MQ=255)
                                                                        ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                           <  1:2884566/100‑1 (MQ=255)
                                                                        ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                           <  2:378493/100‑1 (MQ=255)
                                                                        ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                           <  2:1238329/100‑1 (MQ=255)
                                                                        ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                           >  2:1571663/1‑100 (MQ=255)
                                                                         ggttggtaataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                          >  1:198679/11‑100 (MQ=21)
                                                                           cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcacatt                        >  2:553452/1‑98 (MQ=255)
                                                                           cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                         >  2:211998/1‑99 (MQ=255)
                                                                           cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGt                        >  1:4002794/1‑100 (MQ=255)
                                                                             gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                      <  2:2633378/100‑1 (MQ=255)
                                                                               ggTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAttt                    <  1:584183/100‑1 (MQ=255)
                                                                               ggTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAttt                    >  1:3390276/1‑100 (MQ=255)
                                                                                gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAATa                       >  1:1295118/1‑96 (MQ=255)
                                                                                gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGCAAAACACCTTTCCACAGTAt                      >  1:1787819/1‑97 (MQ=255)
                                                                                aataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                   >  1:3966114/4‑100 (MQ=37)
                                                                                    aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt               >  1:2794631/1‑100 (MQ=37)
                                                                                    aaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt               <  1:1579539/100‑1 (MQ=255)
                                                                                     atgTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt              <  2:2499947/98‑1 (MQ=255)
                                                                                     aaGTCTGGCTCAGTATAGGGAGCCAGACTTTTCTTACCAACCGCAATTAATTGGCCAGCACCGCCTTATGGGAAAACTATTTTTCACAGCATCTGCAtt               >  2:2810701/1‑99 (MQ=255)
                                                                                      aGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg             <  1:3901033/100‑1 (MQ=255)
                                                                                        tCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg           >  1:2461960/1‑100 (MQ=38)
                                                                                         ctgcgccATTATAGGGAGCCAGACTTTTGTTACCAAACGCAATTTATTGGTCAGGACCCCACTAAGCGACAGCGGTTTTTCCCAGTATTTGCATTTGa            >  2:1775791/7‑98 (MQ=255)
                                                                                         cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc          >  1:1978292/1‑100 (MQ=255)
                                                                                         cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATGATGGGAAAACTCTTTTTCACAGTATTTGAATTTGAgc          >  1:3439753/1‑100 (MQ=255)
                                                                                           ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc        >  1:3189277/1‑100 (MQ=255)
                                                                                            gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg       <  1:1483023/100‑1 (MQ=255)
                                                                                             cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa      >  1:3837472/1‑100 (MQ=255)
                                                                                             cTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCACCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGa      <  2:361342/100‑1 (MQ=255)
                                                                                              tCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg       <  2:1262887/98‑1 (MQ=255)
                                                                                                 |                                                                                                   
CCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATATC  >  NC_000913/4470855‑4471051

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