Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 931,952 C→T V33V (GTG→GTA trxB ← thioredoxin reductase, FAD/NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913931,9520CT89.3% 227.4 / 16.7 84V33V (GTG→GTAtrxBthioredoxin reductase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (5/4);  new base T (34/41);  total (39/45)
Fisher's exact test for biased strand distribution p-value = 7.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TTAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCC  >  NC_000913/931855‑932048
                                                                                                 |                                                                                                
ttAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTAGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGtac                                                                                                <  2:1443047/100‑1 (MQ=255)
 taataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGtaca                                                                                               >  1:217557/1‑100 (MQ=255)
 taataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGtaca                                                                                               >  1:1294100/1‑100 (MQ=255)
  aataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcag                                                                                              >  2:3654/1‑100 (MQ=255)
   ataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcagg                                                                                             >  1:271004/1‑100 (MQ=255)
   ataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcagg                                                                                             >  2:1858949/1‑100 (MQ=255)
   ataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcagg                                                                                             <  1:1659367/100‑1 (MQ=255)
    taaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggt                                                                                            <  1:1145785/100‑1 (MQ=255)
     aaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggtt                                                                                           <  1:421109/100‑1 (MQ=255)
        ggACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgca                                                                                        <  2:1068793/100‑1 (MQ=255)
          aCCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcagg                                                                                      <  1:891551/100‑1 (MQ=255)
          aCCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcagg                                                                                      <  1:795389/100‑1 (MQ=255)
          aCCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcagg                                                                                      <  1:1094973/100‑1 (MQ=255)
            cGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggtt                                                                                    <  1:245895/100‑1 (MQ=255)
             ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggttg                                                                                   >  1:1882826/1‑100 (MQ=255)
             ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggtt                                                                                    >  1:776566/1‑99 (MQ=255)
              gTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTgg                                                                                  <  1:809657/100‑1 (MQ=255)
              gTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTgg                                                                                  <  2:1830424/100‑1 (MQ=255)
                  gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcg                                                                              >  1:1934228/1‑100 (MQ=255)
                  gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGACCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcg                                                                              >  2:1818739/1‑100 (MQ=255)
                  gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGGCAGTTGGCCGCTTTGTTCCATACCGGTAATCAGTACAGGTTGCATGTTGgcgcg                                                                              >  1:1362701/1‑100 (MQ=255)
                      gTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgcc                                                                          >  2:399510/1‑100 (MQ=255)
                      gTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgcc                                                                          >  2:105716/1‑100 (MQ=255)
                         tGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCcgcg                                                                       >  1:961147/1‑100 (MQ=255)
                         tGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCAAGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCcgcg                                                                       <  2:704834/100‑1 (MQ=255)
                          ggATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCcgcg                                                                       >  1:1239642/1‑99 (MQ=255)
                          ggATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGt                                                                      >  2:455664/1‑100 (MQ=255)
                           gATCGCCAGGCCAGTTTTCCCCTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                     >  1:1187507/1‑100 (MQ=255)
                              cGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAc                                                                  >  2:815311/1‑100 (MQ=255)
                              cGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAc                                                                  <  2:1657482/100‑1 (MQ=255)
                              cGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAc                                                                  >  2:670523/1‑100 (MQ=255)
                              cGCCAGGCCAGTTTTCCACGTCCGTTGTGGTGGTCAGTTTGTCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAc                                                                  <  2:647613/100‑1 (MQ=255)
                               gccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAca                                                                 >  2:1425547/1‑100 (MQ=255)
                               gccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAca                                                                 >  2:216192/1‑100 (MQ=255)
                                ccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcag                                                                <  2:144188/100‑1 (MQ=255)
                                ccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcag                                                                <  2:912288/100‑1 (MQ=255)
                                ccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcag                                                                >  1:1016039/1‑100 (MQ=255)
                                   ggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagcag                                                             >  2:751981/1‑100 (MQ=255)
                                   ggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACagcagc                                                            >  1:747162/1‑100 (MQ=255)
                                     ccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCg                                                           <  2:559299/100‑1 (MQ=255)
                                      cagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCgg                                                          <  1:1063878/100‑1 (MQ=255)
                                      cagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCgg                                                          <  2:1180078/100‑1 (MQ=255)
                                        gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                        <  2:946123/100‑1 (MQ=255)
                                        gTTTTCCACTTCCGTGGTGGTGGTCAGTTGACCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                        <  1:1818728/100‑1 (MQ=255)
                                         ttttCCACTTCCGTGGTGGTGTTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCATGTTGGCGCGCGCCGCGTAGACagcagc                                                            <  1:209987/95‑1 (MQ=255)
                                         ttttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtgt                                                       >  1:1862254/1‑100 (MQ=255)
                                            tCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATc                                                    >  1:179359/1‑100 (MQ=255)
                                              cACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTTTTTTTTTTTTTTTAGTACAGGGGGCAGGTTGGCGCGCGCGTGGTAGAGAGCTGCGGTGTAt                                                     >  1:1250633/1‑97 (MQ=255)
                                              cACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATccc                                                  >  1:1823414/1‑100 (MQ=255)
                                              cACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATccc                                                  <  2:999321/100‑1 (MQ=255)
                                                cTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgc                                                >  1:462108/1‑100 (MQ=255)
                                                cTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgc                                                >  2:906756/1‑100 (MQ=255)
                                                     gtggtggtgGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGc                                           >  2:456939/1‑100 (MQ=255)
                                                        gtggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTg                                        <  1:902689/100‑1 (MQ=255)
                                                        gtggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTg                                        <  1:744760/100‑1 (MQ=255)
                                                          ggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGaa                                      >  2:148528/1‑100 (MQ=255)
                                                            tggtCGGTTGGCCGCCTTTTTCCAGGCCGGTAATCGGTACAGGTTGCAGGTTGGCGCGTGTCGCGTAGACAGCGGCGGTGTATCCCGCCGGGCCTGAAcc                                    <  1:1066869/100‑1 (MQ=255)
                                                             ggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAccc                                   >  2:527365/1‑100 (MQ=255)
                                                               tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                 <  2:817932/100‑1 (MQ=255)
                                                               tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                 <  2:60644/100‑1 (MQ=255)
                                                               tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                 <  1:1876408/100‑1 (MQ=255)
                                                                 aGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGa                               >  2:1789191/1‑100 (MQ=255)
                                                                 aGTTGGCCGCCTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGAt                              >  2:120936/1‑100 (MQ=255)
                                                                     ggCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAg                           >  1:1370912/1‑100 (MQ=255)
                                                                       ccgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCa                         >  2:1665181/1‑100 (MQ=255)
                                                                       ccgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCa                         <  2:1199309/100‑1 (MQ=255)
                                                                       ccgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCa                         >  2:579982/1‑100 (MQ=255)
                                                                        cgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCTGCCGGGCCTGAACCCAGGATAAGCAg                        <  2:1187509/100‑1 (MQ=255)
                                                                           cTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGttt                     <  2:1254891/100‑1 (MQ=255)
                                                                            tttttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTa                    <  1:399508/100‑1 (MQ=255)
                                                                            tttttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTa                    >  1:1724986/1‑100 (MQ=255)
                                                                             ttttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTAc                   <  1:1858937/100‑1 (MQ=255)
                                                                             ttttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTAc                   <  2:1005138/100‑1 (MQ=255)
                                                                              tttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACt                  <  2:136485/100‑1 (MQ=255)
                                                                               ttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtg                 <  1:1207109/100‑1 (MQ=255)
                                                                               ttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtg                 >  1:1497765/1‑100 (MQ=255)
                                                                                 ccATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtgtg               <  2:1563712/100‑1 (MQ=255)
                                                                                 ccATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtgtg               <  1:815310/100‑1 (MQ=255)
                                                                                  cATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTgtgt              <  2:1823425/100‑1 (MQ=255)
                                                                                   aTGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGtt             <  2:1650049/100‑1 (MQ=255)
                                                                                   aTGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGtt             <  1:906755/100‑1 (MQ=255)
                                                                                       cGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCGGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTgg          <  1:811350/99‑1 (MQ=255)
                                                                                       cGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGt         <  2:1642406/100‑1 (MQ=255)
                                                                                         gTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCg       <  1:1410246/100‑1 (MQ=255)
                                                                                          tAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGt      >  1:719137/1‑100 (MQ=255)
                                                                                             tCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc   <  2:462110/100‑1 (MQ=255)
                                                                                              cAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGccc  >  1:833695/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TTAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCCC  >  NC_000913/931855‑932048

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