Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,889,702 A→G 32.4% V133A (GTG→GCG)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,7020AG32.4% 99.0 / 82.1 96V133A (GTG→GCG) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base A (29/36);  new base G (16/15);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.80e-01

CTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTG  >  NC_000913/2889605‑2889800
                                                                                                 |                                                                                                  
cTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAcc                                                                                                  <  2:1114674/100‑1 (MQ=255)
    gACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGtt                                                                                              <  2:9171/100‑1 (MQ=255)
     aCTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGGAGCATCTGGTGCACCGGttt                                                                                             <  2:1365397/100‑1 (MQ=255)
        gatctgggTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc                                                                                          <  1:1127792/95‑1 (MQ=255)
        cGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAc                                                                                          >  2:2150241/1‑100 (MQ=255)
         gTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                         <  1:89635/100‑1 (MQ=255)
         gTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                         <  1:616559/100‑1 (MQ=255)
         gTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACg                                                                                         <  1:2071543/100‑1 (MQ=255)
             gggTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGAcagatc                                                                                     >  2:1127789/1‑95 (MQ=255)
             gggTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCt                                                                                     <  1:1097626/100‑1 (MQ=255)
              ggTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCtt                                                                                    >  1:876133/1‑100 (MQ=255)
               gTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCttt                                                                                   >  2:898951/1‑100 (MQ=255)
               gTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCttt                                                                                   >  1:1286358/1‑100 (MQ=255)
                ttCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCtttt                                                                                  <  1:955029/100‑1 (MQ=255)
                  cGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTc                                                                                <  1:959721/100‑1 (MQ=255)
                   gAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCa                                                                               <  2:1449372/100‑1 (MQ=255)
                    aGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAg                                                                              >  1:1888628/1‑100 (MQ=255)
                    aGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAg                                                                              >  1:1296948/1‑100 (MQ=255)
                    aGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAg                                                                              >  2:777113/1‑100 (MQ=255)
                      gTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGaa                                                                            >  2:366448/1‑100 (MQ=255)
                        gCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAAt                                                                           >  2:1272147/1‑99 (MQ=255)
                          agatTACGTTACTGTTCATTTCGTTAGTTGGCGGCAGCGCATCTAGACCTAGCGAGTGCAGCGGCTTGTGCACCGGTTTGACGTTCGTTTTCAGAATGc                                                                         <  2:93774/99‑1 (MQ=255)
                          agagTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCTAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGcc                                                                        >  2:1845690/1‑100 (MQ=255)
                           gagTACGTTACGGTTCTTGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCg                                                                       >  1:2205828/1‑100 (MQ=255)
                            agTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGt                                                                      <  2:2168485/100‑1 (MQ=255)
                            agTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGt                                                                      <  2:363676/100‑1 (MQ=255)
                             gTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTg                                                                     <  2:1878683/100‑1 (MQ=255)
                             gTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTg                                                                     <  2:1593603/100‑1 (MQ=255)
                             gTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTg                                                                     <  1:67401/100‑1 (MQ=255)
                               aCGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTgg                                                                    >  2:1671825/1‑99 (MQ=255)
                                 gTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAAc                                                                 <  1:1773695/100‑1 (MQ=255)
                                  ttACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACt                                                                >  2:849929/1‑100 (MQ=255)
                                  ttACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTAGAACt                                                                >  1:744596/1‑100 (MQ=255)
                                  ttACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGGGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACt                                                                <  1:471711/100‑1 (MQ=255)
                                   tACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTg                                                               <  1:2015195/100‑1 (MQ=255)
                                     cGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa                                                             >  2:587937/1‑100 (MQ=255)
                                     cGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa                                                             >  2:2189815/1‑100 (MQ=255)
                                     cGGTTCATGTCGTTAGCTGTCGCCAGCACATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGaa                                                             <  2:1182884/100‑1 (MQ=255)
                                       gTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAAc                                                           <  1:1876145/100‑1 (MQ=255)
                                       gTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAAc                                                           <  2:1676781/100‑1 (MQ=255)
                                          cATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTc                                                        >  1:1953491/1‑100 (MQ=255)
                                            tGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTg                                                      <  2:1601557/100‑1 (MQ=255)
                                             gTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTgg                                                     >  1:350685/1‑100 (MQ=255)
                                              tCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGc                                                    <  1:1833365/100‑1 (MQ=255)
                                                gTTAGCTGTCGCCAGTGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCgg                                                  <  1:2217321/100‑1 (MQ=255)
                                                  tAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCCGGCGGtt                                                <  1:1170545/100‑1 (MQ=255)
                                                   aGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCAGCGGGTTGACGTTCGTTTTCAGAATGCCGTGGAAATGAAACGTCTGGCGGtt                                                <  2:2023170/99‑1 (MQ=255)
                                                    gCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTgg                                              <  1:1447675/100‑1 (MQ=255)
                                                     cTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTAACGTTCTTTTTCAGGATGACGTGGAACTGAAACAGCTGTCGGTTg                                               >  2:2185481/1‑98 (MQ=255)
                                                      tGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTc                                            >  2:1630622/1‑100 (MQ=255)
                                                      tGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTgg                                              >  1:291593/1‑98 (MQ=255)
                                                       gTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCa                                           >  2:765988/1‑100 (MQ=255)
                                                       gTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCa                                           >  2:893644/1‑100 (MQ=255)
                                                         cGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAgg                                         >  2:2148417/1‑100 (MQ=255)
                                                         cGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAgg                                         >  1:2167386/1‑100 (MQ=255)
                                                           ccAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCg                                       >  1:545864/1‑100 (MQ=255)
                                                             aGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGaa                                     >  2:976651/1‑100 (MQ=255)
                                                               cgcATCAAGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtg                                   <  2:2227533/100‑1 (MQ=255)
                                                               cgcATCAAGACCGACCGAGTGCAGCAACTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtg                                   >  1:1365095/1‑100 (MQ=255)
                                                                gcATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgc                                  <  2:1716270/100‑1 (MQ=255)
                                                                 cATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgct                                 >  1:706240/1‑100 (MQ=255)
                                                                  aTCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAAtgctg                                >  2:424031/1‑100 (MQ=255)
                                                                    cAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGcc                              <  2:776261/100‑1 (MQ=255)
                                                                     aaGACCGACCGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCa                             >  2:1696529/1‑100 (MQ=255)
                                                                     aaGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCa                             <  2:1286352/100‑1 (MQ=255)
                                                                      aGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCAt                            <  1:976653/100‑1 (MQ=255)
                                                                      aGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCAt                            <  1:1693856/100‑1 (MQ=255)
                                                                       gaccgtccGAGTGAAGCGTTTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGATCTGAATTGTCTGGGGGTTGGTGAGGCGAATGCTGCCATa                           <  2:2018041/100‑1 (MQ=255)
                                                                       gaccgaccGAGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATa                           >  2:1288382/1‑100 (MQ=255)
                                                                       gaccgaccGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATa                           <  2:1604947/100‑1 (MQ=255)
                                                                         ccgaccgaGTGCAGCATCTGGTGCGCCGGTTTGACGATCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGa                         >  1:1578353/1‑100 (MQ=255)
                                                                         ccgaccgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGa                         <  2:876132/100‑1 (MQ=255)
                                                                            accgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATgg                      <  1:1592355/100‑1 (MQ=255)
                                                                              cgaGTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGgtg                    <  2:413450/100‑1 (MQ=255)
                                                                              cgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGgtg                    >  1:1028859/1‑100 (MQ=255)
                                                                              cgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGgtg                    <  1:1845702/100‑1 (MQ=255)
                                                                                 gTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGttt                 <  1:555418/100‑1 (MQ=255)
                                                                                 gTGCAGCATCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCCGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGtt                  >  2:1216156/1‑99 (MQ=255)
                                                                                 gTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGttt                 >  1:1118907/1‑100 (MQ=255)
                                                                                  tGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGtttt                >  1:860032/1‑100 (MQ=255)
                                                                                  tGCAGCATCTGGTGAACCGGTTTGACGTTATTTTTCAGAATGCAGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGttt                 <  2:659549/99‑1 (MQ=255)
                                                                                   gcagcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTGTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTc               >  2:698366/1‑100 (MQ=255)
                                                                                   gcagcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTc               >  2:2072255/1‑100 (MQ=255)
                                                                                      gcaTGTGGTGCACCGGTTTGACGTTCATTTTCATAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc            <  1:2179676/100‑1 (MQ=255)
                                                                                      gcaTCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc            <  1:1936973/100‑1 (MQ=255)
                                                                                      gcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCAcc            <  1:1613378/100‑1 (MQ=255)
                                                                                        aTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCgg          >  1:1402839/1‑100 (MQ=255)
                                                                                        aTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCgg          <  1:1030316/100‑1 (MQ=255)
                                                                                         tCTGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGc         >  2:893296/1‑100 (MQ=255)
                                                                                         tCTGGTGCGCCGGTTTGACGTTCGTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCgg          <  2:182443/99‑1 (MQ=255)
                                                                                         tCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGc         >  1:2024479/1‑100 (MQ=255)
                                                                                         tCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGc         <  2:2163269/100‑1 (MQ=255)
                                                                                           tGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCaa       >  1:1855973/1‑100 (MQ=255)
                                                                                           tGGTGCGCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCaa       >  2:91362/1‑100 (MQ=255)
                                                                                           tGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCaa       <  1:2150256/100‑1 (MQ=255)
                                                                                             gTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAAt     <  2:836751/100‑1 (MQ=255)
                                                                                                 aCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTg  <  1:777113/99‑1 (MQ=255)
                                                                                                 aCCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGCCTTGCGGTTGGTCAGGCGAATGCTGCCATAGATTGTTTTTTCACCGGCAAATTTg  <  2:1022366/99‑1 (MQ=255)
                                                                                                 |                                                                                                  
CTGCGACTCGTAAGGGTTCGAGGTGCAGAGTACGTTACGGTTCATGTCGTTAGCTGTCGCCAGCGCATCAAGACCGACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTG  >  NC_000913/2889605‑2889800

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: