Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,625 A→T 5.4% L425H (CTC→CAC)  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,6250AT 5.4% 384.5 / 8.7 130L425H (CTC→CAC) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (71/52);  new base T (3/4);  total (74/56)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 133 total aligned reads displayed.

CGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGC  >  NC_000913/2890476‑2890774
                                                                                                                                                     |                                                                                                                                                     
      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)
           aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                            >  1:1651305/1‑150 (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                                                                                                                                     >  1:26595/1‑150 (MQ=255)
                  cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGcaccac                                                                                                                                     >  2:1130866/1‑150 (MQ=255)
                     ggAACTGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCgg                                                                                                                                  <  2:1221039/150‑1 (MQ=255)
                     ggAACGGCGCAATACCGGTGCCTGGGCCAATCATAAGCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCgg                                                                                                                                  <  1:1044286/150‑1 (MQ=255)
                          ggCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACgg                                                                                                                             >  1:735081/1‑150 (MQ=255)
                           gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACggg                                                                                                                            >  1:1707931/1‑150 (MQ=255)
                               aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGcgcg                                                                                                                        >  2:445125/1‑150 (MQ=255)
                                aTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCgg                                                                                                                       <  2:696097/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)
                                             ggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtcgt                                                                                                          >  2:30891/1‑150 (MQ=255)
                                             ggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtcgt                                                                                                          >  1:1889981/1‑150 (MQ=255)
                                              gCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGCGGCCCCCCCGGCCCGGGCGCGGCCTTCCACGTCGTa                                                                                                         >  1:1434756/1‑150 (MQ=255)
                                              gCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTa                                                                                                         >  1:1130361/1‑150 (MQ=255)
                                                cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                                       <  1:459100/150‑1 (MQ=255)
                                                cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                                       <  1:497483/150‑1 (MQ=255)
                                                cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGGTCGATAAATACGCGGACTTCGCCCGCTTCTTCCACGCGGTCACCGAGGAAGCTGGAGGCACCCCCGGCCCCGGCGCGGCCTTCCACGTCGTAAc                                                                                                       >  2:1264165/1‑150 (MQ=255)
                                                    cATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                                   >  1:1497428/1‑150 (MQ=255)
                                                     aTAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcacc                                                                                                  >  1:877673/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)
                                                                 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)
                                                                             gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCg                                                                          >  1:1782667/1‑150 (MQ=255)
                                                                              ctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGt                                                                         <  2:1956529/150‑1 (MQ=255)
                                                                                ggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtc                                                                       >  1:1850932/1‑150 (MQ=255)
                                                                                ggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGGCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtc                                                                       <  2:172461/150‑1 (MQ=255)
                                                                                  cAGGCGGAAGATTTCGTTATGTTCGATAAATACGCGTACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCAGCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCTCACCAACGGTGACGTGTACTTCGTtctc                                                                     <  1:281966/150‑1 (MQ=255)
                                                                                   aGGCGGAAGGTATCGTTATGATCGATAAATACGCTGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                                    <  2:330459/150‑1 (MQ=255)
                                                                                    ggCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCgc                                                                   >  1:871509/1‑149 (MQ=255)
                                                                                      cGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACt                                                                 >  2:1285452/1‑150 (MQ=255)
                                                                                           gTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGc                                                            >  2:1762795/1‑150 (MQ=255)
                                                                                            ttATCGTTATGTTCGATAAATACGCGGTCTTCGCCCTCTTCTTCCACTCGGTCATCTAGGAATATGGAGTCACAACCGTCACGGGCGCGGCCTTCCACGTCGTAACGCAGCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGcc                                                           >  1:1393931/1‑150 (MQ=255)
                                                                                              aTCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCCCGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTg                                                         >  2:1118525/1‑150 (MQ=255)
                                                                                               tCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGc                                                        >  2:1245332/1‑150 (MQ=255)
                                                                                                     tGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgaggag                                                  >  2:768897/1‑150 (MQ=255)
                                                                                                      gtgcGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTACCGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                                                 <  1:1010738/147‑1 (MQ=255)
                                                                                                      gatCGATAAATACGCGGACTACGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCTCGGCCTTACACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGaggagg                                                 >  1:1386764/3‑150 (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)
                                                                                                                 aCGCGGACTTCGCTCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAc                                      >  2:772719/1‑150 (MQ=255)
                                                                                                                  cgcgGCATTCGCCCTCTTCTTCCACTCGGTCAGCGATGAAGCTGGACGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACTCACCACACCAACGGTGACGTGTACTTCGTTCTCGTCTTCCGCCTGCGAGGAGGCGATGGAATACa                                     <  2:309716/150‑1 (MQ=255)
                                                                                                                  cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                     <  2:1511956/150‑1 (MQ=255)
                                                                                                                  cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCCCGTCGTACCGCACCACACCACCGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAcc                                     >  1:140546/1‑149 (MQ=255)
                                                                                                                  cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                     >  1:1914948/1‑150 (MQ=255)
                                                                                                                  cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                     >  1:638442/1‑150 (MQ=255)
                                                                                                                   gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAg                                    >  1:1073136/1‑150 (MQ=255)
                                                                                                                   gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAg                                    >  2:730212/1‑150 (MQ=255)
                                                                                                                    cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa                                   <  2:1908688/150‑1 (MQ=255)
                                                                                                                        cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTACGCCTGCGAGGAGGCCTAGGAATACAGAcgcg                               >  2:195125/1‑150 (MQ=255)
                                                                                                                          tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGCTGGAAGACAGACgcggc                             >  1:317387/1‑150 (MQ=255)
                                                                                                                          tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggg                             >  1:736733/1‑149 (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)
                                                                                                                                       ccACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgca                <  1:1069844/150‑1 (MQ=255)
                                                                                                                                         aCGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcagc              <  1:1130866/150‑1 (MQ=255)
                                                                                                                                          cgcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTTTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcagca             <  2:1881745/150‑1 (MQ=255)
                                                                                                                                           gcgGTCAGCGAGGAAGCTGGAGGCACCTCCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGcagcag            <  1:636418/150‑1 (MQ=255)
                                                                                                                                           gcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGAGAGGAGGCGATGGAATACAGACGCGGCGTCATCGGGCGcagcag            <  2:1620266/150‑1 (MQ=255)
                                                                                                                                                cAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGattaa       >  2:518771/1‑150 (MQ=255)
                                                                                                                                                  gCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCAGCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGattaatt     <  2:2016672/150‑1 (MQ=255)
                                                                                                                                                  gCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACCCCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGatttatt     >  2:1038411/1‑150 (MQ=255)
                                                                                                                                                  gCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGattaatt     >  2:1512754/1‑150 (MQ=255)
                                                                                                                                                   cGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCTCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGattaatta    >  2:1846958/1‑150 (MQ=255)
                                                                                                                                                     aGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCGTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAgc  <  1:1032629/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGC  >  NC_000913/2890476‑2890774

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: