Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,172,442 C→A 8.5% intergenic (‑231/‑88) 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,4420CA 8.5% 406.5 / 14.2 123intergenic (‑231/‑88)ygiZ/mdaBconserved inner membrane protein YgiZ/NADPH:quinone oxidoreductase MdaB
Reads supporting (aligned to +/- strand):  ref base C (48/64);  new base A (7/4);  total (55/68)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 129 total aligned reads displayed.

TATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAA  >  NC_000913/3172293‑3172591
                                                                                                                                                     |                                                                                                                                                     
tatTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTc                                                                                                                                                       <  1:1253874/150‑1 (MQ=255)
 atTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCg                                                                                                                                                      <  2:45611/150‑1 (MQ=255)
 atTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCg                                                                                                                                                      <  1:474096/150‑1 (MQ=255)
 atTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCg                                                                                                                                                      <  2:844970/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)
   tttCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCa                                                                                                                                                    <  2:233268/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)
              aaaaaaaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGc                                                                                                                                         <  2:1186917/150‑1 (MQ=255)
               aaaaaaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCGAATTTAGCAAAAGTCtgatct                                                                                                                                        <  2:367214/150‑6 (MQ=255)
                aaaaaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGcac                                                                                                                                       >  1:210693/1‑150 (MQ=255)
                   aaCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACAtt                                                                                                                                    >  1:1289327/1‑150 (MQ=255)
                    aCCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACAttt                                                                                                                                   <  2:1129679/150‑1 (MQ=255)
                     ccGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACAtttt                                                                                                                                  <  1:1068488/150‑1 (MQ=255)
                     ccGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCGAATTTAGCAAAAGTCtgatc                                                                                                                                         <  1:655464/143‑5 (MQ=255)
                         aCCAAACCTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTa                                                                                                                              >  2:475264/1‑150 (MQ=255)
                             aaCGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAAttt                                                                                                                          >  1:979066/1‑150 (MQ=255)
                                      attattCCGGTGCGAGTATGACTTTTGAGGAATGGCTATCAGTATTGTCATTTTCAGACAATATTTATCCTGCATCGGTGAGGCAGAGTAAGCTGAGACTTTTGCTAAATTCGGAAATGACTTTGCACATTTTGCTAATTTCACCGTAcc                                                                                                                 <  2:506188/150‑1 (MQ=255)
                                      attattCCGGTGCGAGTAAGAATTATGATGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCTGTGAGTGAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTTCACATTTTGCTAATTTCACCGTAcc                                                                                                                 <  2:54175/150‑1 (MQ=255)
                                      attattCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTAcc                                                                                                                 <  1:1103098/150‑1 (MQ=255)
                                        tattCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGc                                                                                                               <  1:370664/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)
                                                        aGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACtata                                                                                               >  2:1284159/1‑150 (MQ=255)
                                                        aGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACtat                                                                                                >  1:555335/1‑149 (MQ=255)
                                                        aGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACtat                                                                                                <  2:555335/149‑1 (MQ=255)
                                                         taaTTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAg                                                                                              <  1:1591022/149‑1 (MQ=255)
                                                         gAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAg                                                                                              >  2:1111923/1‑150 (MQ=255)
                                                         gAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAg                                                                                              <  1:945394/150‑1 (MQ=255)
                                                               tGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGt                                                                                                                    <  1:19096/122‑1 (MQ=255)
                                                               tGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGt                                                                                                                    >  2:19096/1‑122 (MQ=255)
                                                               tGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCa                                                                                        <  2:1426017/150‑1 (MQ=255)
                                                               tGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCa                                                                                        >  1:491380/1‑150 (MQ=255)
                                                                gAGGAATGGCTATGAGTATTGTCATTTTCAGAAAATTTTTATCCTGCATCGTTGAGTCATAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCaa                                                                                       <  2:1049457/150‑1 (MQ=255)
                                                                   gAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACg                                                                                    <  2:939949/150‑1 (MQ=255)
                                                                   gAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACg                                                                                    <  2:1399257/150‑1 (MQ=255)
                                                                            tCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCtt                                                                           <  2:194734/150‑1 (MQ=255)
                                                                              tgTATTGTCATTTTCAGAAAATATTTATACTTCATCGTTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCTCAAAAGACTTTGCTCATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGAAAACGTCTCACCTTGa                                                                         <  2:221809/149‑1 (MQ=255)
                                                                              aGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGa                                                                         >  1:369377/1‑150 (MQ=255)
                                                                               gTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAg                                                                        >  2:797539/1‑150 (MQ=255)
                                                                                  ttGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGtt                                                                     >  1:1494691/1‑150 (MQ=255)
                                                                                    gTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCAcc                                                                             >  1:888724/1‑140 (MQ=255)
                                                                                    gTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTaa                                                                   >  2:1207862/1‑150 (MQ=255)
                                                                                          ttCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATg                                                             <  2:283250/150‑1 (MQ=255)
                                                                                           tCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGa                                                            <  2:432920/150‑1 (MQ=255)
                                                                                           tCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGa                                                            <  2:1603338/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)
                                                                                               aaaaTATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCaa                                                        <  1:1309585/150‑1 (MQ=255)
                                                                                                   tatTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATc                                                    >  1:168277/1‑150 (MQ=255)
                                                                                                    atTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGAAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATcc                                                   >  2:505054/1‑150 (MQ=255)
                                                                                                       tATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGa                                                >  2:216602/1‑150 (MQ=255)
                                                                                                        aTCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGAt                                               <  1:700943/150‑1 (MQ=255)
                                                                                                           cTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTAt                                            >  2:341343/1‑150 (MQ=255)
                                                                                                             gCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCa                                          >  2:1220618/1‑150 (MQ=255)
                                                                                                                  ggTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGgc                                     >  2:723669/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)
                                                                                                                             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)
                                                                                                                                    aagtcTTTTGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCAc                   >  2:1025983/5‑150 (MQ=255)
                                                                                                                                    aagtcTTTTGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGAACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGACCAc                   >  2:1455558/5‑150 (MQ=255)
                                                                                                                                     agtcTTTTGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCATGATTATCAACGGCGCGAAAAAATTCGACCACt                  >  1:627603/4‑150 (MQ=255)
                                                                                                                                      gtcTTTTGCGAATTTAGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTc                 >  2:719963/3‑150 (MQ=255)
                                                                                                                                        cTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCa               >  2:479177/1‑150 (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)
                                                                                                                                               tAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCa        <  1:1220618/150‑1 (MQ=255)
                                                                                                                                                 aaTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCANc      <  1:425687/150‑1 (MQ=255)
                                                                                                                                                 aaTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAAc      >  1:773137/1‑150 (MQ=255)
                                                                                                                                                 aaTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAAc      <  1:431643/150‑1 (MQ=255)
                                                                                                                                                  aTTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAATTCGCCCACTCCAATGGTCAACTg    <  1:20321/150‑1 (MQ=255)
                                                                                                                                                     cGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGaa  >  2:1510886/1‑150 (MQ=255)
                                                                                                                                                     cGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGaa  >  2:678381/1‑150 (MQ=255)
                                                                                                                                                     cGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGaa  >  1:1246285/1‑150 (MQ=255)
                                                                                                                                                     cGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGaa  >  1:1240506/1‑150 (MQ=255)
                                                                                                                                                     cGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGGCAACTGaa  <  2:262504/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TATTTTCAATCCCCAAAAAAACCGAACCAAACGTGTAAATTATTCCGGTGCGAGTAAGAATTATGAGGAATGGCTATCAGTATTGTCATTTTCAGAAAATATTTATCCTGCATCGGTGAGTCAGAGTAAGATCAGACTTTTGCTAAATTCGCAAAAGACTTTGCACATTTTGCTAATTTCACCGTACCGCTCTGTGACGTACTATAGTCGGCAAACGTCTCACCTTGAGGTTAAAAAATGAGCAACATCCTGATTATCAACGGCGCGAAAAAATTCGCCCACTCCAATGGTCAACTGAA  >  NC_000913/3172293‑3172591

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: