Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,470,951 A→C 15.6% intergenic (‑38/+35) 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,9510AC15.6% 523.5 / 82.8 210intergenic (‑38/+35)ridA/pyrIenamine/imine deaminase, reaction intermediate detoxification/aspartate carbamoyltransferase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base A (90/87);  new base C (19/14);  total (109/101)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

Only 100 of 212 total aligned reads displayed.

ACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATAT  >  NC_000913/4470852‑4471050
                                                                                                   |                                                                                                   
 cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAt                                                                                                    <  2:2429470/100‑1 (MQ=255)
 cGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAt                                                                                                    <  1:4283904/100‑1 (MQ=255)
  gCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCAtt                                                                                                   >  2:4323663/1‑100 (MQ=255)
      tGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAg                                                                                               <  1:4091992/100‑1 (MQ=255)
       gTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAgg                                                                                              <  2:3393507/100‑1 (MQ=255)
        tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggg                                                                                             >  1:3854/1‑100 (MQ=255)
        tACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAggg                                                                                             >  1:3563687/1‑100 (MQ=255)
          cGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAg                                                                                           <  2:3262391/100‑1 (MQ=255)
           gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAGAGTCTGGCTCATTATAGGGAg                                                                                           >  1:1027362/1‑99 (MQ=255)
           gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGc                                                                                          <  2:4026931/100‑1 (MQ=255)
           gTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGc                                                                                          >  1:2123135/1‑100 (MQ=255)
            tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTCAAAAGACTGGCTCATTATAgg                                                                                              >  1:668706/1‑95 (MQ=255)
            tAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGcc                                                                                         <  2:2209629/100‑1 (MQ=255)
              aGGACCGCTACCTTCCGGTTCATTTTCCGTCGCCATAGGTTTGCTCATGATTTCTCCTTTATTACAGCGGTATAAAGTCTGACTCATTGTTGAGAGc                                                                                          >  1:783374/1‑97 (MQ=255)
              aGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAg                                                                                       >  2:1987223/1‑100 (MQ=255)
               ggACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGa                                                                                      >  2:907416/1‑100 (MQ=255)
                gACCGATAGCTGCCGTTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  2:4331315/100‑1 (MQ=37)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  2:133828/100‑1 (MQ=38)
                gACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGAc                                                                                     <  1:2436606/100‑1 (MQ=38)
                 aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACt                                                                                    <  2:3980892/100‑1 (MQ=38)
                 aCCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACt                                                                                    <  1:2739939/100‑1 (MQ=255)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtt                                                                                   <  1:2862319/100‑1 (MQ=38)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                   >  2:2557021/1‑100 (MQ=255)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                   >  2:1666506/1‑100 (MQ=255)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGGAAAAAGTCTGGCTCATGATAGGGAGCCAGACt                                                                                    >  2:3313632/1‑99 (MQ=255)
                  ccGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTATGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtt                                                                                   >  2:3107677/1‑100 (MQ=255)
                   cGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACttt                                                                                  <  1:2986269/100‑1 (MQ=37)
                   cGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACttt                                                                                  <  2:3885748/100‑1 (MQ=37)
                    gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                 >  1:4226761/1‑100 (MQ=255)
                    gATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACtttt                                                                                 <  2:1318066/100‑1 (MQ=255)
                       aGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATAATGAGCCAGACtttttac                                                                              <  1:1960552/100‑4 (MQ=37)
                        gCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttatta                                                                             >  2:489210/1‑100 (MQ=255)
                         cTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTttattan                                                                            <  1:523722/100‑2 (MQ=255)
                          tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                           >  2:888388/1‑100 (MQ=255)
                          tGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTAcc                                                                           >  1:1149335/1‑100 (MQ=255)
                            ccGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTGATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCaa                                                                         <  1:614222/100‑1 (MQ=255)
                              ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCCCTATCATGAGCCAGACtttttaccgctgta                                                                       >  2:668708/1‑90 (MQ=21)
                              ggTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                       >  2:1836022/1‑100 (MQ=255)
                                tGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGc                                                                     >  2:1471979/1‑100 (MQ=255)
                                 gCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCa                                                                    >  2:3330527/1‑100 (MQ=255)
                                 gCATTTTCCGTCGCGATAGTTTTGCTCATGATTGCTCCTTTATTACAGCGGTAAAAAGTCTGCCTCATTATAGGGAGCCAGACTTTTATTACCAAcc                                                                       >  1:2598770/1‑97 (MQ=255)
                                   aTTTTCCGTCGCGATAGTTTTGCTCATGATATCACCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTATCAACCGCaat                                                                  >  1:2890031/1‑100 (MQ=255)
                                         cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTGATTATAGGGATCCAGACTTTTATTACCATCCGCAATTAATTg                                                            <  2:2967704/100‑1 (MQ=255)
                                         cGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTg                                                            >  2:302495/1‑100 (MQ=255)
                                            cgcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                         <  2:2372959/100‑1 (MQ=255)
                                            cgcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGcc                                                         <  2:2811150/100‑1 (MQ=255)
                                            cgcgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACAACCGCAATTAATTGGCCa                                                        <  1:1165083/100‑1 (MQ=255)
                                              cgATAGTTTTGCTCATGATTTCTCCTTTATTACATCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                       <  1:1587848/100‑1 (MQ=255)
                                              cgATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAg                                                       <  2:1571852/100‑1 (MQ=255)
                                               gATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGc                                                      >  1:3662970/1‑100 (MQ=255)
                                                   gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                  >  2:1217338/1‑100 (MQ=255)
                                                   gTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcacca                                                  >  2:589896/1‑100 (MQ=255)
                                                     tttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGcaccac                                                 >  1:3440687/1‑99 (MQ=255)
                                                      ttGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAt                                               <  1:3835865/100‑1 (MQ=255)
                                                       tGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACAtt                                              >  1:962625/1‑100 (MQ=255)
                                                        gCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTa                                             <  2:3150566/100‑1 (MQ=255)
                                                         cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                            >  1:1856806/1‑100 (MQ=255)
                                                         cTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTAt                                            >  2:3841302/1‑100 (MQ=255)
                                                          tCATTATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATg                                           <  2:1856812/100‑1 (MQ=255)
                                                           cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                          >  2:4270325/1‑100 (MQ=255)
                                                           cATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                          >  1:2364587/1‑100 (MQ=255)
                                                           cATGATTTCTCCTTGATTACAGCGGTAAAAAGTCTGGCTCATTATATGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATgg                                          <  1:1602878/100‑1 (MQ=255)
                                                            aTGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATggg                                         <  2:2501180/100‑1 (MQ=255)
                                                             tGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGa                                        <  2:4011508/100‑1 (MQ=255)
                                                               aTTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGaaa                                      >  2:653418/1‑100 (MQ=255)
                                                                 ttCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAAc                                    <  2:2252590/100‑1 (MQ=255)
                                                                   ctcCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAActc                                  >  1:4324086/1‑100 (MQ=255)
                                                                     ccTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtt                                >  1:4203984/1‑100 (MQ=255)
                                                                       tttATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCtttt                              >  1:1090767/1‑100 (MQ=255)
                                                                           ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                          <  2:4304369/100‑1 (MQ=255)
                                                                           ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                          <  1:1264840/100‑1 (MQ=255)
                                                                           ttaCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcac                          >  1:3435710/1‑100 (MQ=255)
                                                                            taCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                         <  1:59186/100‑1 (MQ=255)
                                                                            taCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                         <  1:3886543/100‑1 (MQ=255)
                                                                            taCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTcaca                         >  1:306072/1‑100 (MQ=255)
                                                                             gttggtaataaaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                        >  2:2421135/10‑100 (MQ=21)
                                                                             aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                        >  2:790464/1‑100 (MQ=255)
                                                                             aCAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAg                        >  2:1730737/1‑100 (MQ=255)
                                                                              cAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGt                       <  2:3425566/100‑1 (MQ=255)
                                                                                gCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTAt                     <  1:3811867/100‑1 (MQ=255)
                                                                                   gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                  <  2:2883672/100‑1 (MQ=255)
                                                                                   gTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTg                  <  2:2656637/100‑1 (MQ=255)
                                                                                    tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                 <  2:3493535/100‑1 (MQ=255)
                                                                                    tAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGc                 <  2:3703759/100‑1 (MQ=255)
                                                                                     aaaaaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCa                <  2:3615249/100‑1 (MQ=255)
                                                                                       aaaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTAACAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAtt              >  1:2338247/1‑100 (MQ=37)
                                                                                        aaGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAACTAATTGGCAAGCCCCACATAATAGGCAAACTCTTTTTCACAGTATTTGCAtt              >  1:36661/1‑99 (MQ=21)
                                                                                        aaGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt             >  2:2093004/1‑100 (MQ=255)
                                                                                         aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCAttt             >  1:2083294/1‑99 (MQ=38)
                                                                                         aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTg            >  1:2855697/1‑100 (MQ=38)
                                                                                         aGTCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTACATTTg            >  1:573023/1‑100 (MQ=37)
                                                                                           tCTGGCTCCCTATAATGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg          >  1:269506/1‑100 (MQ=38)
                                                                                            cTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgc         >  2:2902702/1‑100 (MQ=255)
                                                                                             tGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTACTTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcc        <  2:1027366/100‑2 (MQ=255)
                                                                                              ggCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAgcgc       <  2:1721076/100‑1 (MQ=255)
                                                                                               gCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGcgcg      >  2:4194785/1‑100 (MQ=255)
                                                                                               gCTCATTATAGGGAGCCACACTTTTATTACCAACAGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAg          <  1:1994485/96‑1 (MQ=255)
                                                                                                cTCATTATAGGGAGCCAGACTTTTATTACAAACCACAATTAATAAAAAAGCAGCACAGTATGGGAAAACGCTTGTTCACAGTATTTGCGTGTGAGCGCGa     >  1:1458169/1‑100 (MQ=255)
                                                                                                  cATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGata   >  1:1338005/1‑100 (MQ=255)
                                                                                                   aTTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGatat  <  2:2927483/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
ACGCCCTGTACGTAAGGACCGATAGCTGCCGGTGCATTTTCCGTCGCGATAGTTTTGCTCATGATTTCTCCTTTATTACAGCGGTAAAAAGTCTGGCTCATTATAGGGAGCCAGACTTTTATTACCAACCGCAATTAATTGGCCAGCACCACATTATGGGAAAACTCTTTTTCACAGTATTTGCATTTGAGCGCGATAT  >  NC_000913/4470852‑4471050

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: