Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,621 A→C 5.5% A426A (GCT→GCG 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,6210AC 5.5% 362.0 / 9.5 127A426A (GCT→GCGcysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (69/50);  new base C (3/4);  total (72/55)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.06e-01

Only 100 of 127 total aligned reads displayed.

CGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAA  >  NC_000913/2890473‑2890769
                                                                                                                                                    |                                                                                                                                                    
cggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAg                                                                                                                                                     <  2:674115/150‑1 (MQ=255)
  gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCg                                                                                                                                                   >  1:1511956/1‑150 (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)
              aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                       >  1:1651305/1‑150 (MQ=255)
              aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                       <  2:537120/150‑1 (MQ=255)
               tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaagaggg                                                                                                                                      <  2:1837386/150‑8 (MQ=255)
               tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAgg                                                                                                                                      >  2:675813/1‑150 (MQ=255)
                aaaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTGCTGGATTGGCTGGCATGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGc                                                                                                                                     <  2:290889/150‑1 (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)
                        ggAACGGCGCAATACCGGTGCCTGGGCCAATCATAAGCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCgg                                                                                                                             <  1:1044286/150‑1 (MQ=255)
                          aaCGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCa                                                                                                                           <  2:4821/150‑1 (MQ=255)
                             ggCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACgg                                                                                                                        >  1:735081/1‑150 (MQ=255)
                              gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACggg                                                                                                                       >  1:1908688/1‑150 (MQ=255)
                              gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACggg                                                                                                                       >  1:1707931/1‑150 (MQ=255)
                               cgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCCCGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGgc                                                                                                                      >  1:1802468/1‑150 (MQ=255)
                                  aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGcgcg                                                                                                                   >  2:445125/1‑150 (MQ=255)
                                   aTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCgg                                                                                                                  <  2:696097/150‑1 (MQ=255)
                                   aTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCgg                                                                                                                  <  1:1618313/150‑1 (MQ=255)
                                       cgctGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCtt                                                                                                              <  1:1533854/147‑1 (MQ=255)
                                           gCCTGGGCCACTCAACATCACCGGGGTTTCTGGCTTGGCTGGCAGGCGGAAGATATCGTTATGTTCGATAAATCCGCGGACTTCGCCCTCTTCTTCCACGCGCTCTGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAc                                                                                                          <  1:78776/150‑1 (MQ=255)
                                                ggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtcgt                                                                                                     >  2:30891/1‑150 (MQ=255)
                                                 gCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGCGGCCCCCCCGGCCCGGGCGCGGCCTTCCACGTCGTa                                                                                                    >  1:1434756/1‑150 (MQ=255)
                                                 gCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTa                                                                                                    >  1:1130361/1‑150 (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)
                                                          aaTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccac                                                                                           <  1:1764582/150‑1 (MQ=255)
                                                             cACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacacc                                                                                        <  2:1059794/150‑1 (MQ=255)
                                                                 ggggTTTCTGGATTGGTTGGCAGGGGGAAGTTATCGTTATGTTCTATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGGACCACACCAACg                                                                                    <  2:1249303/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)
                                                                                      aGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGGGTACTTCGTTCTCg                                                               <  2:1319268/150‑1 (MQ=255)
                                                                                      aGGCGGAAGGTATCGTTATGATCGATAAATACGCTGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCg                                                               <  2:330459/150‑1 (MQ=255)
                                                                                       ggCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCgc                                                              >  1:871509/1‑149 (MQ=255)
                                                                                        gCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGAc                                                             <  2:1731857/150‑1 (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)
                                                                                                    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)
                                                                                                                    aCGCGGACTTCGCTCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAc                                 >  2:772719/1‑150 (MQ=255)
                                                                                                                    aCGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAc                                 <  1:1747846/150‑1 (MQ=255)
                                                                                                                     cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                <  2:1511956/150‑1 (MQ=255)
                                                                                                                     cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCCCGTCGTACCGCACCACACCACCGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATAcc                                >  1:140546/1‑149 (MQ=255)
                                                                                                                     cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                <  2:1156716/150‑1 (MQ=255)
                                                                                                                     cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                >  1:638442/1‑150 (MQ=255)
                                                                                                                     cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                                >  1:1914948/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)
                                                                                                                          aCTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTGTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc                           >  1:1106455/1‑150 (MQ=255)
                                                                                                                           cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTACGCCTGCGAGGAGGCCTAGGAATACAGAcgcg                          >  2:195125/1‑150 (MQ=255)
                                                                                                                             tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGCTGGAAGACAGACgcggc                        >  1:317387/1‑150 (MQ=255)
                                                                                                                             tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggc                        <  2:960487/150‑1 (MQ=255)
                                                                                                                             tCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggc                        >  2:74439/1‑150 (MQ=255)
                                                                                                                               gCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTACCGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGCCGCGGCGt                      >  2:392065/1‑150 (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)
                                                                                                                                     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)
                                                                                                                                             cgcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTTTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcagca        <  2:1881745/150‑1 (MQ=255)
                                                                                                                                             cgcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcagca        >  2:654203/1‑150 (MQ=255)
                                                                                                                                              gcgGTCAGCGAGGAAGCTGGAGGCACCTCCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGcagcag       <  1:636418/150‑1 (MQ=255)
                                                                                                                                              gcgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGAGAGGAGGCGATGGAATACAGACGCGGCGTCATCGGGCGcagcag       <  2:1620266/150‑1 (MQ=255)
                                                                                                                                               cgGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGa      >  2:1068641/1‑150 (MQ=255)
                                                                                                                                                   cAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGattaa  >  2:518771/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAA  >  NC_000913/2890473‑2890769

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: