Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,679,434:1 +C 8.7% intergenic (‑67/+30) yphF ← / ← yphG putative sugar ABC transporter periplasmic binding protein/DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,679,4341.C 8.7% 561.9 / 53.0 138intergenic (‑67/+30)yphF/yphGputative sugar ABC transporter periplasmic binding protein/DUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base . (53/73);  new base C (5/7);  total (58/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.03e-01

Only 100 of 141 total aligned reads displayed.

GCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGA  >  NC_000913/2679288‑2679582
                                                                                                                                                   |                                                                                                                                                    
gCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTg                                                                                                                                                    >  3:455158/1‑149 (MQ=255)
 cAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGc                                                                                                                                                   <  2:92371/149‑1 (MQ=255)
   gCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCg                                                                                                                                                 >  4:475972/1‑149 (MQ=255)
     ccTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCt                                                                                                                                               <  2:870534/149‑1 (MQ=255)
      cTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTc                                                                                                                                              >  3:100395/1‑149 (MQ=255)
      cTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTc                                                                                                                                              <  1:1435432/149‑1 (MQ=255)
       tGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCa                                                                                                                                             <  4:549564/149‑1 (MQ=255)
         gCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTg                                                                                                                                                    >  3:14805/1‑140 (MQ=255)
         gCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTg                                                                                                                                                    <  4:14805/140‑1 (MQ=255)
           aaaTAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAata                                                                                                                                         <  3:661725/149‑1 (MQ=255)
              tAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTa                                                                                                                                      <  1:1035700/149‑1 (MQ=255)
              tAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTa                                                                                                                                      <  2:1039382/149‑1 (MQ=255)
                gTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACa                                                                                                                                    <  3:70480/149‑1 (MQ=255)
                gTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACa                                                                                                                                    <  4:1023252/149‑1 (MQ=255)
                gTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACa                                                                                                                                    <  1:976757/149‑1 (MQ=255)
                gTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACa                                                                                                                                    <  4:313482/149‑1 (MQ=255)
                gTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACa                                                                                                                                    <  1:1185101/149‑1 (MQ=255)
                  gcgcTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGt                                                                                                                                  <  2:1451400/149‑1 (MQ=255)
                      tACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCaaa                                                                                                                                           <  2:1639595/136‑1 (MQ=255)
                      tACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAAt                                                                                                                              >  3:1625664/1‑149 (MQ=255)
                       aCCTAATAGCGTAGCCATTAATAATAGGTTAGGGATTGGGCTAATTTTTGTAGGCATGGAGCCTCCTGTAGGGTGTTTATTAACAACGGCTTATTGTAATTATTTTGTGATGGGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATa                                                                                                                             >  2:342109/1‑149 (MQ=255)
                          tAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCTTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAggg                                                                                                                          >  1:1401121/1‑149 (MQ=255)
                            aTAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGtt                                                                                                                        <  3:545906/149‑1 (MQ=255)
                              aGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGNTTTa                                                                                                                      <  3:1324925/149‑1 (MQ=255)
                              aGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTa                                                                                                                      <  1:1489794/149‑1 (MQ=255)
                              aGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTa                                                                                                                      <  1:1287130/149‑1 (MQ=255)
                               gCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAg                                                                                                                     <  2:733170/149‑1 (MQ=255)
                                     ccATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAAcg                                                                                                               >  3:1476089/1‑149 (MQ=255)
                                       aTTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTc                                                                                                                                              >  2:1702438/1‑116 (MQ=255)
                                        ttAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCgg                                                                                                            <  1:185029/149‑1 (MQ=255)
                                               taGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTaaaaa                                                                                                     <  3:123948/149‑1 (MQ=255)
                                               taGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTaaaaa                                                                                                     <  3:1449915/149‑1 (MQ=255)
                                                 gATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCCGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGCTAATAGGGTTTAGTGAACGCGGTTAAAAATa                                                                                                   >  3:1016426/1‑149 (MQ=255)
                                                    tACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATaaaa                                                                                                <  3:475972/149‑1 (MQ=255)
                                                    tACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATaaaa                                                                                                <  3:1623900/149‑1 (MQ=255)
                                                    tACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATaaaa                                                                                                <  2:1450668/149‑1 (MQ=255)
                                                       gggTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGc                                                                                             >  4:593301/1‑149 (MQ=255)
                                                       gggTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGc                                                                                             <  2:915818/149‑1 (MQ=255)
                                                         gttgttCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCAt                                                                                           >  1:1683069/1‑149 (MQ=255)
                                                          ttgttCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATc                                                                                          <  2:1686034/149‑1 (MQ=255)
                                                             ttCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAatat                                                                                                                                        >  1:1715312/1‑100 (MQ=255)
                                                              tctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACgg                                                                                      >  4:625899/1‑149 (MQ=255)
                                                               ctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggt                                                                                     >  3:860638/1‑149 (MQ=255)
                                                               ctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggt                                                                                     >  1:466449/1‑149 (MQ=255)
                                                               ctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCa                                                                                                                                             <  2:1560980/93‑1 (MQ=255)
                                                               ctcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCACGCGCTGCCGCTCa                                                                                                                                             >  1:1560980/1‑93 (MQ=255)
                                                                tcgtttttGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtg                                                                                    >  4:1546831/4‑149 (MQ=255)
                                                                tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTTATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtg                                                                                    >  2:1707183/1‑149 (MQ=255)
                                                                tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtg                                                                                    <  1:1707183/149‑1 (MQ=255)
                                                                tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtg                                                                                    >  4:418918/1‑149 (MQ=255)
                                                                tcATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtg                                                                                    >  1:1670063/1‑149 (MQ=255)
                                                                  aTTTTTGTAGTCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTAAAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtggt                                                                                  <  4:1682500/149‑1 (MQ=255)
                                                                  aTTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGCTGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtggt                                                                                  >  2:296581/1‑149 (MQ=255)
                                                                  aTTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtggt                                                                                  <  2:1683069/149‑1 (MQ=255)
                                                                   tttttGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGgtggtg                                                                                 >  1:1289124/1‑149 (MQ=255)
                                                                    ttttGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGa                                                                                <  2:534928/149‑1 (MQ=255)
                                                                        gTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataa                                                                            >  2:366246/1‑149 (MQ=255)
                                                                         tAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataat                                                                           <  2:1125711/149‑1 (MQ=255)
                                                                         tAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataat                                                                           <  3:418918/149‑1 (MQ=255)
                                                                         tAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataat                                                                           <  2:212147/149‑1 (MQ=255)
                                                                            gCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAgg                                                                                              <  4:673369/127‑1 (MQ=255)
                                                                             cATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCg                                                                       <  2:1664788/149‑1 (MQ=255)
                                                                             cATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCg                                                                       <  4:567819/149‑1 (MQ=255)
                                                                                 gagCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCttt                                                                   >  2:1402820/1‑149 (MQ=255)
                                                                                   gCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCg                                                                 <  2:585838/149‑1 (MQ=255)
                                                                                    cctcctGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCgg                                                                <  4:12956/149‑1 (MQ=255)
                                                                                         tGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAg                                                                                                                     >  2:1517043/1‑91 (MQ=255)
                                                                                         tGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAg                                                                                                                     <  1:1517043/91‑1 (MQ=255)
                                                                                           tAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtggt                                                                                  >  4:138347/1‑124 (MQ=255)
                                                                                           tAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACggtggt                                                                                  <  3:138347/124‑1 (MQ=255)
                                                                                           tAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGcc                                                         >  4:1615078/1‑149 (MQ=255)
                                                                                             gggTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCgg                                                       >  1:922821/1‑149 (MQ=255)
                                                                                               gTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAt                                                     <  2:1472848/149‑1 (MQ=255)
                                                                                               gTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAt                                                     <  1:1575546/149‑1 (MQ=255)
                                                                                               gTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAt                                                     <  2:1289124/149‑1 (MQ=255)
                                                                                                tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCa                                                                                            <  2:1243579/109‑1 (MQ=255)
                                                                                                tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataat                                                                           >  2:469850/1‑126 (MQ=255)
                                                                                                tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataat                                                                           <  1:469850/126‑1 (MQ=255)
                                                                                                tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAtt                                                    >  3:778302/1‑149 (MQ=255)
                                                                                                tttttATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGAtt                                                    >  3:1676047/1‑149 (MQ=255)
                                                                                                      ttaACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACgg                                                                                      <  1:297325/109‑1 (MQ=255)
                                                                                                       taACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGcc                                             >  1:1015400/1‑149 (MQ=255)
                                                                                                         acaacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGt                                           >  1:902747/1‑149 (MQ=255)
                                                                                                           aacGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCa                                         >  1:179891/1‑149 (MQ=255)
                                                                                                            acGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCaa                                        <  1:1402820/149‑1 (MQ=255)
                                                                                                                cTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCa                                    <  4:1625664/149‑1 (MQ=255)
                                                                                                                cTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCa                                    <  2:466449/149‑1 (MQ=255)
                                                                                                                          ttattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAg                          >  4:267032/1‑149 (MQ=255)
                                                                                                                          ttattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAg                          >  3:896259/1‑149 (MQ=255)
                                                                                                                           tattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                         <  2:195137/149‑1 (MQ=255)
                                                                                                                           tattTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGa                         <  1:222671/149‑1 (MQ=255)
                                                                                                                                 gtgATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAGTTCAGCGCGAGCAGACTCata                   >  1:194228/1‑149 (MQ=255)
                                                                                                                                  tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTaaaa                                                                                                      <  3:630399/65‑1 (MQ=255)
                                                                                                                                  atttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGaataat                                                                           <  3:310108/89‑1 (MQ=38)
                                                                                                                                      gaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcg              >  4:1450030/1‑149 (MQ=255)
                                                                                                                                       aGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCgg                                                       >  3:767483/1‑107 (MQ=255)
                                                                                                                                       aGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCgg                                                       <  4:767483/107‑1 (MQ=255)
                                                                                                                                             aGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc       <  3:1225174/149‑1 (MQ=255)
                                                                                                                                               cgcgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTccc     <  1:1498845/149‑1 (MQ=255)
                                                                                                                                                  gTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGa  <  1:555630/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCAGCCCTGGCAAATAGTGCGCTACCTAATAGCGTAGCCATTAATAATAGATTACGGGTTGTTCTCATTTTTGTAGGCATAGAGCCTCCTGTAGGGTTTTTATTAACAACGGCTTATTCTAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGA  >  NC_000913/2679288‑2679582

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: