Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 181,274 C→G 6.1% H131D (CAC→GAC) ‡ degP → serine endoprotease (protease Do), membrane‑associated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913181,2740CG 6.1% 211.4 / 5.1 77H131D (CAC→GAC) ‡degPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base C (46/23);  new base G (2/3);  total (50/27)
Fisher's exact test for biased strand distribution p-value = 3.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

GTTCTGCCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCG  >  NC_000913/181126‑181420
                                                                                                                                                    |                                                                                                                                                  
gagcTGCCAGGAAGGTTCTCCGTTCCAGAGCTCTCGGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCAAGGCGCTGGGTTCCGGCGGCATCATTGATGCCGATAAAGGCTATGTCGTCACCAAGAAcc                                                                                                                                                    <  2:1406485/146‑1 (MQ=255)
      ccAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttg                                                                                                                                              >  1:1520445/1‑149 (MQ=255)
       cAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGCTGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                             >  2:56460/1‑149 (MQ=255)
         ggAATGTTCTCCGGTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGttgttg                                                                                                                                           <  2:1250336/149‑1 (MQ=255)
         ggAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGttgttg                                                                                                                                           >  2:413595/1‑149 (MQ=255)
           aaGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttggt                                                                                                                                         <  1:448563/149‑3 (MQ=255)
                tctcCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAAcgc                                                                                                                                    >  2:477365/1‑149 (MQ=255)
                 ctcCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGTTGTTGGTGAcgac                                                                                                                                   <  2:689029/149‑3 (MQ=255)
                  tcCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGa                                                                                                                                  >  2:987331/1‑149 (MQ=255)
                   ccGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGGTGTTGATAACGCGAc                                                                                                                                 >  1:688664/1‑149 (MQ=255)
                    cGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACg                                                                                                                                >  1:1204700/1‑149 (MQ=255)
                    cGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACg                                                                                                                                >  2:761212/1‑149 (MQ=255)
                    cGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACg                                                                                                                                >  2:1217523/1‑149 (MQ=255)
                    cGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTACGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACACCAACGTTGTTGATAACGCGACg                                                                                                                                >  1:805606/1‑149 (MQ=255)
                     gTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACgg                                                                                                                               >  2:1022381/1‑149 (MQ=255)
                       tCCAGAGCGCTGCGTTCTGCCAGGGTCGCCAGGGCGGAAATCGGGGCGTCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGGCATCATTGATGCGGATAAAGGCTATGGCGTCACCAACAACCACGTTGTTGATAACGCGACGGTc                                                                                                                             <  2:78222/149‑1 (MQ=255)
                             gCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTaaa                                                                                                                       <  2:810501/149‑1 (MQ=255)
                              ctctCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTTATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAg                                                                                                                      <  2:1460724/149‑1 (MQ=255)
                              ctctCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTCGGTTCCGGCGTCAGCATTGATGCCGATAAAGTCTATGTCGTCACCAACAACCACGGTGTTGATACAGCGACGGTCATTaaat                                                                                                                      >  1:1026312/1‑148 (MQ=255)
                               tctcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCAGGTTGTTGATAACGCGACGGGCATTAAAGt                                                                                                                     >  1:394083/1‑149 (MQ=255)
                                ctcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGtt                                                                                                                    >  1:1113023/1‑149 (MQ=255)
                                 tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGATATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTc                                                                                                                   <  1:661413/149‑1 (MQ=255)
                                 tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTc                                                                                                                   >  1:209246/1‑149 (MQ=255)
                                 tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTCGCGTCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCTCGTTGTTGATAACGCGACGGTCATTAAAGTTc                                                                                                                   >  1:874600/1‑149 (MQ=255)
                                   cGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCaa                                                                                                                 <  2:142990/149‑1 (MQ=255)
                                   cGTGCTGCCAGGGTGTCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTACGGCGGCATCATTGATGCCGCTAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCaa                                                                                                                 <  2:148453/149‑1 (MQ=255)
                                          ccAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCg                                                                                                          >  1:256251/1‑149 (MQ=255)
                                             gggTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATg                                                                                                       <  1:814401/149‑1 (MQ=255)
                                                    cAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGCCGGTCATTAAAGTTCAACTGAGCGATGGCCGTaa                                                                                                >  1:1184809/1‑149 (MQ=255)
                                                    cAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACCTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTaa                                                                                                >  2:24426/1‑149 (MQ=255)
                                                    cAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGGTCCGGCGACATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACAACGTTGTTGATAACGCGCCGGTCATTAAAGTTCAACTGAGCGATGGCCGTaa                                                                                                >  1:1188200/1‑149 (MQ=255)
                                                     aGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCTCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAg                                                                                               <  1:299374/149‑1 (MQ=255)
                                                      gggCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGt                                                                                              >  1:624086/1‑149 (MQ=255)
                                                           gTAATGGTGGCGGACAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGAc                                                                                         <  2:1007636/149‑1 (MQ=255)
                                                                             aaCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCaaa                                                                       >  2:27998/1‑149 (MQ=255)
                                                                              aCAGAAATACATGGCGCTGGGTTCCGGCGTCATCATTGAAGCCGATAAAGGCTATGTCTTCACCAACAACAACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGTGCTATGGCCGTAAGTTCGACGCGTAGATGGTTGGCAAAg                                                                      <  2:452302/149‑1 (MQ=255)
                                                                                aGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCCACAACCACGTTGTTGATAACGCGCCGGTCATTAAAGTTCACCGGAGCGATGGCCATAAGTTCGACGCGAAGATGGTTGGCAAAGAt                                                                    >  1:398282/1‑149 (MQ=255)
                                                                                aGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGAt                                                                    >  1:1508714/1‑149 (MQ=255)
                                                                                  aaaTTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACACACACGTTGTTGATACCGCGACGGTCATTAAAGTTCAACTTAGCGATGGCCGTAAGGTCGACGCGATGATGGTTGGCAAAGATcc                                                                  >  1:1033880/1‑149 (MQ=255)
                                                                                  aaaTTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATcc                                                                  <  1:761212/149‑1 (MQ=255)
                                                                                   aaTTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCcg                                                                 <  2:22547/149‑1 (MQ=255)
                                                                                    aTTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCcgc                                                                <  1:1235209/149‑1 (MQ=255)
                                                                                      tCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGAGGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGGTGGCAAAGATCCgcgc                                                              >  2:106632/1‑149 (MQ=255)
                                                                                             gcTGGGTTCCGGCGTCATCATTGATGCCGATAAAGTCTATGTCGTCAACAACAACGACGTTGTTGATTACGCTACGGTCATTAAAGTTCAACTGAGCGATTGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGata                                                       <  1:46742/149‑1 (MQ=255)
                                                                                               tGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAATTACAACTGAGCGATGGCAGTAAGTTCGAAGCGAAGATTGTTGGCAAAGATCCGCGCTCTGatata                                                     >  1:856229/1‑148 (MQ=255)
                                                                                                gggTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATcg                                                    <  1:1538501/149‑1 (MQ=255)
                                                                                                  gTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATcgcg                                                  <  2:1520445/149‑1 (MQ=255)
                                                                                                       ggcggcATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGAt                                             <  2:552445/149‑1 (MQ=255)
                                                                                                       ggCGTCATCATTGATGCCGATAAATGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGAt                                             <  1:322330/149‑1 (MQ=255)
                                                                                                       ggCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCCCGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGAt                                             <  1:241490/149‑1 (MQ=255)
                                                                                                          gTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCa                                          >  1:870527/1‑149 (MQ=255)
                                                                                                          gTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCa                                          >  2:955458/1‑149 (MQ=255)
                                                                                                          gTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTAGACGCGAAGATGGTTGGCAAAGATCCTCGATCTTATATCGCTCTGATCCa                                          >  2:371427/1‑149 (MQ=255)
                                                                                                             atcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAAt                                       <  1:1022381/149‑1 (MQ=255)
                                                                                                              tcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATc                                      >  1:1237735/1‑149 (MQ=255)
                                                                                                              tcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTAGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCAGCGCTCTGATATCGAGCTGATCAAAATc                                      >  1:306764/1‑149 (MQ=255)
                                                                                                               catTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATcc                                     >  1:837242/1‑149 (MQ=255)
                                                                                                               catTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATcc                                     >  2:775568/1‑149 (MQ=255)
                                                                                                               catTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATcc                                     <  2:533828/149‑1 (MQ=255)
                                                                                                                 ttGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAg                                   >  1:1026750/1‑149 (MQ=255)
                                                                                                                    aTGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCTCGCTGATCCAAATACAGAAc                                >  2:1466271/1‑149 (MQ=255)
                                                                                                                          aTAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATACCGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGaaa                          >  1:157165/1‑149 (MQ=255)
                                                                                                                            aaaGGCTATGTCGTCACCAACACCCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGAGCCAAATCCAGAACCCGaaaaa                        >  1:1556208/1‑149 (MQ=255)
                                                                                                                            aaaGGCTATGTCGTCACCAACAACCACGTTGTGGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGaaaaa                        <  1:710581/149‑1 (MQ=255)
                                                                                                                             aaGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAAc                       <  2:209246/149‑1 (MQ=255)
                                                                                                                               ggCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCt                     >  1:1469082/1‑149 (MQ=255)
                                                                                                                                  tATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGAc                  >  1:431011/1‑149 (MQ=255)
                                                                                                                                   aTGTCGTCACCATCAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCACCTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATACTCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGAcc                 >  2:302207/1‑149 (MQ=255)
                                                                                                                                    tGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCg                >  1:1175903/1‑149 (MQ=255)
                                                                                                                                      tcgtcACCAACAACCACGTTGTTGATAAAGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCTACGCGAAGAAGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCa              >  2:762787/1‑149 (MQ=255)
                                                                                                                                       cgttatcaacaaCGTGGTTGTTGGTGACGCGCCGGTCATTAAAGTTCACCTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAACTCCAGCCCCCGAAAGACCTGACCGCaa             >  1:957855/7‑149 (MQ=255)
                                                                                                                                        gtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAAt            >  2:742334/1‑149 (MQ=255)
                                                                                                                                         ttatcaacaaCGTGGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTTGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGACCCCGAAACACCTGCCCGCAAtt           >  1:625533/5‑149 (MQ=255)
                                                                                                                                           aCCAACAACCACGTTGATTATAACCCGACGGTCATTAAAGATCAACTGAGCGATGGCAGTAAGTTCGACGCGAAGAAGGATGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTaa         >  2:949403/1‑149 (MQ=255)
                                                                                                                                            ccAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAg        >  1:940738/1‑149 (MQ=255)
                                                                                                                                             caacaaCCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGa       >  1:641860/1‑149 (MQ=255)
                                                                                                                                                 aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGc   >  2:1244779/1‑149 (MQ=255)
                                                                                                                                                  acCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCg  <  2:624086/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTTCTGCCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCG  >  NC_000913/181126‑181420

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: ‑