Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,172,439 A→T 7.6% intergenic (‑228/‑91) ygiZ ← / → mdaB conserved inner membrane protein YgiZ/NADPH:quinone oxidoreductase MdaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,172,4390AT 7.6% 373.2 / 15.0 119intergenic (‑228/‑91)ygiZ/mdaBconserved inner membrane protein YgiZ/NADPH:quinone oxidoreductase MdaB
Reads supporting (aligned to +/- strand):  ref base A (47/63);  new base T (6/3);  total (53/66)
Fisher's exact test for biased strand distribution p-value = 1.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 124 total aligned reads displayed.

CATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTG  >  NC_000913/3172291‑3172589
                                                                                                                                                    |                                                                                                                                                      
cATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAAt                                                                                                                                                       <  1:872506/150‑1 (MQ=255)
  tatTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTc                                                                                                                                                     <  1:1253874/150‑1 (MQ=255)
   atTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCg                                                                                                                                                    <  2:45611/150‑1 (MQ=255)
    ttttCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGc                                                                                                                                                   <  2:438022/150‑1 (MQ=255)
    ttttCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGc                                                                                                                                                   >  1:1426017/1‑150 (MQ=255)
     tttCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCa                                                                                                                                                  <  2:806138/150‑1 (MQ=255)
        cAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCaaaa                                                                                                                                               <  2:752970/150‑1 (MQ=255)
         aaTCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAg                                                                                                                                              <  1:371313/150‑1 (MQ=255)
           tCCCCAAAAAAACCGAAGCAAACCTTTAAATTATTCCGGTGCGACTAATAATTATTATTAATGGCTATGAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCTGTGATTCATAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGAc                                                                                                                                            <  2:532268/150‑1 (MQ=255)
           tCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCGAATTTAGCAAAAgtc                                                                                                                                            <  2:236820/150‑3 (MQ=255)
            ccccAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCg                                                                                                                                                    <  2:934803/141‑1 (MQ=255)
            ccccAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCg                                                                                                                                                    >  1:934803/1‑141 (MQ=255)
             cccAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTTCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGAAtt                                                                                                                                          >  1:146430/1‑150 (MQ=255)
             cccAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACtt                                                                                                                                          <  2:938200/150‑1 (MQ=255)
             cccAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGCAAGATCAGACTTTTGCTAAATTCGCAAAAGACtt                                                                                                                                          <  1:635661/150‑1 (MQ=255)
             cccAAAAAAAACGAACCAAACGTGTAAATTAATCCGGTGATAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACtt                                                                                                                                          >  2:763457/1‑150 (MQ=255)
                aaaaaaaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGc                                                                                                                                       <  2:1186917/150‑1 (MQ=255)
                  aaaaaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGcac                                                                                                                                     >  1:210693/1‑150 (MQ=255)
                     aaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACAtt                                                                                                                                  >  1:1289327/1‑150 (MQ=255)
                     aaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACAtt                                                                                                                                  >  1:392659/1‑150 (MQ=255)
                      aCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACAttt                                                                                                                                 <  2:1129679/150‑1 (MQ=255)
                       ccGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACAtttt                                                                                                                                <  1:1068488/150‑1 (MQ=255)
                       ccGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCGAATTTAGCAAAAGTCtgatc                                                                                                                                       >  2:655464/1‑139 (MQ=255)
                       ccGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCGAATTTAGCAAAAGTCtgatc                                                                                                                                       <  1:655464/143‑5 (MQ=255)
                           aCCAAACCTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTa                                                                                                                            >  2:475264/1‑150 (MQ=255)
                             cAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGCAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAAt                                                                                                                          <  2:1615338/150‑1 (MQ=255)
                               aaCGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAAttt                                                                                                                        >  1:979066/1‑150 (MQ=255)
                                    atAAATTATTCCGCTGCGAGTAAGAATTATGAGGAATGGATATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGACTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGAACATTTTGCTAATTTAACCg                                                                                                                   >  1:54175/2‑150 (MQ=255)
                                        attattCCGGTGCGAGTATGACTTTTGAGGAATGGCTATCAGTATTGTCATTTTCAGACAATATTTATCCTGCATCGGTGAGGCAGAGTAAGCTGAGACTTTTGCTAAATTCGGAAATGACTTTGCACATTTTGCTAATTTCACCGTAcc                                                                                                               <  2:506188/150‑1 (MQ=255)
                                        attattCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTAcc                                                                                                               <  1:1103098/150‑1 (MQ=255)
                                          tattCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGc                                                                                                             <  1:370664/150‑1 (MQ=255)
                                            ttCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGctc                                                                                                           <  1:915555/150‑1 (MQ=255)
                                            ttCCGGTGCGAGTAAGAATTATGAGGAATGGATATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGctc                                                                                                           >  2:443845/1‑150 (MQ=255)
                                                ggTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCtgtg                                                                                                       >  1:1603338/1‑150 (MQ=255)
                                                       gTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACt                                                                                                >  2:1309585/1‑150 (MQ=255)
                                                         aaGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGATAAATTATCAAAATACTTTGAAAATTTTTCTAATTTAAACTTACCGCTATGTGAAGTACtat                                                                                              <  2:126265/150‑1 (MQ=255)
                                                          aGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACtata                                                                                             <  1:1649585/150‑1 (MQ=255)
                                                          aGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACtat                                                                                              >  1:555335/1‑149 (MQ=255)
                                                          aGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACtat                                                                                              <  2:555335/149‑1 (MQ=255)
                                                           taaTTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAg                                                                                            <  1:1591022/149‑1 (MQ=255)
                                                           gAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAg                                                                                            <  1:945394/150‑1 (MQ=255)
                                                           gAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAg                                                                                            >  2:1111923/1‑150 (MQ=255)
                                                             aTTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTc                                                                                          >  2:1253578/1‑150 (MQ=255)
                                                                 tGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGt                                                                                                                  <  1:19096/122‑1 (MQ=255)
                                                                 tGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCa                                                                                      >  1:491380/1‑150 (MQ=255)
                                                                  gAGGAATGGCTATGAGTATTGTCATTTTCAGAAAATTTTTATCCTGCATCGTTGAGTCATAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCaa                                                                                     <  2:1049457/150‑1 (MQ=255)
                                                                     gAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACg                                                                                  <  2:1399257/150‑1 (MQ=255)
                                                                     gAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACg                                                                                  <  2:939949/150‑1 (MQ=255)
                                                                             aTCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTAc                                                                                                 <  1:1135330/127‑1 (MQ=255)
                                                                              tCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCtt                                                                         <  2:194734/150‑1 (MQ=255)
                                                                                tgTATTGTCATTTTCAGAAAATATTTATACTTCATCGTTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCTCAAAAGACTTTGCTCATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGAAAACGTCTCACCTTGa                                                                       <  2:221809/149‑1 (MQ=255)
                                                                                 gTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAg                                                                      >  2:797539/1‑150 (MQ=255)
                                                                                    ttGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGtt                                                                   >  1:1494691/1‑150 (MQ=255)
                                                                                      gTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCAcc                                                                           <  2:888724/140‑1 (MQ=255)
                                                                                      gTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTaa                                                                 >  2:1207862/1‑150 (MQ=255)
                                                                                        cATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTaaaa                                                               <  2:271173/150‑1 (MQ=255)
                                                                                         aTTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTaaaaa                                                              >  1:1319611/1‑150 (MQ=255)
                                                                                            ttCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATg                                                           <  2:283250/150‑1 (MQ=255)
                                                                                             tCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGa                                                          <  2:432920/150‑1 (MQ=255)
                                                                                              cAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAg                                                         <  1:663524/150‑1 (MQ=255)
                                                                                              cAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAg                                                         <  1:1337085/150‑1 (MQ=255)
                                                                                              cAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAg                                                         >  2:1601441/1‑150 (MQ=255)
                                                                                              cAGAAAATATTTATCCTGCATCGGTCAGTCAGAGTAAGATCAGACTTTTGCTAAATTAGCAAAAGACTTTGCACATTTTGCTAAATTCACCTTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACATTGAGGTTAAAAAATGac                                                         >  1:1640977/1‑149 (MQ=255)
                                                                                               aGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGc                                                        <  2:1321199/150‑1 (MQ=255)
                                                                                                         tATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGa                                              >  2:216602/1‑150 (MQ=255)
                                                                                                          aTCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGAt                                             <  1:1253578/150‑1 (MQ=255)
                                                                                                          aTCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGAt                                             <  1:700943/150‑1 (MQ=255)
                                                                                                             cTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTAt                                          >  2:341343/1‑150 (MQ=255)
                                                                                                               gCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCa                                        >  2:1220618/1‑150 (MQ=255)
                                                                                                                    ggTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGgc                                   >  1:1483075/1‑150 (MQ=255)
                                                                                                                      tGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGgcgc                                 <  1:100641/150‑1 (MQ=255)
                                                                                                                      tGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGgcgc                                 <  2:1289327/150‑1 (MQ=255)
                                                                                                                          tCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCCCCGTCCGGGTGTGTGACGTCCTATAGTCGGCAAACGTGTCACGTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGaaa                             >  1:757512/1‑150 (MQ=255)
                                                                                                                             gagTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGaaaaaa                          >  2:220597/1‑150 (MQ=255)
                                                                                                                               gTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAAtt                        >  2:538097/1‑150 (MQ=255)
                                                                                                                               gTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAAtt                        >  2:508421/1‑150 (MQ=255)
                                                                                                                                tAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTc                       <  2:392659/150‑1 (MQ=255)
                                                                                                                                 aaGATCAGACTTTTGCTAAATTCGCAATAGACTTTGCTCATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATATTCTTCAAACGTCTCACCTTGATGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAATTTCg                      <  2:1640977/150‑1 (MQ=255)
                                                                                                                                   gATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGcc                    <  2:630908/150‑1 (MQ=255)
                                                                                                                                    aTCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGccc                   <  2:1572263/150‑1 (MQ=255)
                                                                                                                                     tCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGaaa                             <  2:1667024/139‑1 (MQ=255)
                                                                                                                                     tCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGaaa                             >  1:1667024/1‑139 (MQ=255)
                                                                                                                                      aagtcTTTTGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCAc                 >  2:1025983/5‑150 (MQ=255)
                                                                                                                                      aagtcTTTTGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGAACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGACCAc                 >  2:1455558/5‑150 (MQ=255)
                                                                                                                                       aGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACt                <  1:264245/150‑1 (MQ=255)
                                                                                                                                        gtcTTTTGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTc               >  2:719963/3‑150 (MQ=255)
                                                                                                                                         aCTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTcc              <  2:1503681/150‑1 (MQ=255)
                                                                                                                                          cTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCa             >  2:479177/1‑150 (MQ=255)
                                                                                                                                           ttttGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGGGAAAAAATTCTCCCACTCCaa            <  2:491380/150‑1 (MQ=255)
                                                                                                                                            tttGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAAt           >  1:1531226/1‑150 (MQ=255)
                                                                                                                                            tttGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAAt           >  1:1296072/1‑150 (MQ=255)
                                                                                                                                             ttGCTAAATTCGCAAAAGACTTTGCACATTTTTCTAATTTCACCGTACCGCTCTTTGACGTACTATAGTCCGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTTATTATCAACGGCGCTAAAAAATTCGCCCACTCCAATg          <  1:691051/150‑1 (MQ=255)
                                                                                                                                              tGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATgg         >  2:1656792/1‑150 (MQ=255)
                                                                                                                                              tGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATgg         >  1:617807/1‑150 (MQ=255)
                                                                                                                                                 tAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCa      <  1:1111923/150‑1 (MQ=255)
                                                                                                                                                  aaaTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCaa     >  2:380418/1‑150 (MQ=255)
                                                                                                                                                   aaTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCANc    <  1:425687/150‑1 (MQ=255)
                                                                                                                                                   aaTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAAc    <  1:431643/150‑1 (MQ=255)
                                                                                                                                                   aaTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAAc    >  1:773137/1‑150 (MQ=255)
                                                                                                                                                    aTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCG‑AAAAATTCGCCCACTCCAATGGTCAACTg  <  1:20321/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CATATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTG  >  NC_000913/3172291‑3172589

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑