Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,580,806 T→C 5.4% intergenic (‑16/+36) ydeM ← / ← ydeN putative YdeN‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,580,8060TC 5.4% 223.9 / 3.1 74intergenic (‑16/+36)ydeM/ydeNputative YdeN‑specific sulfatase‑maturating enzyme/putative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (34/36);  new base C (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCATCTCTTT  >  NC_000913/1580662‑1580955
                                                                                                                                                |                                                                                                                                                     
gatgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAACAAggg                                                                                                                                                  >  2:93958/1‑149 (MQ=255)
gatgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAACAAggg                                                                                                                                                  >  2:14054/1‑149 (MQ=255)
 atgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAACAAgggg                                                                                                                                                 <  1:14054/150‑3 (MQ=255)
 atgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAACAAgggg                                                                                                                                                 <  1:93958/150‑3 (MQ=255)
  tgaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGAAATAACATTTGAGATTACATTGAAAACTGGAGGCCTTGTCTGTAACGTGCATATCTAACTCGCTCAATAAggcgg                                                                                                                                                >  2:716458/1‑150 (MQ=255)
   gaACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAggcggc                                                                                                                                               >  2:633973/1‑150 (MQ=255)
   acaCTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAggcggc                                                                                                                                               >  2:344387/3‑150 (MQ=255)
    aaCTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGgcggcg                                                                                                                                              >  1:758373/1‑150 (MQ=255)
        ctTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGaaa                                                                                                                                          >  2:246255/1‑150 (MQ=255)
        ctTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGaaa                                                                                                                                          >  2:658416/1‑150 (MQ=255)
        ctTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGaaa                                                                                                                                          >  2:796318/1‑150 (MQ=255)
             agagTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATcc                                                                                                                                     >  1:1337127/1‑150 (MQ=255)
               agTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTATATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATccgc                                                                                                                                   >  1:317305/1‑150 (MQ=255)
               agTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATccgc                                                                                                                                   >  1:1276428/1‑150 (MQ=255)
                 tGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATccgccg                                                                                                                                 >  2:274936/1‑150 (MQ=255)
                   cTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCa                                                                                                                               <  2:657692/150‑1 (MQ=255)
                   cTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCa                                                                                                                               >  2:529661/1‑150 (MQ=255)
                         accatccatTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAgg                                                                                                                         <  2:324861/149‑1 (MQ=255)
                            atccatTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGttt                                                                                                                      >  2:340347/1‑150 (MQ=255)
                            atccatTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTCAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGttt                                                                                                                      <  1:902499/150‑1 (MQ=255)
                             tccatTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTa                                                                                                                     >  1:775597/1‑150 (MQ=255)
                               catTTTTCATGTGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTTCATTGAAAACTGGATGGCTTGGCTGTAACGTGCAACGCTATCTAGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGt                                                                                                                   <  2:426066/150‑1 (MQ=255)
                                atTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCTGCGGAAAAATCCGCCGCATGAAGGTTTAGtt                                                                                                                  >  1:1161465/1‑150 (MQ=255)
                                     tCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATttc                                                                                                             >  2:633280/1‑150 (MQ=255)
                                      cATGAGTAAACTGCGCCTCTTTTTCAAGCTAAAAACAGTAATCACATTTGAGATGACATTGAAAACTGGACGGCTTGGCTGTAACTTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCAGGTAGGTTTAGTTATttcc                                                                                                            >  2:809555/1‑149 (MQ=255)
                                      cATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATttcg                                                                                                            >  1:1206934/1‑150 (MQ=255)
                                         gAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATttcgctt                                                                                                         >  2:451851/1‑150 (MQ=255)
                                         gAGTAAACTGCGACTCCTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATttcgctt                                                                                                         >  1:494526/1‑150 (MQ=255)
                                             ataCTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGTAGGGCTTGGCTGTAACGTGCATCGATATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTtcgcttcgct                                                                                                     <  1:281321/148‑1 (MQ=255)
                                              aaCTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTcgcttcgctt                                                                                                    >  1:801087/1‑150 (MQ=255)
                                                   cGACTCATTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGATAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGc                                                                                               <  1:242318/150‑1 (MQ=255)
                                                    gACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCt                                                                                              >  2:544372/1‑150 (MQ=255)
                                                      ctctTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCttt                                                                                            >  2:55044/1‑150 (MQ=255)
                                                       tctTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCTCCGAATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTc                                                                                           >  1:383993/1‑150 (MQ=255)
                                                          ttttCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTg                                                                                        <  2:537181/150‑1 (MQ=255)
                                                            ttCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGat                                                                                      >  2:132767/1‑150 (MQ=255)
                                                              cAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGatat                                                                                    <  1:339989/150‑1 (MQ=255)
                                                              cAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGatat                                                                                    <  1:302885/150‑1 (MQ=255)
                                                               aaGGTAAAAACAGTAATGACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCTCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATAtt                                                                                   <  1:555553/150‑1 (MQ=255)
                                                                   tAAAAACAGTAATCACTTTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTGAATAATGCGGCGGAAAAATCCGCCGCATGAATGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTa                                                                               <  2:697089/150‑1 (MQ=255)
                                                                   tAAAAACAGTAATCACATTTGATATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTa                                                                               <  1:80488/150‑1 (MQ=255)
                                                                      aaaCAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTCCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGATGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAAc                                                                            >  2:210436/1‑150 (MQ=255)
                                                                       aaCAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACt                                                                           <  1:274936/150‑1 (MQ=255)
                                                                         taGTAATCACATTTGATATTATATAGTAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTGTCTCGCTCAATAAGGCGGCGGAAAAATCCGACGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTtc                                                                         <  1:1307746/149‑1 (MQ=255)
                                                                            tAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTcc                                                                      <  2:810790/150‑1 (MQ=255)
                                                                              aTCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAATCTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGAAATTGTTAAACTTATCCTg                                                                    >  2:1105796/1‑150 (MQ=255)
                                                                                 acaTTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGAtt                                                                 <  2:625763/150‑1 (MQ=255)
                                                                                           ttACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTa                                                       <  1:292574/150‑1 (MQ=255)
                                                                                           ttACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTa                                                       <  2:1265369/150‑1 (MQ=255)
                                                                                             aCATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAg                                                     <  1:340347/150‑1 (MQ=255)
                                                                                                 tGAAAACTGGAGGGCTTGGCTGTAACGTGCATTGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGt                                                 >  2:946924/1‑150 (MQ=255)
                                                                                                  gAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTGAGTGCTTTCTTGATATTGTTTAACTTCTCCTGATTTACCTCTCTAAGCGGtg                                                >  1:974234/1‑150 (MQ=255)
                                                                                                      aCTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCAGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGtggctg                                            >  2:274366/1‑150 (MQ=255)
                                                                                                       cTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGtggctgg                                           <  2:331761/150‑1 (MQ=255)
                                                                                                        tGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGtggctggc                                          <  1:769078/150‑1 (MQ=255)
                                                                                                        tGGAGGGCGTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGtggctggc                                          <  1:659836/150‑1 (MQ=255)
                                                                                                           cgggCTTGGCTGTTACGGGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGGTTAGTTTTTTCCCTTCGCTTAGAGCATGCTTGATATTGATAAACTTCTCCTGAATTACCTCGCTAAGCGGTGGCTGgctgc                                       <  2:68160/149‑1 (MQ=255)
                                                                                                               cTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTc                                   <  1:633973/150‑1 (MQ=255)
                                                                                                                ttGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCg                                  <  1:796318/150‑1 (MQ=255)
                                                                                                                ttGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCg                                  <  1:344387/150‑1 (MQ=255)
                                                                                                                ttGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCg                                  <  1:658416/150‑1 (MQ=255)
                                                                                                                ttGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCg                                  <  1:246255/150‑1 (MQ=255)
                                                                                                                       tAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAAct                           <  2:1097815/150‑1 (MQ=255)
                                                                                                                             gCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATTAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGAGGTGGCTGGCTGCTGTCGATAAACTCTCTTa                     <  2:1194728/150‑1 (MQ=255)
                                                                                                                              cATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTTTTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTAc                    <  1:490775/150‑1 (MQ=255)
                                                                                                                               aTCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTAcc                   <  1:1123475/150‑1 (MQ=255)
                                                                                                                               aTCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCAGGAAGGTTTAGTTATTTCGCTTCGCTTAGTGATTTCTTGATATTGTTAAACTTCTCCAGATTTACCTCGCTAAGCGGTGGCAGTCTGCTGTCGATAAACTCTCTTaac                   >  1:258593/1‑148 (MQ=255)
                                                                                                                                  gCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACg                <  2:445178/150‑1 (MQ=255)
                                                                                                                                     aTCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCt             >  1:1018990/1‑150 (MQ=255)
                                                                                                                                     aTCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCt             <  1:655894/150‑1 (MQ=255)
                                                                                                                                        tcGCTCAATAAGTCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGc          <  1:633280/150‑1 (MQ=255)
                                                                                                                                        tcGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGc          >  2:267714/1‑150 (MQ=255)
                                                                                                                                         cGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCa         <  2:318446/150‑1 (MQ=255)
                                                                                                                                              aataatGCGGCGGAAAAATCCGCCGGATGAATGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGTTATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCATctct    <  2:614565/150‑1 (MQ=255)
                                                                                                                                               aTAAGGCGGCGGAAAAATCCGCCGCAAGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCATCTCtt   >  2:498328/1‑150 (MQ=255)
                                                                                                                                                tAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCATCTCttt  <  2:758373/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GATGAACTCTTTCAGAGTGCTGTCATCCATCCATTTTTCATGAGTAAACTGCGACTCTTTTTCAAGGTAAAAACAGTAATCACATTTGAGATTACATTGAAAACTGGAGGGCTTGGCTGTAACGTGCATCGCTATCTCGCTCAATAAGGCGGCGGAAAAATCCGCCGCATGAAGGTTTAGTTATTTCGCTTCGCTTAGTGCTTTCTTGATATTGTTAAACTTCTCCTGATTTACCTCGCTAAGCGGTGGCTGGCTGCTGTCGATAAACTCTCTTACCACGCCTTGCATCTCTTT  >  NC_000913/1580662‑1580955

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