Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,585,907 C→T 5.6% P285P (CCG→CCA hsdR ← endonuclease R Type I restriction enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,585,9070CT 5.6% 227.7 / 2.4 71P285P (CCG→CCAhsdRendonuclease R Type I restriction enzyme
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/2);  total (34/37)
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.73e-01

AAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTG  >  NC_000913/4585761‑4586052
                                                                                                                                                  |                                                                                                                                                 
aaaCATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAACCAGCGCATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCgg                                                                                                                                                 >  4:158490/1‑149 (MQ=255)
aaaCATTTACAGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGCTATTGTTACGTTTCGCCTCTACAACCGCGATGGGGTTGAGGCCGACGAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCGGCTTTGCCGGTctg                                                                                                                                                 >  6:18800/1‑147 (MQ=255)
   cATTTACTGTAGCGATACGACTCATTGAGCCTGGAGGGAACTTCGATATTGTTACGTTTCGCATCTACCAGCGCGATGGGTTTGATGCCGACAAACAGCACATAATCCGAAAAGCCCTTATTACCCGTTTCATCTTTTCCGGTCGGCCa                                                                                                                                              <  8:224740/149‑1 (MQ=255)
     tttACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCAtt                                                                                                                                            <  6:226799/149‑1 (MQ=255)
     tttACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCAtt                                                                                                                                            >  4:155356/1‑149 (MQ=255)
                gATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCt                                                                                                                                 <  3:155356/149‑1 (MQ=255)
                 aTACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCtt                                                                                                                                <  2:176071/149‑1 (MQ=255)
                    cGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTa                                                                                                                              <  6:193332/148‑1 (MQ=255)
                     gACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTAtt                                                                                                                            >  8:232902/1‑149 (MQ=255)
                     gACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGAGGGGTTTGAGGCCGACAAACATCACATAATCCGCAAAGCCCTGCTTACCCGTTTCATCTTTTCCGGTCGTCCATTCGGCAATGGCTTTAtt                                                                                                                            >  4:155780/1‑149 (MQ=255)
                        tCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCt                                                                                                                                 >  3:150321/1‑141 (MQ=255)
                        tCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCt                                                                                                                                 <  4:150321/141‑1 (MQ=255)
                          aTTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGc                                                                                                                       >  8:130965/1‑149 (MQ=255)
                             gAGCCTGGCGGGAACGGCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGCGGCCGACAAACAGCACATAATCCGCACAGCCCTGATTCCCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCg                                                                                                                     <  3:155780/148‑1 (MQ=255)
                                 cTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTc                                                                                                                >  1:212670/1‑149 (MQ=255)
                                  tGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCAGACAAACAGCACATAATCCGCCAAGCCCTGATTAGCCGTTTCATCTTTTCCTGTCGGCCATTCGGCAATGGCTTTATTGACTCCGGGTTcc                                                                                                               >  6:224598/1‑149 (MQ=255)
                                      gggAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGAc                                                                                                           <  7:232902/149‑1 (MQ=255)
                                            gTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcgc                                                                                                     >  8:67629/1‑149 (MQ=255)
                                            gTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcgc                                                                                                     <  7:151362/149‑1 (MQ=255)
                                            gTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTgcgc                                                                                                     <  2:212670/149‑1 (MQ=255)
                                                atatTGTTACGTTTCGCCTCTACAACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCttt                                                                                                 >  8:66181/1‑149 (MQ=255)
                                                   ttGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGGCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTgg                                                                                               <  8:48752/148‑1 (MQ=255)
                                                   ttGTTACGATTCGCCTCTTCCACCGCGATGGGGTTGATGCCGACAATCAGCACATAATCCGCAAAGCCCTGAATACCCGTTTCATCTTTACCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGTTTTCCGGACGTGCGCCTTTGga                                                                                              <  7:161493/149‑1 (MQ=255)
                                                     gTTACGCTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATACTCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGgag                                                                                             <  6:206467/148‑1 (MQ=255)
                                                      ttACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTg                                                                                                >  8:50046/1‑144 (MQ=255)
                                                      ttACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTg                                                                                                <  7:50046/144‑1 (MQ=255)
                                                      ttACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGgag                                                                                             >  6:228467/1‑147 (MQ=255)
                                                      ttACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGgag                                                                                             >  2:157831/1‑147 (MQ=255)
                                                        aCGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAgc                                                                                         >  6:193240/1‑149 (MQ=255)
                                                         cGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAgcg                                                                                        <  1:29714/149‑1 (MQ=255)
                                                           tttCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCa                                                                                      <  3:89433/149‑1 (MQ=255)
                                                             tCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAg                                                                                     <  8:225702/148‑1 (MQ=255)
                                                                ccTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGtt                                                                                 >  2:89026/1‑149 (MQ=255)
                                                                   ctACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAACGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTg                                                                              <  1:240672/149‑1 (MQ=255)
                                                                    tACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTTGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGc                                                                             >  4:188324/1‑149 (MQ=255)
                                                                    tACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTTGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGc                                                                             >  2:142900/1‑149 (MQ=255)
                                                                    tACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGtttt                                                                               <  7:130965/147‑1 (MQ=255)
                                                                    tACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCAATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGtttt                                                                               <  2:149814/147‑1 (MQ=255)
                                                                      ccaccGCGATGGGTTTGAGGCCGACAAACAGCCCATAATCCGCAAAGCCCGGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTg                                                                           >  6:82955/1‑149 (MQ=255)
                                                                      ccaccGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGaa                                                                                           >  7:209577/1‑133 (MQ=255)
                                                                      ccaccGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGaa                                                                                           <  8:209577/133‑1 (MQ=255)
                                                                      ccaccGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTg                                                                              <  6:30424/146‑1 (MQ=255)
                                                                      ccaccGCGATGGGTTTGAGGACGCCAAACAGCACATAATCCGCAAAGCCCTGAATACCCGTTTCAGCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGGGCGCCTGTGGAGACGCGCAGGGTTTTGCTg                                                                           <  5:18800/149‑1 (MQ=255)
                                                                                   tttGAGGCCGACAAACAGCACATACTGCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCTGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGc                                                              <  6:166068/149‑1 (MQ=255)
                                                                                   tttGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTGATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGc                                                              >  7:167198/1‑149 (MQ=255)
                                                                                                aaCAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAg                                                 >  8:33027/1‑149 (MQ=255)
                                                                                                 aCAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGTCGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGt                                                <  3:45834/149‑1 (MQ=255)
                                                                                                  cAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGtt                                               >  3:75327/1‑149 (MQ=255)
                                                                                                     cacaTAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGCTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTgcg                                            >  5:63984/1‑149 (MQ=255)
                                                                                                     cacaTAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCAATCGGCAATCGCTTTATTGACGCCTGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTgcg                                            >  1:146401/1‑149 (MQ=255)
                                                                                                       caTAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTg                                                                                                >  3:32911/1‑95 (MQ=255)
                                                                                                       caTAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTg                                                                                                <  4:32911/95‑1 (MQ=255)
                                                                                                         tAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCatc                                        <  7:66181/149‑1 (MQ=255)
                                                                                                         tAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCat                                         >  8:161443/1‑148 (MQ=255)
                                                                                                           aTCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGTTTCCGGACGTGCGCCTTTGGAGAAGCTCAGGGTTTTCCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCatcaa                                      >  5:99141/1‑149 (MQ=255)
                                                                                                            tCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCatcaat                                     <  8:167198/149‑1 (MQ=255)
                                                                                                               gCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGtt                                               >  6:194030/1‑136 (MQ=255)
                                                                                                               gCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGtt                                               <  5:194030/136‑1 (MQ=255)
                                                                                                               gCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAg                                  >  6:139665/1‑149 (MQ=255)
                                                                                                               gCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAg                                  >  3:91316/1‑149 (MQ=255)
                                                                                                                    gCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATGCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGc                             <  1:230534/149‑1 (MQ=255)
                                                                                                                      ccTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGa                           <  7:161443/149‑1 (MQ=255)
                                                                                                                       cTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCCGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGAc                          >  5:69375/1‑149 (MQ=255)
                                                                                                                              cccGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCtcttc                   <  5:82955/149‑1 (MQ=255)
                                                                                                                                  tttCATCTTTTCCGGTTGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCg               <  1:142900/149‑1 (MQ=255)
                                                                                                                                  tttCATCTTTTCCGGTTGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCg               <  3:188324/149‑1 (MQ=255)
                                                                                                                                       tCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAg          >  6:138818/1‑149 (MQ=255)
                                                                                                                                           ttCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAgg         <  5:193240/146‑1 (MQ=255)
                                                                                                                                           ttCCGGGCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGtt       <  5:228467/148‑1 (MQ=255)
                                                                                                                                               ggTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTg                                                                              <  4:159110/73‑1 (MQ=255)
                                                                                                                                               ggTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAgtg  >  6:153577/1‑149 (MQ=255)
                                                                                                                                               ggTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGGTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTg                                                                              >  3:159110/1‑73 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AAACATTTACTGTAGCGATACGACTCATTGAGCCTGGCGGGAACGTCGATATTGTTACGTTTCGCCTCTACCACCGCGATGGGTTTGAGGCCGACAAACAGCACATAATCCGCAAAGCCCTGATTACCCGTTTCATCTTTTCCGGTCGGCCATTCGGCAATGGCTTTATTGACGCCGGGTTCCGGACGTGCGCCTTTGGAGAAGCGCAGGGTTTTGCTGTCGGCCTGCCAGCCTGCTTTACGCAGTTGCGCATCAATCAGGAAGCGACTCTCTTCTTCGCTAAGGTTGAGTG  >  NC_000913/4585761‑4586052

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