Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,376,822 A→T 15.9% intergenic (+123/‑53) ecnB → / → sugE entericidin B membrane lipoprotein/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,376,8220AT15.9% 451.2 / 73.6 185intergenic (+123/‑53)ecnB/sugEentericidin B membrane lipoprotein/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (79/75);  new base T (13/16);  total (92/93)
Fisher's exact test for biased strand distribution p-value = 5.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.31e-01

Only 100 of 195 total aligned reads displayed.

GTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATG  >  NC_000913/4376725‑4376921
                                                                                                 |                                                                                                   
gTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGAcc                                                                                                   <  1:2187547/100‑1 (MQ=255)
gTGCCGTTTGTTTCTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGAcc                                                                                                   <  2:4179182/100‑1 (MQ=255)
 tGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCt                                                                                                  <  1:1137908/100‑1 (MQ=255)
  gCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTg                                                                                                 >  1:2668364/1‑100 (MQ=255)
    cGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCa                                                                                               <  2:69377/100‑1 (MQ=255)
       ttgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaag                                                                                            <  1:2122004/100‑2 (MQ=255)
        tgtttTTTCTTCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACg                                                                                           <  2:4009501/100‑1 (MQ=255)
        tgtttTTTCTGCGATCGTCACAAAGGTAATAGTTGAAATTCACCTGCCACGTGGCAAAATATCTGTTCAACCATCATCTTTGCATGTCGGCCTGCaaagc                                                                                           <  2:3922580/100‑3 (MQ=255)
        tgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagc                                                                                           <  2:817686/100‑3 (MQ=255)
        tgtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACg                                                                                           >  1:3148920/1‑100 (MQ=255)
         gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagct                                                                                          <  2:2312653/100‑4 (MQ=255)
         gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGc                                                                                          >  1:727955/1‑100 (MQ=255)
         gtttTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGc                                                                                          >  1:4088401/1‑100 (MQ=255)
          ttttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                                                                         >  2:1526819/1‑100 (MQ=255)
          ttttttCTGAGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGcc                                                                                         <  1:3421542/100‑1 (MQ=255)
           tttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagctga                                                                                        <  2:1708131/100‑6 (MQ=255)
           tttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCct                                                                                        >  1:3377146/1‑100 (MQ=255)
            ttttCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagctgat                                                                                       <  2:1817336/100‑7 (MQ=39)
              ttCTGCGATAGTTACAAAGGTAATAGTTTAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagctgatgg                                                                                     <  2:3731156/100‑9 (MQ=38)
                cTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGTCGTCCTGCaaagctgatggtt                                                                                   <  2:2351700/100‑11 (MQ=38)
                  gCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCa                                                                                 >  1:529768/1‑100 (MQ=255)
                  gCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCa                                                                                 >  1:3529761/1‑100 (MQ=255)
                   cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCCGGCCGACCTGCAAAAGCCTCTTTTCa                                                                                 >  1:3500455/1‑99 (MQ=255)
                   cGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCAc                                                                                >  1:1494251/1‑100 (MQ=255)
                     aTAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCg                                                                              >  1:1048999/1‑100 (MQ=255)
                      tAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgg                                                                             <  1:1263648/100‑1 (MQ=255)
                        gTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg                                                                           <  2:2672515/100‑1 (MQ=255)
                        gTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg                                                                           >  1:2250469/1‑100 (MQ=255)
                        gTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCgggg                                                                           >  1:2673487/1‑100 (MQ=255)
                          cacaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGAc                                                                         >  2:2134752/1‑100 (MQ=255)
                           acaAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACg                                                                        >  1:3827562/1‑100 (MQ=255)
                            caAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACgg                                                                       >  2:156879/1‑100 (MQ=255)
                             aaaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGc                                                                      <  1:2714007/100‑1 (MQ=255)
                              aaGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGCTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGcc                                                                     <  1:2628902/100‑1 (MQ=255)
                                   aaTAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCagt                                                                >  2:2545635/1‑98 (MQ=255)
                                   aaTAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAt                                                                <  1:3102268/100‑1 (MQ=255)
                                   aaTAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAAt                                                                <  2:2909836/100‑1 (MQ=255)
                                     tAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtc                                                              <  2:169487/100‑1 (MQ=255)
                                      aGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtct                                                             >  1:3621016/1‑100 (MQ=255)
                                      aGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtct                                                             <  1:448037/100‑1 (MQ=255)
                                      aGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtct                                                             >  1:1300071/1‑100 (MQ=255)
                                       gTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACTACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtctc                                                            <  1:3018409/100‑1 (MQ=255)
                                       gTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATtctc                                                            >  2:1612057/1‑100 (MQ=255)
                                        ttGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTcc                                                           >  2:2069380/1‑100 (MQ=255)
                                         gatctttCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCg                                                          <  2:1267819/95‑1 (MQ=255)
                                            aaTTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAg                                                       >  2:334266/1‑100 (MQ=255)
                                                 cccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGa                                                  >  1:2020592/1‑100 (MQ=255)
                                                  ccTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCAGat                                                 >  2:3528190/1‑100 (MQ=255)
                                                        aCCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTcc                                           <  2:4148021/100‑1 (MQ=255)
                                                         ccTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCTCCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTcc                                           <  2:3003904/99‑1 (MQ=255)
                                                           tGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTgg                                        >  2:3444892/1‑100 (MQ=255)
                                                             gCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAt                                      <  1:3887722/100‑1 (MQ=255)
                                                             gCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAt                                      <  1:2006031/100‑1 (MQ=255)
                                                             gCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAt                                      >  1:4035562/1‑100 (MQ=255)
                                                              cAAGATAGCCGTTCAACCATCAGCGCTGCAGGCCGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTATCCGGAGCCTGATATGTCCTGGAt                                      <  2:3659224/99‑1 (MQ=255)
                                                              cAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAtt                                     >  2:748551/1‑100 (MQ=255)
                                                              cAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGAtt                                     <  1:739363/100‑1 (MQ=255)
                                                                 aaTATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATc                                  <  1:1728408/100‑1 (MQ=255)
                                                                  atatCCGTGCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCt                                 <  2:1775507/100‑1 (MQ=255)
                                                                    atCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAGCGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                               >  2:2229743/1‑100 (MQ=255)
                                                                    atCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTa                               <  2:1951672/100‑1 (MQ=255)
                                                                     tCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAg                              <  2:1764843/100‑1 (MQ=255)
                                                                     tCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAg                              <  2:578068/100‑1 (MQ=255)
                                                                       cGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGtt                            >  1:4154433/1‑100 (MQ=255)
                                                                       cGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGtt                            >  1:1098922/1‑100 (MQ=255)
                                                                       cGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGtt                            >  2:3987916/1‑100 (MQ=255)
                                                                        gTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTa                           >  2:407600/1‑100 (MQ=255)
                                                                         ttCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTAt                          <  2:2794656/100‑1 (MQ=255)
                                                                           cAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTg                        >  2:870622/1‑100 (MQ=255)
                                                                            aaCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                       >  2:3146225/1‑100 (MQ=255)
                                                                            aaCCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGc                       <  1:2401449/100‑1 (MQ=255)
                                                                               cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTgg                    <  1:3728450/100‑1 (MQ=255)
                                                                               cATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTgg                    <  2:3274015/100‑1 (MQ=255)
                                                                                 tCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTc                  <  1:2738137/100‑1 (MQ=255)
                                                                                  ggcgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTct                 >  2:187389/5‑100 (MQ=255)
                                                                                  ggcgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTct                 >  1:3733664/5‑100 (MQ=255)
                                                                                    gCTTTGCAGGACGACCTGCAAACGTCTCTTTTCACTGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               <  1:396447/100‑1 (MQ=255)
                                                                                    cgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               >  1:1370788/3‑100 (MQ=255)
                                                                                    cgtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgc               >  1:3094025/3‑100 (MQ=255)
                                                                                     gtttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctg                >  2:918878/2‑98 (MQ=255)
                                                                                      tttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTctgctg             <  2:33767/100‑1 (MQ=255)
                                                                                       ttGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTgg            <  1:187389/100‑1 (MQ=255)
                                                                                       ttGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGCTATGTCCTGGATTATCTTAGTTATTGCTGGTctgctga            <  2:3134167/100‑2 (MQ=255)
                                                                                        tGCAGGTCGTCCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGa           >  1:2854188/1‑100 (MQ=255)
                                                                                        tGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGa           <  2:735060/100‑1 (MQ=255)
                                                                                         gCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa          >  1:930866/1‑100 (MQ=255)
                                                                                         gCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGaa          >  1:887672/1‑100 (MQ=255)
                                                                                          cAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAg         >  2:3120570/1‑100 (MQ=255)
                                                                                          cAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAg         >  1:1542737/1‑100 (MQ=255)
                                                                                           aGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGt        <  1:4172843/100‑1 (MQ=255)
                                                                                            ggACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg       >  1:3596333/1‑100 (MQ=255)
                                                                                            ggACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg       >  1:2355088/1‑100 (MQ=255)
                                                                                            ggACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg       >  1:1465181/1‑100 (MQ=255)
                                                                                            ggACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTg       <  1:1050681/100‑1 (MQ=255)
                                                                                               cgacCTGCAAACGCTTCTTTTCACCGTGGACGGTCCCAATTCTCCGGAGCCTGTTATGTCCTGGATTATCGTATTTATTGCTGGTCTGCTGGAAGTGGTa    >  2:1659273/1‑100 (MQ=255)
                                                                                               cgacCTGCAAACGCCTCTTTTCACCGGGGATGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTa    >  2:3792834/1‑100 (MQ=255)
                                                                                               cgacCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCGGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTgg      <  2:3572010/98‑1 (MQ=255)
                                                                                                 acCTGCAACCGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATg  <  2:1050727/100‑1 (MQ=255)
                                                                                                 acCTGCAAACGCCTCTTTTCACCGGGGACGGCCTCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATg  >  2:3038326/1‑100 (MQ=255)
                                                                                                 acCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATg  <  1:613902/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GTGCCGTTTGTTTTTTCTGCGATAGTCACAAAGGTAATAGTTGAAATTCCCCTGCCACCTGGCAAAATATCCGTTCAACCATCAGCTTTGCAGGACGACCTGCAAACGCCTCTTTTCACCGGGGACGGCCCCAATTCTCCGGAGCCTGATATGTCCTGGATTATCTTAGTTATTGCTGGTCTGCTGGAAGTGGTATG  >  NC_000913/4376725‑4376921

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