Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 181,276 C→G 5.6% 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 5.6% 512.6 / 16.4 143H131Q (CAC→CAG) ‡degPserine endoprotease (protease Do), membrane‑associated
Reads supporting (aligned to +/- strand):  ref base C (70/65);  new base G (4/4);  total (74/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

Only 100 of 143 total aligned reads displayed.

TGCCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGC  >  NC_000913/181130‑181419
                                                                                                                                                  |                                                                                                                                               
   cAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAATCACgttgt                                                                                                                                            <  4:84775/149‑1 (MQ=255)
    aGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgtt                                                                                                                                           >  2:1447624/1‑149 (MQ=255)
     ggAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGttgttg                                                                                                                                          >  3:116110/1‑149 (MQ=255)
       aaGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttggt                                                                                                                                        <  4:465321/149‑3 (MQ=255)
                 gTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACgg                                                                                                                              <  2:1368512/149‑1 (MQ=255)
                  ttCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                            >  1:141343/1‑134 (MQ=255)
                     cAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAt                                                                                                                          <  1:1447624/149‑1 (MQ=255)
                      agagCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAtt                                                                                                                         >  3:1464829/1‑149 (MQ=255)
                      agagCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCAtt                                                                                                                         >  4:545932/1‑149 (MQ=255)
                             tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttg                                                                                                                                          <  2:550052/125‑1 (MQ=255)
                             tcCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttg                                                                                                                                          >  1:550052/1‑125 (MQ=255)
                                 ttCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACt                                                                                                              >  3:1358112/1‑149 (MQ=255)
                                        aGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGAt                                                                                                       <  1:1205639/149‑1 (MQ=255)
                                        aGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGAt                                                                                                       <  4:116110/149‑1 (MQ=255)
                                         gggTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATg                                                                                                      <  4:1358112/149‑1 (MQ=255)
                                         gggTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATg                                                                                                      >  3:1073998/1‑149 (MQ=255)
                                           gTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTaaa                                                                                                                      >  2:612877/1‑131 (MQ=255)
                                              gCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                            <  3:688462/106‑1 (MQ=255)
                                              gCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                            >  4:688462/1‑106 (MQ=255)
                                                cAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                            <  2:48683/104‑1 (MQ=255)
                                                cAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACgttgt                                                                                                                                            >  1:48683/1‑104 (MQ=255)
                                                cAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTaa                                                                                               <  4:1073998/149‑1 (MQ=255)
                                                  gggCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGAt                                                                                                       <  2:238305/139‑1 (MQ=255)
                                                  gggCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGAt                                                                                                       >  1:238305/1‑139 (MQ=255)
                                                   ggCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTaa                                                                                               >  4:1375525/1‑146 (MQ=255)
                                                   ggCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTaa                                                                                               <  3:1375525/146‑1 (MQ=255)
                                                                gcggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGAt                                                                                                       >  1:117931/1‑125 (MQ=255)
                                                                gcggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGAt                                                                                                       <  2:117931/125‑1 (MQ=255)
                                                                gcggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATg                                                                               <  2:23092/149‑1 (MQ=255)
                                                                  ggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                                                             >  4:130979/1‑149 (MQ=255)
                                                                  ggcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                                                             >  2:917635/1‑149 (MQ=255)
                                                                   gcCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGtt                                                                            <  3:545932/149‑1 (MQ=255)
                                                                     cAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTgg                                                                          >  4:330402/1‑149 (MQ=255)
                                                                      aGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTaaa                                                                                                                      >  4:320275/1‑104 (MQ=255)
                                                                      aGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTaaa                                                                                                                      <  3:320275/104‑1 (MQ=255)
                                                                      aGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGc                                                                         >  1:841412/1‑149 (MQ=255)
                                                                      aGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTAGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGc                                                                         <  4:1453763/149‑1 (MQ=255)
                                                                         aaCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCaaa                                                                      >  4:822506/1‑149 (MQ=255)
                                                                         aaCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCaaa                                                                      <  3:130979/149‑1 (MQ=255)
                                                                          aCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttgg                                                                                                                                         <  2:71059/81‑2 (MQ=255)
                                                                          aCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttgg                                                                                                                                         >  1:71059/1‑80 (MQ=255)
                                                                          aCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATg                                                                               <  1:89176/139‑1 (MQ=255)
                                                                          aCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATg                                                                               >  2:89176/1‑139 (MQ=255)
                                                                          aCAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAg                                                                     <  2:841412/149‑1 (MQ=255)
                                                                            aGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGAt                                                                   >  3:1319884/1‑149 (MQ=255)
                                                                            aGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGAt                                                                   >  1:756518/1‑149 (MQ=255)
                                                                             gAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATg                                                                               <  2:145023/136‑1 (MQ=255)
                                                                             gAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATg                                                                               >  1:145023/1‑136 (MQ=255)
                                                                                       gcgcTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTAAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGa                                                        <  4:1464829/149‑1 (MQ=255)
                                                                                           tGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATcc                                                                 >  3:1112293/1‑136 (MQ=255)
                                                                                           tGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATcc                                                                 <  4:1112293/136‑1 (MQ=255)
                                                                                               ttCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGAcgc                                                                                      >  1:248902/1‑111 (MQ=255)
                                                                                                    gCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATc                                           <  3:751503/149‑1 (MQ=255)
                                                                                                       tcatcaTTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGtt                                                                            <  2:608546/113‑1 (MQ=255)
                                                                                                       tcatcaTTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGtt                                                                            >  1:608546/1‑113 (MQ=255)
                                                                                                       tcatcaTTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCaa                                        >  3:160833/1‑149 (MQ=255)
                                                                                                         atcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAAt                                      <  2:756518/149‑1 (MQ=255)
                                                                                                         atcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAAt                                      >  4:748151/1‑149 (MQ=255)
                                                                                                         atcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAAt                                      >  4:162573/1‑149 (MQ=255)
                                                                                                          tcatTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGAcgc                                                                                      >  4:1322149/1‑100 (MQ=255)
                                                                                                           catTGATGCCGATAAAGGCTATGTCGTCACCAACAACGTGGttgttg                                                                                                                                          <  1:9730/47‑1 (MQ=37)
                                                                                                            atTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCa                                   >  3:38693/1‑149 (MQ=255)
                                                                                                             ttGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGc                                                                         <  1:87445/110‑1 (MQ=255)
                                                                                                             ttGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGc                                                                         >  2:87445/1‑110 (MQ=255)
                                                                                                             ttGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATcc                                                                 >  2:1277289/1‑118 (MQ=255)
                                                                                                             ttGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATcc                                                                 <  1:1277289/118‑1 (MQ=255)
                                                                                                                aTGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAAc                               <  4:1319884/149‑1 (MQ=255)
                                                                                                                 tGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAAt                                      >  1:971774/1‑141 (MQ=255)
                                                                                                                 tGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAAt                                      <  2:971774/141‑1 (MQ=255)
                                                                                                                      aTAAAGGGTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGCTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGaaa                         >  4:919980/1‑149 (MQ=255)
                                                                                                                      aTAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGctct                                                          >  1:1159153/1‑116 (MQ=255)
                                                                                                                      aTAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGctct                                                          <  2:1159153/116‑1 (MQ=255)
                                                                                                                       tAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGaaaa                        <  4:38693/149‑1 (MQ=255)
                                                                                                                        aaaGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGaaaaa                       <  3:535152/149‑1 (MQ=255)
                                                                                                                         aaGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAAc                      <  3:822506/149‑1 (MQ=255)
                                                                                                                         aaGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAAc                      <  4:160833/149‑1 (MQ=255)
                                                                                                                          aGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAAcc                     >  3:331894/1‑149 (MQ=255)
                                                                                                                              tATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTg                                                         >  3:355363/1‑109 (MQ=255)
                                                                                                                              tATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGAc                 >  3:249360/1‑149 (MQ=255)
                                                                                                                                 gtcgtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCaaa                                       >  4:909792/1‑124 (MQ=255)
                                                                                                                                 gtcgtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCaaa                                       <  3:909792/124‑1 (MQ=255)
                                                                                                                                  tcgtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCa             >  4:16754/1‑149 (MQ=255)
                                                                                                                                    gtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGa                  >  1:1278564/1‑142 (MQ=255)
                                                                                                                                    gtcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGa                  <  2:1278564/142‑1 (MQ=255)
                                                                                                                                     tcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGa                                                                                    <  1:480226/75‑1 (MQ=255)
                                                                                                                                     tcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGa                                                                                    >  2:480226/1‑75 (MQ=255)
                                                                                                                                     tcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCa                                   >  1:911237/1‑124 (MQ=255)
                                                                                                                                     tcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCa                                   <  3:1442986/124‑1 (MQ=255)
                                                                                                                                     tcACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAAtt          >  3:1005125/1‑149 (MQ=255)
                                                                                                                                       aCCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGa                  >  1:1212952/1‑139 (MQ=255)
                                                                                                                                        tcaacaaCGTGGTTGTTGGTAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAAtt          >  1:116755/2‑146 (MQ=255)
                                                                                                                                        ccAACAACCACGTTGTTGATAAGGCGACGGTCATTAAAGTTCAAATGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                                                             >  2:532605/1‑79 (MQ=255)
                                                                                                                                        ccAACAACCACGTTGTTGATAACGCGACGGTCTTTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGt                                                                             <  1:532605/79‑1 (MQ=255)
                                                                                                                                         caacaaCCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGa      >  2:419979/1‑149 (MQ=255)
                                                                                                                                          aacaacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCgcgc                                                <  3:718919/106‑1 (MQ=255)
                                                                                                                                             aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGa                  >  4:625691/1‑133 (MQ=255)
                                                                                                                                             aacCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGc  >  1:826421/1‑149 (MQ=255)
                                                                                                                                                cACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATcc                                          >  4:1218847/1‑106 (MQ=255)
                                                                                                                                                cACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATcc                                          <  3:1218847/106‑1 (MQ=255)
                                                                                                                                                cACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGa                  >  4:749295/1‑130 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TGCCAGGAAGGTTCTCCGTTCCAGAGCTCTCCGTTCTGCCAGGGTGGCCAGGGCGGTAATGGTGGCGGCCAGCAACAGAAATTCATGGCGCTGGGTTCCGGCGTCATCATTGATGCCGATAAAGGCTATGTCGTCACCAACAACCACGTTGTTGATAACGCGACGGTCATTAAAGTTCAACTGAGCGATGGCCGTAAGTTCGACGCGAAGATGGTTGGCAAAGATCCGCGCTCTGATATCGCGCTGATCCAAATCCAGAACCCGAAAAACCTGACCGCAATTAAGATGGC  >  NC_000913/181130‑181419

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