Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,521 G→A 100% P460S (CCA→TCA)  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,5210GA97.3% 527.0 / 1.0 166P460S (CCA→TCA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (1/4);  new base A (87/74);  total (88/78)
Fisher's exact test for biased strand distribution p-value = 1.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

AAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGG  >  NC_000913/2890424‑2890618
                                                                                                 |                                                                                                 
aaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGc                                                                                                 >  2:616841/1‑100 (MQ=255)
 aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                                <  2:3678102/100‑1 (MQ=255)
  gTGCGGCTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCa                                                                                               <  2:1505853/100‑1 (MQ=255)
   tGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCaa                                                                                              >  1:996865/1‑100 (MQ=255)
    gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAAt                                                                                             >  1:2910006/1‑100 (MQ=255)
    gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAAt                                                                                             >  2:570747/1‑100 (MQ=255)
    gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAAt                                                                                             >  1:812885/1‑100 (MQ=255)
     cGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATc                                                                                            >  1:1452679/1‑100 (MQ=255)
      ggATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCa                                                                                           <  2:2989031/100‑1 (MQ=255)
      ggATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCa                                                                                           <  1:1918889/100‑1 (MQ=255)
       gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTGCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCAt                                                                                          <  2:3156405/100‑1 (MQ=255)
            ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCa                                                                                     >  1:3549605/1‑100 (MQ=255)
            ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCa                                                                                     >  1:1530234/1‑100 (MQ=255)
              aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                   <  2:2480118/100‑1 (MQ=255)
               aaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCg                                                                                  <  2:233007/100‑1 (MQ=255)
                   aCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGt                                                                              <  1:2715629/100‑1 (MQ=255)
                      gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGttt                                                                            >  1:589186/1‑99 (MQ=255)
                      gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTc                                                                           >  1:2119468/1‑100 (MQ=255)
                        cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTg                                                                         <  2:1227993/100‑1 (MQ=255)
                         aGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTgg                                                                        <  1:3199275/100‑1 (MQ=255)
                          gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTgg                                                                        >  2:421033/1‑99 (MQ=255)
                          gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGa                                                                       >  2:827711/1‑100 (MQ=255)
                          gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGa                                                                       >  2:136991/1‑100 (MQ=255)
                           tttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATAACGGTGGCTGAGCCAATCATAATCACCGGGGTTTCTGGAt                                                                      >  1:1244696/1‑100 (MQ=255)
                            ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTTCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                     >  2:2524290/1‑100 (MQ=255)
                            ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                     >  2:177065/1‑100 (MQ=255)
                            ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGAGCCTGAGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                     <  1:913229/100‑1 (MQ=255)
                             tttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtg                                                                    <  1:1111560/100‑1 (MQ=255)
                             tttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtg                                                                    <  2:961704/100‑1 (MQ=255)
                                aCCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggct                                                                 <  2:2958773/100‑1 (MQ=255)
                                 ccTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggctg                                                                <  2:1951425/100‑1 (MQ=255)
                                    ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                             <  2:1380785/100‑1 (MQ=255)
                                    ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                             <  2:2506218/100‑1 (MQ=255)
                                     gcgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCan                                                            <  2:1061289/100‑2 (MQ=255)
                                     gcgcTTAGTGCGCGGCGCGTTGCTGCGTATCGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCCCCTGTGTTTCTGGATTGGCTGGCAg                                                            <  2:2219623/100‑1 (MQ=255)
                                          tcgtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGa                                                       <  2:2380662/100‑1 (MQ=255)
                                            gtcgGCGGCGCGTTGCTGCATAAAGTCGCGGAAGGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                     <  2:3799281/100‑1 (MQ=255)
                                            gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                     >  2:124547/1‑100 (MQ=255)
                                            gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                     >  1:585892/1‑100 (MQ=255)
                                             tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                    <  1:377209/100‑1 (MQ=255)
                                             tcgGCGGCGCGGTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCTTGAGCCAATCATAATCGCCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                    <  1:3236986/100‑1 (MQ=255)
                                               ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCGTAATCACCGGGGTTTCTGGCTTGGCTGGCAGGCGGACGTTa                                                  >  1:1633511/1‑100 (MQ=255)
                                               ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTa                                                  <  1:1154398/100‑1 (MQ=255)
                                                gcggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTAt                                                 <  1:3402503/100‑1 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGCCGGAAGTTATCGt                                              <  1:2565741/100‑1 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                              >  1:477011/1‑100 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                              <  2:3643628/100‑1 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                              <  2:2756396/100‑1 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                              <  2:1550190/100‑1 (MQ=255)
                                                   gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                              <  1:3027462/100‑1 (MQ=255)
                                                    cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGtt                                             >  2:1450/1‑100 (MQ=255)
                                                     gcgTTGCTGCATAAAGGCGCGGAACGGGGCAATACCGGTGCCTGAGCCAATCATAAGCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTAGCGtt                                             <  2:589483/99‑1 (MQ=255)
                                                       gtcgctgcATAAAGGCGGTGAACGGCGGAATGCCGGCGCCGGAGCCAAGCATAATCTTCGGAGTTACTGGAGTGGCTGGCAGGGGGAAGTTATCGTTa                                            <  2:2250617/95‑1 (MQ=255)
                                                        ttGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGt                                         >  1:2453066/1‑100 (MQ=255)
                                                          gctgcATAAAGGCGCGGAACGGCGCACTACCGGTGCCTGAGCCAATCATAATCACCGGGGATTCTGGATAGGCTGGCAGGCGGAAGTTATCGTCATGTTc                                       >  2:679480/1‑100 (MQ=255)
                                                          gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                       >  2:840237/1‑100 (MQ=255)
                                                          gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                       >  1:2812452/1‑100 (MQ=255)
                                                          gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                       <  1:2720564/100‑1 (MQ=255)
                                                           ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCg                                      >  2:935829/1‑100 (MQ=255)
                                                           ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCg                                      >  1:2278362/1‑100 (MQ=255)
                                                            tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                     >  2:2887864/1‑100 (MQ=255)
                                                            tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                     >  2:796753/1‑100 (MQ=255)
                                                            tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                     >  2:807835/1‑100 (MQ=255)
                                                             gcATAAAAGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                    >  1:1999769/1‑100 (MQ=255)
                                                               aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATaa                                  >  1:2520595/1‑100 (MQ=255)
                                                                tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGGTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATaaa                                 >  1:1405497/1‑100 (MQ=255)
                                                                 aaaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAAt                                <  2:3772985/100‑1 (MQ=255)
                                                                 aaaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAAt                                <  1:2556586/100‑1 (MQ=255)
                                                                  aaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                               >  1:3584549/1‑100 (MQ=255)
                                                                     gcgcGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgc                            <  1:437607/100‑1 (MQ=255)
                                                                      cgcgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                           <  2:2151504/100‑1 (MQ=255)
                                                                       gcgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTGCTGGATTAGCTGGCAGGCGGAAGTTAGCGTTATGTTCGATAAATACGCgg                          >  2:2677163/1‑100 (MQ=255)
                                                                        cgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGa                         >  1:1838120/1‑100 (MQ=255)
                                                                         ggAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTGTCGTTATGTTCGATAAATACGCGGAc                        <  2:2878112/100‑1 (MQ=255)
                                                                         ggAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGAc                        <  2:1897775/100‑1 (MQ=255)
                                                                            aCGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTc                     <  1:1332864/100‑1 (MQ=255)
                                                                             cGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCg                    >  2:1624674/1‑100 (MQ=255)
                                                                             cGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCg                    <  1:808137/100‑1 (MQ=255)
                                                                              ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGTGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                   >  2:2031816/1‑100 (MQ=255)
                                                                              ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                   <  1:1493322/100‑1 (MQ=255)
                                                                              ggCGCAAGTGCGGTGCCTGGGGCGATCATAATCAGCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACtt                      <  2:3376082/97‑1 (MQ=255)
                                                                               gcgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTCATGTTCGATCAATACGCGGACTTCGcc                  >  1:3291335/1‑100 (MQ=255)
                                                                                cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                 <  1:2881760/100‑1 (MQ=255)
                                                                                cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                 >  1:2082395/1‑100 (MQ=255)
                                                                                  naaTACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGATCTGTTCGATAAATACGCGGACGTCGcc                  >  2:1545053/2‑97 (MQ=255)
                                                                                  cAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctc               <  1:1831197/100‑1 (MQ=255)
                                                                                   aaTACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctct              <  1:2646027/100‑1 (MQ=255)
                                                                                      aCCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttct           <  2:1838110/100‑1 (MQ=255)
                                                                                       ccGGTGCCTGGGCCACTCATAATCACCGGGGCTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttc            >  1:3743676/1‑98 (MQ=255)
                                                                                       ccGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTcttctt          >  2:2191440/1‑100 (MQ=255)
                                                                                       ccGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTACCGCTATGTTCGGTACATACGCGGACTTCCCCCtctt             >  2:3128959/1‑97 (MQ=255)
                                                                                        cGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTcttctt          <  2:3073689/99‑1 (MQ=255)
                                                                                         ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc        >  1:67485/1‑100 (MQ=255)
                                                                                         ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc        >  1:1208358/1‑100 (MQ=255)
                                                                                         ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc        >  2:348502/1‑100 (MQ=255)
                                                                                          gTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCa       <  1:3492453/100‑1 (MQ=255)
                                                                                          gTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCa       <  2:3399346/100‑1 (MQ=255)
                                                                                            gCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcg     >  1:3101173/1‑100 (MQ=255)
                                                                                              cTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgcg   >  1:1357020/1‑100 (MQ=255)
                                                                                               tGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCgg  <  2:3181429/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
AAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGG  >  NC_000913/2890424‑2890618

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑