Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,521 G→A 68.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,5210GA68.7% 121.4 / 85.4 107P460S (CCA→TCA) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (14/19);  new base A (27/45);  total (42/65)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

Only 100 of 111 total aligned reads displayed.

TAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTC  >  NC_000913/2890422‑2890620
                                                                                                   |                                                                                                   
tAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTgg                                                                                                     >  2:522430/1‑100 (MQ=255)
  aaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGc                                                                                                   <  1:1101758/100‑1 (MQ=255)
  aaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGc                                                                                                   >  2:40599/1‑100 (MQ=255)
   aGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGcc                                                                                                  <  2:2229072/100‑1 (MQ=255)
    gTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCa                                                                                                 >  1:105821/1‑100 (MQ=255)
      gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCaa                                                                                                >  2:609984/1‑99 (MQ=255)
      gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAAt                                                                                               >  1:310215/1‑100 (MQ=255)
      gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAAt                                                                                               <  1:755477/100‑1 (MQ=255)
       cGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATc                                                                                              <  2:383901/100‑1 (MQ=255)
       cGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATc                                                                                              <  1:642705/100‑1 (MQ=255)
         gATTACCAAAGAACAGCCAGTTTTTTTTTTTTGTTTCATCCGCGGACCGTTGCCGCATAAAGGCGCGGAACGGCGCAATACCGGTGGCTGTGCCAATCa                                                                                             >  2:2101912/1‑99 (MQ=255)
         gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCan                                                                                            <  1:66382/100‑2 (MQ=255)
         gATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGACGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCAt                                                                                            >  2:1192148/1‑100 (MQ=255)
          aTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATa                                                                                           <  2:112123/100‑1 (MQ=255)
            tACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAAt                                                                                         <  2:1263686/100‑1 (MQ=255)
             aCCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGCCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATaa                                                                                          <  2:554987/98‑1 (MQ=255)
                aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                     >  1:2070936/1‑100 (MQ=255)
                aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                     >  2:141971/1‑100 (MQ=255)
                aaaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCAcc                                                                                     <  2:579269/100‑1 (MQ=255)
                 aaGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCg                                                                                    <  1:78875/100‑1 (MQ=255)
                  aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCgg                                                                                   <  2:549714/100‑1 (MQ=255)
                   gAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCggg                                                                                  >  2:861239/1‑100 (MQ=255)
                   gAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCggg                                                                                  <  2:360160/100‑1 (MQ=255)
                    aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCgggg                                                                                 >  1:441002/1‑100 (MQ=255)
                        gCCAGTTTTTACTTGGCTCTTCGTCGGCGGCGCGTTGCTGCATAAAGCCGCGGACCGGCGCAATACTGGTGCCTGGGCCAATCCTAATCACCGGGGtctc                                                                             >  2:504545/1‑100 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTc                                                                             >  2:2316893/1‑100 (MQ=255)
                        gCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTc                                                                             >  1:977268/1‑100 (MQ=255)
                          cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTg                                                                           <  1:284098/100‑1 (MQ=255)
                            gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGa                                                                         <  1:2288941/100‑1 (MQ=255)
                            gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGa                                                                         <  1:423464/100‑1 (MQ=255)
                              ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                       >  1:892508/1‑100 (MQ=255)
                              ttttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGAtt                                                                       <  2:1008855/100‑1 (MQ=255)
                               tttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtg                                                                      <  2:854009/100‑1 (MQ=255)
                                ttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTTAGCCAATCATAATCACCGGGGTTTCTGGAt                                                                        <  2:1068813/97‑1 (MQ=255)
                                ttACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATtgg                                                                     <  2:1381476/100‑1 (MQ=255)
                                    cTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAAAGGCGCAATACCGGTGCCTCAGCCAATCATAATCACCGGGGTTTCTGGATtggctgg                                                                 <  2:825681/100‑1 (MQ=255)
                                    cTGGCGATTCGCCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCTCCGGGGTTTCTGGATtggcttg                                                                 <  2:790473/100‑3 (MQ=255)
                                      ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                               >  1:1375313/1‑100 (MQ=255)
                                      ggCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                               >  1:434974/1‑100 (MQ=255)
                                      ggCGCTTCGTCGGCGACGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCa                                                               <  1:1192154/100‑1 (MQ=255)
                                       gcgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAg                                                              <  1:444545/100‑1 (MQ=255)
                                       gcgcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAg                                                              <  1:2013410/100‑1 (MQ=255)
                                         gcTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGc                                                            <  2:2226059/100‑1 (MQ=255)
                                          cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCg                                                           >  2:941824/1‑100 (MQ=255)
                                          cTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCg                                                           <  2:802672/100‑1 (MQ=255)
                                              gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                       <  1:938159/100‑1 (MQ=255)
                                              gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                       <  1:2139592/100‑1 (MQ=255)
                                              gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                       >  2:1246950/1‑100 (MQ=255)
                                              gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGCTTTCTGGATTGGCTGGCAGGCGGAAg                                                       <  1:1451213/100‑1 (MQ=255)
                                              gtcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAAGCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAg                                                       <  2:574525/100‑1 (MQ=255)
                                               tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                      <  1:332872/100‑1 (MQ=255)
                                               tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                      <  1:726851/100‑1 (MQ=255)
                                               tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGt                                                      <  1:1951003/100‑1 (MQ=255)
                                                 ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTa                                                    >  1:1038182/1‑100 (MQ=255)
                                                    ggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCg                                                 <  2:1485618/100‑1 (MQ=255)
                                                      cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGtt                                               <  1:1241262/100‑1 (MQ=255)
                                                        cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTAt                                             <  2:809733/100‑1 (MQ=255)
                                                        cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTAt                                             >  1:1115473/1‑100 (MQ=255)
                                                         gTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATg                                            <  2:498373/100‑1 (MQ=255)
                                                            gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGTGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                         <  2:501429/100‑1 (MQ=255)
                                                            gctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTc                                         >  2:621438/1‑100 (MQ=255)
                                                              tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                       <  1:2141646/100‑1 (MQ=255)
                                                              tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGa                                       >  1:2064241/1‑100 (MQ=255)
                                                                cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCAGAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                     <  1:1317790/100‑1 (MQ=255)
                                                                cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATa                                     <  1:1212061/100‑1 (MQ=255)
                                                                  tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATaaa                                   >  2:93575/1‑100 (MQ=255)
                                                                    aaGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATa                                 <  2:1468511/100‑1 (MQ=255)
                                                                      ggCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcg                               >  1:2089725/1‑100 (MQ=255)
                                                                      ggCGCGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcg                               >  2:1077929/1‑100 (MQ=255)
                                                                       gcgcGGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCGCCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgc                              <  2:1958503/100‑1 (MQ=255)
                                                                        cgcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                             >  1:323840/1‑100 (MQ=255)
                                                                        cgcgGAACGGCGCAATACCGGTGCCTGCGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATAcgcg                             <  1:1081961/100‑1 (MQ=255)
                                                                         gcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCgg                            <  1:861241/100‑1 (MQ=255)
                                                                          cgGAACGGCGCAATACCGGTGGCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGa                           <  1:2197175/100‑1 (MQ=255)
                                                                          cgGAACGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATagatcgg                           >  1:742247/1‑94 (MQ=255)
                                                                            nagcGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACt                         >  2:1144429/4‑100 (MQ=255)
                                                                             aaCGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACtt                        >  2:805246/1‑100 (MQ=255)
                                                                             aaCGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACtt                        <  1:1077934/100‑1 (MQ=255)
                                                                             aaCGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTGCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACtt                        >  1:1770590/1‑100 (MQ=255)
                                                                               cGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCg                      >  2:487778/1‑100 (MQ=255)
                                                                               cGGCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCg                      >  2:585921/1‑100 (MQ=255)
                                                                                ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                     <  2:1203148/100‑1 (MQ=255)
                                                                                ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                     >  1:1886109/1‑100 (MQ=255)
                                                                                ggCGCAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGc                     <  1:736791/100‑1 (MQ=255)
                                                                                 gcgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGcc                    >  1:1487795/1‑100 (MQ=255)
                                                                                 gcgcAATACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGcc                    <  1:1199564/100‑1 (MQ=255)
                                                                                    cAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctc                 >  1:1243432/1‑100 (MQ=255)
                                                                                     aaTACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGCCTTCGCCctct                >  2:1059912/1‑100 (MQ=255)
                                                                                     aaTACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCctct                >  2:2363661/1‑100 (MQ=255)
                                                                                      aTACCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtctt               <  2:2327091/100‑1 (MQ=255)
                                                                                       tACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttc              <  1:1680737/100‑1 (MQ=255)
                                                                                        aCCTGTGCCTTAGCCAATCATAATCACCTGTGTTTCTGGATTGGCTGGCAGGCTGAAGTTATCGTTATGTTCGATAAATACGCTGACTTCGCCCtcttct             <  2:428479/100‑1 (MQ=255)
                                                                                        aCCGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCtcttct             <  2:2328794/100‑1 (MQ=255)
                                                                                         ccGGTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTcttctt            >  2:416118/1‑100 (MQ=255)
                                                                                          cGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCttcttc           <  2:1188775/100‑1 (MQ=255)
                                                                                           ggTGCCTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTcc          >  2:1086536/1‑100 (MQ=255)
                                                                                                cTGAGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgcg     >  2:2082234/1‑100 (MQ=255)
                                                                                                  gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGt   >  1:2167878/1‑100 (MQ=255)
                                                                                                   ggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTc  <  1:14713/100‑1 (MQ=255)
                                                                                                   ggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTc  <  2:451258/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTC  >  NC_000913/2890422‑2890620

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: