Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,679,430 Δ1 bp 8.3% intergenic (‑63/+34) yphF ← / ← yphG putative ABC transporter periplasmic binding protein YphF/DUF5107 domain‑containing protein YphG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,679,4300G. 8.3% 395.1 / 28.0 121intergenic (‑63/+34)yphF/yphGputative ABC transporter periplasmic binding protein YphF/DUF5107 domain‑containing protein YphG
Reads supporting (aligned to +/- strand):  ref base G (62/47);  new base . (5/5);  total (68/53)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

Only 100 of 124 total aligned reads displayed.

TTTTCCGCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCC  >  NC_000913/2679282‑2679579
                                                                                                                                                    |                                                                                                                                                      
    ccGCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGt                                                                                                                                                  >  2:143838/1‑150 (MQ=255)
     cGCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTg                                                                                                                                                 <  2:1286599/150‑1 (MQ=255)
      gCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGc                                                                                                                                                >  1:933828/1‑150 (MQ=255)
         gCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGc                                                                                                                                             <  2:1008313/150‑1 (MQ=255)
         gCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGc                                                                                                                                             >  1:574518/1‑150 (MQ=255)
         gCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGc                                                                                                                                             <  2:82061/150‑1 (MQ=255)
         gCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGc                                                                                                                                             <  2:933828/150‑1 (MQ=255)
            cTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCa                                                                                                                                          >  2:738641/1‑150 (MQ=255)
            cTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCa                                                                                                                                          >  2:603606/1‑150 (MQ=255)
            cTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCa                                                                                                                                          <  1:1198065/150‑1 (MQ=255)
             tGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCaa                                                                                                                                         <  1:603606/150‑1 (MQ=255)
             tGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCat                                                                                                                                         <  1:274759/150‑2 (MQ=255)
             tGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCat                                                                                                                                         <  1:1427395/150‑2 (MQ=255)
                cAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCatcac                                                                                                                                      <  1:628742/150‑5 (MQ=255)
                cAAATAGGGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCCCGCGCTGCCGCTCatcac                                                                                                                                      <  1:1284668/150‑5 (MQ=255)
                    tAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTAc                                                                                                                                  >  2:295740/1‑150 (MQ=255)
                     aGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGACTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACa                                                                                                                                 <  1:916566/150‑1 (MQ=255)
                      gTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAg                                                                                                                                >  1:1394653/1‑150 (MQ=255)
                      gTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCAGCTCAAATATTACAg                                                                                                                                >  1:1097229/1‑150 (MQ=255)
                          gcTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTaa                                                                                                                            >  2:1253159/1‑150 (MQ=255)
                           cTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAAt                                                                                                                           >  2:726912/1‑150 (MQ=255)
                              ccTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAgg                                                                                                                        >  1:1124274/1‑150 (MQ=255)
                               cTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAggg                                                                                                                       >  2:116815/1‑150 (MQ=255)
                                tAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGATCGGCAGCGCGTTCCGCTCAAATATTACAGTTAATAgtgt                                                                                                                      <  2:572331/150‑1 (MQ=255)
                                tAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGt                                                                                                                      >  1:394364/1‑150 (MQ=255)
                                 aaTAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGtt                                                                                                                     <  1:1634497/150‑1 (MQ=255)
                                 aaTAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGtt                                                                                                                     >  2:625076/1‑150 (MQ=255)
                                    aGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAg                                                                                                                  <  2:1097229/150‑1 (MQ=255)
                                     gCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGt                                                                                                                 >  1:1516443/1‑150 (MQ=255)
                                          gCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAAcg                                                                                                            >  2:243971/1‑150 (MQ=255)
                                          gCCATTAATAATAGATTACGGGTTGATCTCATTTTTGTAGGCATAGAGACTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAAATATTTTGTGATGAGCGGCAGCGCGTGAAGCTCAAATATTACAGTTAATAGGGTTTAGTTAAcg                                                                                                            >  2:553823/1‑150 (MQ=255)
                                           ccATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAat                                                                                                                                       >  1:1632716/1‑122 (MQ=255)
                                           ccATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAAcgc                                                                                                           >  2:383961/1‑150 (MQ=255)
                                            cATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAAcgcg                                                                                                          >  2:823802/1‑150 (MQ=255)
                                                  taataGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTaaa                                                                                                    >  2:1993/1‑150 (MQ=255)
                                                   aataGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTaaaa                                                                                                   >  1:308599/1‑150 (MQ=255)
                                                   aataGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTaaaa                                                                                                   >  1:1113886/1‑150 (MQ=255)
                                                          tACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCTCGTGCCTCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAg                                                                                            <  1:1630887/150‑1 (MQ=255)
                                                           aCGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAgg                                                                                           <  1:129561/150‑1 (MQ=255)
                                                             gggTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCa                                                                                         <  1:1009928/150‑1 (MQ=255)
                                                              ggTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCAt                                                                                        >  2:934934/1‑150 (MQ=255)
                                                                   ttCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACgg                                                                                   <  2:614025/150‑1 (MQ=255)
                                                                    tctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggt                                                                                  >  2:1017857/1‑150 (MQ=255)
                                                                    tctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggt                                                                                  >  1:984241/1‑150 (MQ=255)
                                                                     ctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtg                                                                                 >  2:915417/1‑150 (MQ=255)
                                                                      tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtgg                                                                                <  2:857901/150‑1 (MQ=255)
                                                                      tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtgg                                                                                <  2:1394653/150‑1 (MQ=255)
                                                                       cATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtggt                                                                               >  1:689168/1‑150 (MQ=255)
                                                                        aTTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGgtggtg                                                                              <  2:62384/150‑1 (MQ=255)
                                                                         tttttGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGa                                                                             <  2:30568/150‑1 (MQ=255)
                                                                         tattagTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGa                                                                             <  2:896922/145‑1 (MQ=255)
                                                                             tGTAGGCATACAGCCTCCTGTAGGCTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTCCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataa                                                                         >  2:1429210/1‑150 (MQ=255)
                                                                              gTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataat                                                                        >  2:394101/1‑150 (MQ=255)
                                                                                 ggCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGcc                                                                     >  2:326738/1‑150 (MQ=255)
                                                                                 ggCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGcc                                                                     >  1:379293/1‑150 (MQ=255)
                                                                                   cATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGc                                                                   >  2:1151051/1‑150 (MQ=255)
                                                                                   cATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGc                                                                   <  1:295740/150‑1 (MQ=255)
                                                                                      agagCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCttt                                                                >  2:905908/1‑150 (MQ=255)
                                                                                        agCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCg                                                              >  2:307802/1‑150 (MQ=255)
                                                                                         gCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCgg                                                             <  1:1993/150‑1 (MQ=255)
                                                                                          cctcctGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGc                                                            <  1:116815/150‑1 (MQ=255)
                                                                                           ctcctGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGcc                                                           >  1:1667820/1‑150 (MQ=255)
                                                                                               tGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGc                                                       >  1:1413553/1‑150 (MQ=255)
                                                                                                  aGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCgg                                                    >  1:1244780/1‑150 (MQ=255)
                                                                                                   gggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCa                                                                                                                                          >  2:210866/1‑63 (MQ=37)
                                                                                                    ggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAt                                                  <  1:647962/150‑1 (MQ=255)
                                                                                                     gTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAtt                                                 <  1:596948/150‑1 (MQ=255)
                                                                                                       ttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTTGTGAATAATGCCGCTTTCGGCCAAGCCGGATTaa                                               <  2:962101/150‑1 (MQ=255)
                                                                                                        tttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAAt                                              >  1:1658554/1‑150 (MQ=255)
                                                                                                         ttattaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATc                                             <  2:752491/150‑1 (MQ=255)
                                                                                                           attaACAAAGGCTTATTATAATTATTTTTTTATGATCGGCAGCGCGTGCCGCGCAAATATTAAAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAATGCATCACGGTGGTGAATAATGCCGCTTTAGGCCACGCCGGATTAATAGc                                           >  2:877170/1‑150 (MQ=255)
                                                                                                            ttaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGCTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGCCCAAGCCGGATTAATCGcc                                          >  1:409926/1‑150 (MQ=255)
                                                                                                             taACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCg                                         >  2:553971/1‑150 (MQ=255)
                                                                                                               acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTc                                       <  2:1534358/150‑1 (MQ=255)
                                                                                                               acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTc                                       <  2:394364/150‑1 (MQ=255)
                                                                                                                    ggCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTc                                  <  1:307802/150‑1 (MQ=255)
                                                                                                                     gCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATACTGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCa                                 >  2:796126/1‑150 (MQ=255)
                                                                                                                     gCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCa                                 >  1:1610024/1‑150 (MQ=255)
                                                                                                                      cTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAg                                <  2:1128441/150‑1 (MQ=255)
                                                                                                                      cTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAg                                <  1:1151051/150‑1 (MQ=255)
                                                                                                                          ttCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGcgcg                            >  2:1368582/1‑150 (MQ=255)
                                                                                                                            caaaTTATTTTGTGATGAGCGGCAGCGCCTGCCGATCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTTAATAGTGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAg                          >  1:426323/3‑150 (MQ=255)
                                                                                                                            cTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAg                          >  2:687501/1‑150 (MQ=255)
                                                                                                                            cTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGGTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAg                          <  1:93741/150‑1 (MQ=255)
                                                                                                                             tAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGc                         <  1:915417/150‑1 (MQ=255)
                                                                                                                             tAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGc                         <  2:1658554/150‑1 (MQ=255)
                                                                                                                              aaTTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCa                        <  1:967012/150‑1 (MQ=255)
                                                                                                                               attattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAg                       >  1:1044858/1‑150 (MQ=255)
                                                                                                                                 tattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGAc                     <  1:1160694/150‑1 (MQ=255)
                                                                                                                                       tatttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCatat               >  1:624646/5‑150 (MQ=255)
                                                                                                                                          ttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc            >  2:1079470/2‑150 (MQ=255)
                                                                                                                                          ttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc            >  1:200736/2‑150 (MQ=255)
                                                                                                                                          atgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGATTCATATTgc            <  1:257204/150‑1 (MQ=255)
                                                                                                                                           tgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgct           <  2:610016/150‑2 (MQ=255)
                                                                                                                                           tgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcg           >  1:597274/1‑150 (MQ=255)
                                                                                                                                             aGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcgcc         >  2:19220/1‑149 (MQ=255)
                                                                                                                                             aGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCa         <  2:1610024/150‑1 (MQ=255)
                                                                                                                                               cGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc       >  2:940073/1‑150 (MQ=255)
                                                                                                                                                  cAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc    >  1:1294596/1‑150 (MQ=255)
                                                                                                                                                    gcgcGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCGGCGCGAGCAGACTCATATTGCGCAACCTccc  <  2:574518/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTTTCCGCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCC  >  NC_000913/2679282‑2679579

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: