Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,048,429 C→T 7.3% W88* (TGG→TAG) ‡ gfcC ← putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,4290CT 7.3% 218.7 / 6.3 68W88* (TGG→TAG) ‡gfcCputative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base C (33/30);  new base T (3/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

CCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAAC  >  NC_000913/1048281‑1048577
                                                                                                                                                    |                                                                                                                                                    
cccGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTccc                                                                                                                                                      >  1:625175/1‑149 (MQ=255)
   gACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAgg                                                                                                                                                   >  3:132209/1‑149 (MQ=255)
     ccAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaa                                                                                                                                                 <  1:41003/149‑1 (MQ=255)
     ccAACGGGGGATTGCTGTTATCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGaa                                                                                                                                                 <  4:422168/149‑1 (MQ=255)
            gggATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCt                                                                                                                                          >  3:160652/1‑149 (MQ=255)
            gggATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCt                                                                                                                                          <  2:410890/149‑1 (MQ=255)
                       tctcATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACa                                                                                                                               <  1:231788/149‑1 (MQ=255)
                        ctcATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACAt                                                                                                                              >  3:85571/1‑149 (MQ=255)
                         tcATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATg                                                                                                                             >  4:123753/1‑149 (MQ=255)
                           aTCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGtt                                                                                                                           >  4:102322/1‑149 (MQ=255)
                            tCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGa                                                                                                                                 >  1:136823/1‑142 (MQ=255)
                            tCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGa                                                                                                                                 <  2:136823/142‑1 (MQ=255)
                                cTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagcagaa                                                                                                                                           >  4:466835/1‑125 (MQ=255)
                                cTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCTGGGaagcagaa                                                                                                                                           <  3:466835/128‑4 (MQ=255)
                                 tCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTg                                                                                                                          >  3:594285/1‑144 (MQ=255)
                                 tCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTATGCCATGACATGTTg                                                                                                                          <  4:594285/144‑1 (MQ=255)
                                     aCAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCttt                                                                                                                 >  1:145645/1‑149 (MQ=255)
                                       aaaaTCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTaa                                                                                                               <  3:586709/149‑1 (MQ=255)
                                       aaaaTCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTaa                                                                                                               <  1:380164/149‑1 (MQ=255)
                                           tCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAgcgc                                                                                                           <  1:583398/149‑1 (MQ=255)
                                             aGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCtt                                                                                                         <  3:65301/149‑1 (MQ=255)
                                               ggTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTa                                                                                                       <  3:248914/149‑1 (MQ=255)
                                               ggTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCT‑‑TGCTTCCTGGGaagcagaa                                                                                                                                           <  1:293471/111‑4 (MQ=255)
                                               ggTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCT‑‑TGCTTCCTGGGaagcagaa                                                                                                                                           >  2:293471/1‑108 (MQ=255)
                                                  ccAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCt                                                                                                    >  1:402099/1‑149 (MQ=255)
                                                     gTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTc                                                                                                 >  3:333755/1‑149 (MQ=255)
                                                      tttGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTcg                                                                                                >  2:500922/1‑149 (MQ=255)
                                                             ggTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCaga                                                                                         <  1:390928/149‑1 (MQ=255)
                                                               tAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGt                                                                                       <  4:54738/149‑1 (MQ=255)
                                                                aaGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTc                                                                                      <  2:402099/149‑1 (MQ=255)
                                                                   cGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGt                                                                                   <  1:133609/149‑1 (MQ=255)
                                                                     aCCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCa                                                                                 >  1:382499/1‑149 (MQ=255)
                                                                           gATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAgg                                                                           <  4:132209/149‑1 (MQ=255)
                                                                               ttgagattgagTAACTGCTGTCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCAcc                                                                       >  2:500543/1‑149 (MQ=255)
                                                                                tgagattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCa                                                                      <  2:364130/149‑1 (MQ=255)
                                                                                 gagattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAg                                                                     <  1:392225/149‑1 (MQ=255)
                                                                                 gagattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAg                                                                     >  3:243501/1‑149 (MQ=255)
                                                                                   gattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGTTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGc                                                                   >  1:516477/1‑149 (MQ=255)
                                                                                   gattgagTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGc                                                                   >  1:287127/1‑149 (MQ=255)
                                                                                         gTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAgcgc                                                                                                           >  2:599738/1‑103 (MQ=255)
                                                                                         gTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAgcgc                                                                                                           <  1:599738/103‑1 (MQ=255)
                                                                                           aaCTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCagag                                                           <  3:113131/149‑1 (MQ=255)
                                                                                                 tgGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGATGCGATc                                                     <  2:625175/149‑1 (MQ=255)
                                                                                                     gcACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCgcgc                                                 <  1:500922/149‑1 (MQ=255)
                                                                                                               aaTCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGcc                                                                                                                                     >  4:86976/1‑55 (MQ=255)
                                                                                                               aaTCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGcc                                                                                                                                     <  3:86976/55‑1 (MQ=255)
                                                                                                                   gTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTgtgt                                   >  4:550379/1‑149 (MQ=255)
                                                                                                                   gTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTgtgt                                   >  3:481593/1‑149 (MQ=255)
                                                                                                                           cacaTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACt                           <  4:160652/149‑1 (MQ=255)
                                                                                                                            acaTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGt                                                                                   >  2:295472/1‑92 (MQ=255)
                                                                                                                            acaTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGt                                                                                   <  1:295472/92‑1 (MQ=255)
                                                                                                                              aTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGa                        <  3:594669/149‑1 (MQ=255)
                                                                                                                              aTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGa                        <  1:500543/149‑1 (MQ=255)
                                                                                                                              aTCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGa                        >  1:158968/1‑149 (MQ=255)
                                                                                                                               tcgtcgTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGac                       <  2:158968/149‑1 (MQ=255)
                                                                                                                               tcgtcgTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGac                       <  4:360896/149‑1 (MQ=255)
                                                                                                                                cgtcgtCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGATCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGcaca                      >  3:521959/1‑149 (MQ=255)
                                                                                                                                    gtcAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACAt                  >  1:618558/1‑149 (MQ=255)
                                                                                                                                     tcAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACAtt                 >  4:95708/1‑149 (MQ=255)
                                                                                                                                     tcAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACAtt                 >  4:245668/1‑149 (MQ=255)
                                                                                                                                        gCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTc              >  3:137733/1‑149 (MQ=255)
                                                                                                                                         cTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTcc             >  4:273393/1‑149 (MQ=255)
                                                                                                                                            ttGCTTCCTGGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACg          >  4:516345/2‑149 (MQ=255)
                                                                                                                                             tGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACACTTTCCACgg         <  2:516477/149‑1 (MQ=255)
                                                                                                                                                ttCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGc      >  4:425914/1‑149 (MQ=255)
                                                                                                                                                ttCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGc      >  4:546679/1‑149 (MQ=255)
                                                                                                                                                ttCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGc      >  2:613071/1‑149 (MQ=255)
                                                                                                                                                ttCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGc      >  3:259099/1‑149 (MQ=255)
                                                                                                                                                ttCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGGCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGc      >  2:415754/1‑149 (MQ=255)
                                                                                                                                                   ccAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCaa   >  3:398927/1‑149 (MQ=255)
                                                                                                                                                    cAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAAc  >  3:33931/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCCAGGAAGCAAGCTGTGCCATGACATGTTGGTAGTCTTTAAGCGCTTTAGCTTTCGCCGCAGAGTCGGTCAGCAAGGCACCAGGCCACCAGAGGCGATCGCGCAGTTGCGGTTGTGTCACCAACTGGACAACATTTTCCACGGGGCCAAC  >  NC_000913/1048281‑1048577

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