Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,521 G→A 52.7% 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,5210GA52.7% 6.5 / 146.1 110P460S (CCA→TCA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (25/27);  new base A (28/30);  total (53/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 111 total aligned reads displayed.

AGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGT  >  NC_000913/2890425‑2890619
                                                                                                |                                                                                                  
aGTGCGGCTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGcc                                                                                                 >  2:1533490/1‑100 (MQ=255)
aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAATGGTGCGGAACGGCGCAATACCGGTGCCTGGGcc                                                                                                 <  1:1920645/100‑1 (MQ=255)
aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGcc                                                                                                 >  2:2210402/1‑100 (MQ=255)
aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGcc                                                                                                 >  2:471077/1‑100 (MQ=255)
 gTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCa                                                                                                >  2:858584/1‑100 (MQ=255)
     ggATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCa                                                                                            <  2:1939775/100‑1 (MQ=255)
     ggATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCa                                                                                            >  1:385892/1‑100 (MQ=255)
      gATTACCAAAGAACAGCCAGTTTTTACCCGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCAt                                                                                           >  2:2222940/1‑100 (MQ=255)
        ttACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCCATACCGGTGCCTGGGCCACACATaa                                                                                         >  1:103685/1‑100 (MQ=255)
        ttACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATaa                                                                                         <  2:523293/100‑1 (MQ=255)
          aCCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATc                                                                                       >  1:1902466/1‑100 (MQ=255)
          aCCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATc                                                                                       >  1:1212492/1‑100 (MQ=255)
           ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCa                                                                                      >  1:2895510/1‑100 (MQ=255)
           ccAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCa                                                                                      >  2:2369716/1‑100 (MQ=255)
            cAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCAc                                                                                     >  1:1337339/1‑100 (MQ=255)
             aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGGGGCATAAGGGCGCGGTACGTCGCAGTACCGGTGGCTGAGCCAATCATAATCAcc                                                                                    >  1:242122/1‑100 (MQ=255)
             aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                    <  1:50123/100‑1 (MQ=255)
               aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATCCCGGTGCCTGGGCCAATCATAATCACCgg                                                                                  >  1:1934361/1‑100 (MQ=255)
               aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcctg                                                                                  <  1:1995531/100‑3 (MQ=255)
               aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCgg                                                                                  >  2:1172973/1‑100 (MQ=255)
                   cAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGtt                                                                              >  1:2229767/1‑100 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTc                                                                            >  2:113737/1‑100 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCgg                                                                                  >  1:2375771/1‑94 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTc                                                                            >  2:2097968/1‑100 (MQ=255)
                     gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTc                                                                            >  1:2384066/1‑100 (MQ=255)
                      ccAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCt                                                                           >  1:1665535/1‑100 (MQ=255)
                      ccAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCt                                                                           <  1:1132289/100‑1 (MQ=255)
                       cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTg                                                                          >  2:563209/1‑100 (MQ=255)
                       cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTg                                                                          <  2:2064433/100‑1 (MQ=255)
                        aGTTTTTACCTGGCGCTTCGTCGGTGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTgg                                                                         <  2:2645215/100‑1 (MQ=255)
                         gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGa                                                                        <  2:1155977/100‑1 (MQ=255)
                         gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGa                                                                        <  2:2229781/100‑1 (MQ=255)
                          tttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGAt                                                                       <  2:2470355/100‑1 (MQ=255)
                          tttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCACGGAACGGCGCCATACCGGAACCAGGGCCAAACAAAATCACCAGGGTTTCTGGAt                                                                       >  2:2433652/1‑100 (MQ=255)
                           ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                      <  2:2459135/100‑1 (MQ=255)
                            tttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtg                                                                     <  2:343192/100‑1 (MQ=255)
                              tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggc                                                                   <  2:2132746/100‑1 (MQ=255)
                              tACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggc                                                                   <  2:2792455/100‑1 (MQ=255)
                               aCCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggct                                                                  <  2:1910825/100‑1 (MQ=255)
                               aCCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtggct                                                                  <  1:2332634/100‑1 (MQ=255)
                                 cTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctgg                                                                >  2:2726271/1‑100 (MQ=255)
                                  tGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAATGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATtggctggc                                                               <  1:1079555/100‑1 (MQ=255)
                                      gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGc                                                           >  2:2041121/1‑100 (MQ=255)
                                       cTTCGTCGGCGGAGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGACTGAGCCAATCATAATAACCGGGGTTTCTGGAATGGATGGCAGGCg                                                          <  1:1656991/100‑1 (MQ=255)
                                        ttCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCgg                                                         <  2:764258/100‑1 (MQ=255)
                                        ttCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCgg                                                         >  1:2784284/1‑100 (MQ=255)
                                          cgtccGCCGCGCGTTGATGCATAAACGCGCGGTACGGCGCCATACCGGTCCCTGAGCCCATCATAATCACCGGGGTTTCTGCATTGGCTGGCATGCgg                                                         <  1:1757351/98‑1 (MQ=255)
                                           gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGaa                                                       >  2:2322028/1‑99 (MQ=255)
                                           gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                      <  2:713018/100‑1 (MQ=255)
                                            tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                     <  2:654750/100‑1 (MQ=255)
                                             cggcggCGCGTTGCTGCATAAAGACGCAGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                     >  1:291057/1‑99 (MQ=255)
                                              ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTa                                                   >  2:2442217/1‑100 (MQ=255)
                                              ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTa                                                   >  2:1191441/1‑100 (MQ=255)
                                                 ggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCg                                                <  1:2369396/100‑1 (MQ=255)
                                                 ggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCg                                                <  1:114356/100‑1 (MQ=255)
                                                  gcgcGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                               <  2:1937053/100‑1 (MQ=255)
                                                   cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGtt                                              <  1:2505998/100‑1 (MQ=255)
                                                   cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGtt                                              <  2:495284/100‑1 (MQ=255)
                                                    gcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGt                                               <  1:1249366/98‑1 (MQ=255)
                                                    gcgTTGCGGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTa                                             <  1:1415341/100‑1 (MQ=255)
                                                       ttGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGt                                          <  1:2266329/100‑1 (MQ=255)
                                                        tgctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGtt                                         <  1:1927579/100‑1 (MQ=255)
                                                         gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                        <  2:1095565/100‑1 (MQ=255)
                                                          ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCg                                       <  2:2776774/100‑1 (MQ=255)
                                                          ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCg                                       <  1:1565822/100‑1 (MQ=255)
                                                           tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                      <  1:2274521/100‑1 (MQ=255)
                                                           tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                      >  2:1419497/1‑100 (MQ=255)
                                                           tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                      >  2:1401470/1‑100 (MQ=255)
                                                           tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                      >  1:882657/1‑100 (MQ=255)
                                                            gcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGAt                                     <  2:384499/100‑1 (MQ=255)
                                                               tAAAGGCGCGGAACGGCGCAATACCGGTGCCTAAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATaaa                                  >  2:2770024/1‑100 (MQ=255)
                                                                aaaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAAt                                 <  2:2017299/100‑1 (MQ=255)
                                                                 aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                                <  2:879444/100‑1 (MQ=255)
                                                                 aaGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                                >  1:2089391/1‑100 (MQ=255)
                                                                  aGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGCTATCGTTATGTTCGATACAtgc                               >  1:497264/1‑98 (MQ=255)
                                                                    gcgcGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgc                             <  2:2554305/100‑1 (MQ=255)
                                                                     cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                            >  1:1854418/1‑100 (MQ=255)
                                                                     cgcgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                            >  1:445274/1‑100 (MQ=255)
                                                                       cgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGCGGTTCCTGGATTGGCTGGCAGGCGGAAGATACCGTTATGTTCGATAAATACGCgg                           >  1:1712423/1‑99 (MQ=255)
                                                                       cgGAACGGCGCAATACCGGTGACTGGGCCAATCATGATCACAGGGGTTTCTGGATTGGCTGGCACGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                            >  1:1407741/1‑98 (MQ=255)
                                                                        ggAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGAc                         <  2:2296255/100‑1 (MQ=255)
                                                                        ggAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGATTCTGGATTGGCTGGCAGGCGGAAGTTATCGGTCTGTTCGATAAATACGCgg                           >  1:515868/1‑98 (MQ=255)
                                                                         gAACTGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGTGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                            <  1:202283/96‑1 (MQ=255)
                                                                          aaCGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACtt                       <  2:2375786/100‑1 (MQ=255)
                                                                             ggCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                    <  1:1500489/100‑1 (MQ=255)
                                                                             ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                    <  2:2614496/100‑1 (MQ=255)
                                                                              gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGcc                   >  1:2282828/1‑100 (MQ=255)
                                                                               cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                  <  2:2100608/100‑1 (MQ=255)
                                                                               cgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGccc                  <  1:646240/100‑1 (MQ=255)
                                                                                  aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctct               <  2:2187069/100‑1 (MQ=255)
                                                                                     aCCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttct            <  1:1925129/100‑1 (MQ=255)
                                                                                     aCCGGTGCCGGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttct            <  1:1172576/100‑1 (MQ=255)
                                                                                        ggTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc         >  2:2698626/1‑100 (MQ=255)
                                                                                        ggTGCCTGAGCCAGTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCCGAAGTTATCGCTCTGTTCAATAAATACGCGGCCTTCGCCCtcttct            >  1:1805085/1‑97 (MQ=255)
                                                                                        ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc         >  2:1785180/1‑100 (MQ=255)
                                                                                           gCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTCTGTTCGATAACTACCCGGACTTCGCCCTCTTCTTCCa        >  2:2478377/1‑98 (MQ=255)
                                                                                           gCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcg      >  1:1538630/1‑100 (MQ=255)
                                                                                             cTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgcg    >  2:1663976/1‑100 (MQ=255)
                                                                                              tGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCgg   <  1:1681301/100‑1 (MQ=255)
                                                                                               gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGt  <  2:2782508/100‑1 (MQ=255)
                                                                                                |                                                                                                  
AGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGT  >  NC_000913/2890425‑2890619

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: