Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,613 A→T 5.7% V429E (GTG→GAG)  cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,6130AT 5.7% 317.8 / 8.7 124V429E (GTG→GAG) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (64/50);  new base T (3/4);  total (70/54)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-02

Only 100 of 125 total aligned reads displayed.

CGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGC  >  NC_000913/2890467‑2890762
                                                                                                                                                  |                                                                                                                                                     
cgtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgc                                                                                                                                                    >  2:220086/1‑150 (MQ=255)
  tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGGGCCTGGGCCAATCATACTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCgg                                                                                                                                                  <  2:306027/150‑1 (MQ=255)
   cggcggCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCAGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGt                                                                                                                                                 <  2:1527357/150‑1 (MQ=255)
      cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAg                                                                                                                                              <  2:674115/150‑1 (MQ=255)
          gcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAg                                                                                                                                          <  1:1055874/150‑1 (MQ=255)
               gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGCCTTCGCCCTCTTCTTCCCCGCGGTCAGCGAGGAAGc                                                                                                                                     >  1:498764/1‑150 (MQ=255)
               gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAga                                                                                                                                     <  1:1574710/150‑2 (MQ=255)
               gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGc                                                                                                                                     >  2:92940/1‑150 (MQ=255)
                ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGCGGTATCTGGATTGGCTGGCAGGCGGAAGTAAACGTTATCTTCGCTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGCTCAGCGAGGAAGCt                                                                                                                                    >  1:1249303/1‑150 (MQ=255)
                 tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTg                                                                                                                                   >  1:1958677/1‑150 (MQ=255)
                  gcATAACGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTgg                                                                                                                                  <  2:478384/150‑1 (MQ=255)
                  gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCCCGCGGTCAGCGAGGAAGCTgg                                                                                                                                  >  2:1204870/1‑150 (MQ=255)
                   cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTTTCGTGAGGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGa                                                                                                                                 >  1:1255466/1‑150 (MQ=255)
                   cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaagag                                                                                                                                 <  1:754218/150‑6 (MQ=255)
                   cATAAAGGCGCGGAACGGCGCAATACCGGTGCCAGGGCCAATCATAAAAACCGGGTTTTCTGTATTGGCTTGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGAGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGTAAGCTGGa                                                                                                                                 >  2:183740/1‑150 (MQ=255)
                   aaTAAAGGCGCGGAACTGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaagag                                                                                                                                 <  2:1341281/149‑6 (MQ=255)
                    atatAGGCGCGGAACGGCGCAATACCGGTGCCTGTGCCAATCATAATCACCGTGTTTTCTGGCTTGGCTGGCAGGAGTAAGTTATCTTTATGTTCTATAAATACGAGGACTTCGCCATCTTCTTCCACGCGTTCAGCGAGGAAGCTGGAg                                                                                                                                <  1:1091705/150‑1 (MQ=255)
                    aaaaaGGCGCGGAACTGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                <  1:1893939/148‑1 (MQ=255)
                    aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATACTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                >  1:1506962/1‑150 (MQ=255)
                    aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                <  2:537120/150‑1 (MQ=255)
                     tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAgg                                                                                                                               <  2:1705134/150‑1 (MQ=255)
                     tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAgg                                                                                                                               >  2:675813/1‑150 (MQ=255)
                       aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGca                                                                                                                             <  2:1418188/150‑1 (MQ=255)
                       aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATCGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGca                                                                                                                             >  2:1747846/1‑150 (MQ=255)
                         ggCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGcacc                                                                                                                           <  1:1005335/150‑1 (MQ=255)
                           cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGcaccac                                                                                                                         >  2:1130866/1‑150 (MQ=255)
                           cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGcaccac                                                                                                                         >  1:26595/1‑150 (MQ=255)
                              ggAACTGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCgg                                                                                                                      <  2:1221039/150‑1 (MQ=255)
                              ggAACGGCGCAATACCGGTGCCTGGGCCAATCATAAGCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCgg                                                                                                                      <  1:1044286/150‑1 (MQ=255)
                                aaCGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCa                                                                                                                    <  2:4821/150‑1 (MQ=255)
                                  cGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACg                                                                                                                  <  2:913754/150‑1 (MQ=255)
                                    gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACggg                                                                                                                >  1:1707931/1‑150 (MQ=255)
                                     cgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCCCGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGgc                                                                                                               >  1:1802468/1‑150 (MQ=255)
                                        aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGcgcg                                                                                                            >  2:445125/1‑150 (MQ=255)
                                         aTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCgg                                                                                                           <  1:1618313/150‑1 (MQ=255)
                                            ccGGTGCCGGGGCCAATCATAATCACCTGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCTATAAATACGCGGACTTCGCCCTCTTCATCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCt                                                                                                        <  2:478603/150‑1 (MQ=255)
                                             cgctGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCtt                                                                                                       <  1:1533854/147‑1 (MQ=255)
                                             cGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGCGCGCGGCCtt                                                                                                       >  2:579167/1‑150 (MQ=255)
                                               gTGCCTGGTCCAATCATACTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTcc                                                                                                     <  1:1761944/150‑1 (MQ=255)
                                                 gCCTGGGCCACTCAACATCACCGGGGTTTCTGGCTTGGCTGGCAGGCGGAAGATATCGTTATGTTCGATAAATCCGCGGACTTCGCCCTCTTCTTCCACGCGCTCTGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAc                                                                                                   <  1:78776/150‑1 (MQ=255)
                                                      ggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtcgt                                                                                              >  1:1889981/1‑150 (MQ=255)
                                                      ggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtcgt                                                                                              >  2:30891/1‑150 (MQ=255)
                                                       gCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTa                                                                                             >  1:1130361/1‑150 (MQ=255)
                                                         cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                           <  1:459100/150‑1 (MQ=255)
                                                         cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGGTCGATAAATACGCGGACTTCGCCCGCTTCTTCCACGCGGTCACCGAGGAAGCTGGAGGCACCCCCGGCCCCGGCGCGGCCTTCCACGTCGTAAc                                                                                           >  2:1264165/1‑150 (MQ=255)
                                                           aTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTTGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGc                                                                                         >  2:1990069/1‑150 (MQ=255)
                                                             cATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                       >  1:1497428/1‑150 (MQ=255)
                                                               tAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcacca                                                                                     >  2:1426001/1‑150 (MQ=255)
                                                                aaTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccac                                                                                    <  1:1764582/150‑1 (MQ=255)
                                                                   cACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacacc                                                                                 <  2:1059794/150‑1 (MQ=255)
                                                                    aCCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacacca                                                                                >  2:1108290/1‑150 (MQ=255)
                                                                       ggggTTTCTGGATTGGTTGGCAGGGGGAAGTTATCGTTATGTTCTATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGGACCACACCAACg                                                                             <  2:1249303/150‑1 (MQ=255)
                                                                          gTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTg                                                                          >  1:1620266/1‑150 (MQ=255)
                                                                               tGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGCGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGCTGACgtg                                                                     >  1:879376/1‑150 (MQ=255)
                                                                               tGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACgtg                                                                     >  1:1982347/1‑150 (MQ=255)
                                                                                      gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCg                                                              <  2:1760031/150‑1 (MQ=255)
                                                                                         ggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtc                                                           >  1:1850932/1‑150 (MQ=255)
                                                                                            aGGCGGAAGGTATCGTTATGATCGATAAATACGCTGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                        <  2:330459/150‑1 (MQ=255)
                                                                                             ggCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCgc                                                       >  1:871509/1‑149 (MQ=255)
                                                                                              gCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGAc                                                      <  2:1731857/150‑1 (MQ=255)
                                                                                               cGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACt                                                     >  2:1285452/1‑150 (MQ=255)
                                                                                                    gTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGc                                                >  2:1762795/1‑150 (MQ=255)
                                                                                                     ttATCGTTATGTTCGATAAATACGCGGTCTTCGCCCTCTTCTTCCACTCGGTCATCTAGGAATATGGAGTCACAACCGTCACGGGCGCGGCCTTCCACGTCGTAACGCAGCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGcc                                               >  1:1393931/1‑150 (MQ=255)
                                                                                                          gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGATGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGAGACGTGTACTTCGTTCTCGACTTCCGCCTGCgc                                          >  1:454214/1‑149 (MQ=255)
                                                                                                              tGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgaggag                                      >  2:768897/1‑150 (MQ=255)
                                                                                                               gtgcGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTACCGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                                     <  1:1010738/147‑1 (MQ=255)
                                                                                                               gatCGATAAATACGCGGACTACGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCTCGGCCTTACACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                                     >  1:1386764/3‑150 (MQ=255)
                                                                                                               gTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                                     <  2:1506962/150‑1 (MQ=255)
                                                                                                               gTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                                     >  1:485872/1‑150 (MQ=255)
                                                                                                                  cGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                  >  2:926679/1‑150 (MQ=255)
                                                                                                                       aaTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGCCTTCCGCCTGCGAGGAGGCGATGGaa                             >  1:117422/1‑150 (MQ=255)
                                                                                                                         tACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATa                           >  1:1495122/1‑150 (MQ=255)
                                                                                                                          aCGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAc                          <  1:1747846/150‑1 (MQ=255)
                                                                                                                           cgcgGCATTCGCCCTCTTCTTCCACTCGGTCAGCGATGAAGCTGGACGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACTCACCACACCAACGGTGACGTGTACTTCGTTCTCGTCTTCCGCCTGCGAGGAGGCGATGGAATACa                         <  2:309716/150‑1 (MQ=255)
                                                                                                                           cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCCCGTCGTACCGCACCACACCACCGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAcc                         >  1:140546/1‑149 (MQ=255)
                                                                                                                           cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                         >  1:638442/1‑150 (MQ=255)
                                                                                                                           cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                         >  1:1914948/1‑150 (MQ=255)
                                                                                                                           cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                         <  2:1156716/150‑1 (MQ=255)
                                                                                                                            gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAg                        >  2:730212/1‑150 (MQ=255)
                                                                                                                            gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAg                        >  1:1073136/1‑150 (MQ=255)
                                                                                                                             cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa                       <  2:1908688/150‑1 (MQ=255)
                                                                                                                               gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGTCACGGGCGCGGCCTTCCATGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcg                     <  2:1667247/150‑1 (MQ=255)
                                                                                                                               gACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGGGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcg                     <  1:1426001/150‑1 (MQ=255)
                                                                                                                                aCTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTGTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc                    >  1:1106455/1‑150 (MQ=255)
                                                                                                                                 cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTACGCCTGCGAGGAGGCCTAGGAATACAGAcgcg                   >  2:195125/1‑150 (MQ=255)
                                                                                                                                 cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgcg                   <  2:1889981/150‑1 (MQ=255)
                                                                                                                                   tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGCTGGAAGACAGACgcggc                 >  1:317387/1‑150 (MQ=255)
                                                                                                                                   tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggg                 >  1:736733/1‑149 (MQ=255)
                                                                                                                                   tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggc                 >  2:74439/1‑150 (MQ=255)
                                                                                                                                   tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggc                 <  2:960487/150‑1 (MQ=255)
                                                                                                                                      cccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCCCCGGCACGGGCGCGGCCTTCCACGTCGTACCGCACCACACCCACTGTGACGGGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTc              <  2:1255466/150‑1 (MQ=255)
                                                                                                                                       ctccagcttcCTCGCTGACCGCGTGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCCCGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTcc             >  2:1017639/7‑149 (MQ=255)
                                                                                                                                       ccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCa             <  2:26595/150‑1 (MQ=255)
                                                                                                                                         tcttctTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGc           <  2:1713979/150‑1 (MQ=255)
                                                                                                                                          cagcttcCTCGCTGACCGCGTGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCg          >  1:575470/4‑150 (MQ=255)
                                                                                                                                           ttcttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCgg         >  2:764058/1‑150 (MQ=255)
                                                                                                                                           agcttcCTCGCTGACCGCGTGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTACCGCACCACACCAACGGGGTCGTGTCCTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAAGTCCGACGCGGCGTCAGCgg         >  2:274750/3‑150 (MQ=255)
                                                                                                                                            tcttcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGACTGCGTGGAGGCGATGGTATACAGACGCGGCGTCAGCggg        >  1:1263165/1‑150 (MQ=255)
                                                                                                                                               tcCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGgcgc     >  1:1503052/1‑150 (MQ=255)
                                                                                                                                                 cACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcag   >  1:666287/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGC  >  NC_000913/2890467‑2890762

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: