Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,779 A→T 5.5% intergenic (‑19/+49) ycjD ← / ← sapF DUF559 domain‑containing protein YcjD/putrescine ABC exporter ATP binding protein SapF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,351,7790AT 5.5% 524.1 / 12.3 163intergenic (‑19/+49)ycjD/sapFDUF559 domain‑containing protein YcjD/putrescine ABC exporter ATP binding protein SapF
Reads supporting (aligned to +/- strand):  ref base A (79/75);  new base T (7/2);  total (86/77)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

Only 100 of 164 total aligned reads displayed.

CCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCAC  >  NC_000913/1351631‑1351923
                                                                                                                                                    |                                                                                                                                                
ccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTaa                                                                                                                                                 >  1:566458/1‑150 (MQ=255)
 cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTTATTCTGACAAGATTTAAt                                                                                                                                                <  1:1862016/150‑1 (MQ=255)
  aCTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATc                                                                                                                                               >  1:1231096/1‑150 (MQ=255)
   cTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCt                                                                                                                                              <  1:41422/150‑1 (MQ=255)
       aTGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                                                                                                                          >  2:1283271/1‑150 (MQ=255)
        tGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGt                                                                                                                                         <  1:2355518/150‑1 (MQ=255)
         gTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTc                                                                                                                                        >  2:221877/1‑150 (MQ=255)
          ttGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCa                                                                                                                                       <  2:815279/150‑1 (MQ=255)
           tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAg                                                                                                                                      >  1:96811/1‑150 (MQ=255)
           tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAg                                                                                                                                      >  1:2407688/1‑150 (MQ=255)
            gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGa                                                                                                                                     >  2:736861/1‑150 (MQ=255)
            gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGa                                                                                                                                     >  2:1581030/1‑150 (MQ=255)
            gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGa                                                                                                                                     >  2:2001006/1‑150 (MQ=255)
                gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGa                                                                                                                                 >  1:921885/1‑150 (MQ=255)
                 cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                <  1:221877/150‑1 (MQ=255)
                  gAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAgg                                                                                                                               <  1:389252/150‑1 (MQ=38)
                         gAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGaataa                                                                                                                        >  1:273902/1‑150 (MQ=255)
                          aaaTCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGaataat                                                                                                                       >  2:1993350/1‑150 (MQ=255)
                             tCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTg                                                                                                                    >  2:304828/1‑150 (MQ=255)
                             tCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTg                                                                                                                    <  1:1411564/150‑1 (MQ=255)
                             tCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTg                                                                                                                    <  1:1177706/150‑1 (MQ=255)
                                     ccctCGGCTGCGACGATATCGCCAGAGCTGTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAAATTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggagg                                                                                                            <  2:1988416/147‑1 (MQ=255)
                                     aCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggagg                                                                                                            <  1:2236030/150‑1 (MQ=255)
                                      cGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAggagga                                                                                                           >  1:1272063/1‑150 (MQ=255)
                                          ggCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCa                                                                                                       <  2:60150/150‑1 (MQ=255)
                                           gCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAg                                                                                                      >  1:964617/1‑150 (MQ=255)
                                              gCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaa                                                                                                   <  2:1716967/150‑1 (MQ=255)
                                              gCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCAATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaa                                                                                                   <  2:1157686/150‑1 (MQ=255)
                                               cGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaaa                                                                                                  <  1:1205187/150‑1 (MQ=255)
                                                   gATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAg                                                                                              <  1:2051513/150‑1 (MQ=255)
                                                   gATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAg                                                                                              >  2:492879/1‑150 (MQ=255)
                                                    atatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGc                                                                                             <  2:669140/150‑1 (MQ=255)
                                                     tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                            >  2:2211318/1‑150 (MQ=255)
                                                     tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                            >  2:716772/1‑150 (MQ=255)
                                                     tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                            <  1:524514/150‑1 (MQ=255)
                                                        cGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGAta                                                                                         >  1:2365597/1‑149 (MQ=255)
                                                         gCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCt                                                                                        <  2:2125884/150‑1 (MQ=255)
                                                         gCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCt                                                                                        >  2:1123445/1‑150 (MQ=255)
                                                            agagCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTa                                                                                     >  1:1484826/1‑150 (MQ=255)
                                                            agagCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCAATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTa                                                                                     <  1:1646673/150‑1 (MQ=255)
                                                              agCTTTCGTTCCTGCAAACTGAGATTGCGTCTTAAATCACGCGCATTTCATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACg                                                                                   <  2:911857/150‑1 (MQ=255)
                                                               gCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGc                                                                                  >  2:572622/1‑150 (MQ=255)
                                                                 tttCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCa                                                                                >  1:2073419/1‑150 (MQ=255)
                                                                    cGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGc                                                                             >  1:638659/1‑150 (MQ=255)
                                                                    cGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGc                                                                             >  2:1012287/1‑150 (MQ=255)
                                                                         cTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccg                                                                        <  2:196413/150‑1 (MQ=255)
                                                                         cTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccg                                                                        >  1:1122891/1‑150 (MQ=255)
                                                                         cTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccg                                                                        >  2:1027642/1‑150 (MQ=255)
                                                                          tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAGTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccgc                                                                       >  2:2252922/1‑150 (MQ=255)
                                                                          tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccgc                                                                       >  2:541318/1‑150 (MQ=255)
                                                                          tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccgc                                                                       >  2:552324/1‑150 (MQ=255)
                                                                             aaaGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTATGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGt                                                                    <  2:2392598/150‑1 (MQ=255)
                                                                              aaGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTc                                                                   <  1:437432/150‑1 (MQ=255)
                                                                                 tGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAAt                                                                <  2:1786184/150‑1 (MQ=255)
                                                                                    gaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGcc                                                             >  2:18842/1‑150 (MQ=255)
                                                                                      ttGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTc                                                           >  2:120711/1‑150 (MQ=255)
                                                                                      ttGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTc                                                           >  2:948202/1‑150 (MQ=255)
                                                                                        gCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGc                                                         >  1:769137/1‑150 (MQ=255)
                                                                                        gCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGc                                                         <  1:230413/150‑1 (MQ=255)
                                                                                         cGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                        <  1:1016322/150‑1 (MQ=255)
                                                                                           tCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCaa                                                      >  1:2100016/1‑150 (MQ=255)
                                                                                           tCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCaa                                                      >  1:15198/1‑150 (MQ=255)
                                                                                              tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                   >  2:596566/1‑150 (MQ=255)
                                                                                                  tCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGAcc                                               <  2:134897/150‑1 (MQ=255)
                                                                                                     cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                            <  1:652910/150‑1 (MQ=255)
                                                                                                      gcgcATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCa                                           <  1:2482328/150‑1 (MQ=255)
                                                                                                      gcgcATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCa                                           <  1:1575754/150‑1 (MQ=255)
                                                                                                        gcATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAAt                                         <  2:96811/150‑1 (MQ=255)
                                                                                                         cATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATc                                        >  1:637709/1‑150 (MQ=255)
                                                                                                           tttGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAg                                      <  2:2455971/150‑1 (MQ=255)
                                                                                                           tttGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAg                                      <  1:1660656/150‑1 (MQ=255)
                                                                                                             tGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGcc                                    >  2:2212444/1‑150 (MQ=255)
                                                                                                               aTTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGt                                  <  2:2453848/150‑1 (MQ=255)
                                                                                                                tttAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGtt                                 >  2:230673/1‑150 (MQ=255)
                                                                                                                  tAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTg                               >  1:541610/1‑150 (MQ=255)
                                                                                                                    aTTTTATCTATCACCTCATTCTGACAATTTTTAATCTTATGTCAGAATGAGGTGAATAATCTGGAAGTAGGATTCAGAAAATTAGCGATCTTTACGCCAGGCAGCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCATCCGTTTGGt                             <  2:2164392/150‑1 (MQ=255)
                                                                                                                    aTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGt                             >  2:828869/1‑150 (MQ=255)
                                                                                                                     ttttATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                            <  2:1357518/150‑1 (MQ=255)
                                                                                                                         atctatcACCTAATTCTGACAAGATTTAATCTTTTGTAAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGtt                        >  2:688543/1‑150 (MQ=255)
                                                                                                                            tatcACCTCTTTCTGACAAGATTTAATCTTTTGTCAGAATGAGTTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCAt                     <  1:20698/150‑1 (MQ=255)
                                                                                                                            tatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCAt                     <  2:1516093/150‑1 (MQ=255)
                                                                                                                            tatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCAt                     <  1:1079595/150‑1 (MQ=255)
                                                                                                                              tcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                   <  1:1045492/150‑1 (MQ=255)
                                                                                                                              tcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                   <  2:712164/150‑1 (MQ=255)
                                                                                                                                aCCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAg                 >  1:1661610/1‑150 (MQ=38)
                                                                                                                                 ccTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTAAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgc                >  1:389765/1‑150 (MQ=255)
                                                                                                                                   tCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcgg              >  1:1546102/1‑150 (MQ=255)
                                                                                                                                   tCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcgg              >  1:1326965/1‑150 (MQ=38)
                                                                                                                                   tCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcgg              >  1:1970086/1‑150 (MQ=38)
                                                                                                                                   tCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcgg              >  1:1752576/1‑150 (MQ=38)
                                                                                                                                   tCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGATAATAATCTGGAAGTAGAATTAAGAAAATTATAGATCTTTACAACACGCATCAGCCGTAAATAAATAACCAAAATCACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcgg              >  2:1551729/1‑150 (MQ=16)
                                                                                                                                    cATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcggc             >  1:380231/1‑150 (MQ=38)
                                                                                                                                     aTTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcggcg            <  1:835990/150‑1 (MQ=255)
                                                                                                                                        cTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATg         <  1:2212444/150‑1 (MQ=255)
                                                                                                                                        cTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATg         >  1:4807/1‑150 (MQ=255)
                                                                                                                                          gACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGcc       >  2:26836/1‑150 (MQ=255)
                                                                                                                                          gACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGcc       <  2:1703622/150‑1 (MQ=255)
                                                                                                                                           aCAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCa      <  1:1331491/150‑1 (MQ=255)
                                                                                                                                              aGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGca   <  1:2211318/150‑1 (MQ=255)
                                                                                                                                              aGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGca   <  2:1069623/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
CCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCAC  >  NC_000913/1351631‑1351923

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