Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 181,276 C→G 6.3% H131Q (CAC→CAG) ‡ 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,2760CG 6.3% 203.7 / 5.3 77H131Q (CAC→CAG) ‡degPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base C (47/25);  new base G (3/2);  total (50/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

CCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCG  >  NC_000913/181132‑181420
                                                                                                                                                |                                                                                                                                                
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                                                                                                                                >  2:761212/1‑149 (MQ=255)
              cGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACg                                                                                                                                >  1:1204700/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                                     >  2:775568/1‑149 (MQ=255)
                                                                                                         catTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATcc                                     <  2:533828/149‑1 (MQ=255)
                                                                                                         catTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATcc                                     >  1:837242/1‑149 (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)
                                                                                                                                                |                                                                                                                                                
CCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGCG  >  NC_000913/181132‑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: