Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 931,952 C→T 100% 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,9520CT97.0% 313.9 / ‑0.2 98V33V (GTG→GTAtrxBthioredoxin reductase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (46/49);  total (47/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCC  >  NC_000913/931856‑932047
                                                                                                |                                                                                               
taataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGtaca                                                                                              >  2:372806/1‑100 (MQ=255)
 aataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcag                                                                                             >  2:391089/1‑100 (MQ=255)
 aataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcag                                                                                             >  1:1475646/1‑100 (MQ=255)
 aataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcag                                                                                             <  2:1977189/100‑1 (MQ=255)
 aataaCGGACCGGCCAGATCGTTTGGAACGCCAGGCCAGTATCACACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGtaca                                                                                              >  2:1049855/1‑99 (MQ=255)
  ataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcagg                                                                                            >  1:1512934/1‑100 (MQ=255)
   taaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggt                                                                                           <  2:952378/100‑1 (MQ=255)
   taaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggt                                                                                           >  2:2556470/1‑100 (MQ=255)
      cGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgc                                                                                        >  2:208161/1‑100 (MQ=255)
      cGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgc                                                                                        >  1:1519269/1‑100 (MQ=255)
       ggACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgca                                                                                       <  2:1667602/100‑1 (MQ=255)
       ggACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGCGGTGGTGGTCAGTTGGCCGCCTTTTTCCCTGCCGGTAATCAGTAcaggttgca                                                                                       >  2:894925/1‑100 (MQ=255)
         aCCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcagg                                                                                     <  1:25755/100‑1 (MQ=255)
          ctggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggt                                                                                    <  1:1231321/98‑1 (MQ=255)
          ccGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggt                                                                                    >  1:1752117/1‑100 (MQ=255)
          ccGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggt                                                                                    >  1:2173751/1‑100 (MQ=255)
            ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTGCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggtt                                                                                   >  2:1573305/1‑99 (MQ=255)
            ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggttg                                                                                  >  2:2255779/1‑100 (MQ=255)
            ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggttg                                                                                  >  1:425426/1‑100 (MQ=255)
            ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggtag                                                                                  >  2:1231313/1‑98 (MQ=255)
             gTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTgg                                                                                 <  1:1583623/100‑1 (MQ=255)
             gTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTgg                                                                                 <  1:268376/100‑1 (MQ=255)
               cAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTACATGCCGGTAATCAGTACAGGTCGCAGGTTGgcg                                                                               >  1:211364/1‑100 (MQ=255)
                aGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgc                                                                              <  2:1545937/100‑1 (MQ=255)
                 gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcg                                                                             <  2:2479415/100‑1 (MQ=255)
                 gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcg                                                                             >  1:1771379/1‑100 (MQ=255)
                 gATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcg                                                                             >  1:880570/1‑100 (MQ=255)
                   tCGTTTGGATCGCCAGGCCAGTTTTCTACTTCCGGGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcgcg                                                                           <  2:1860472/100‑1 (MQ=255)
                   tCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcgcg                                                                           >  1:569294/1‑100 (MQ=255)
                   tCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGgcgcgcg                                                                           <  1:768074/100‑1 (MQ=255)
                    cGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcgcgc                                                                          <  2:2024723/100‑1 (MQ=255)
                     gTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgcc                                                                         >  2:1916110/1‑100 (MQ=255)
                       ttGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGGCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgccgc                                                                       <  2:1373032/100‑1 (MQ=255)
                        tGGATCGCCGGGCCAGTTTTACACTTCCGTGGTAGTGGTCAGTTGGCCGCCTTTTTCCACGACGGTAATCAGTACAGGCTGCAGGTTGGCGCGCGCcgcg                                                                      <  2:1260052/100‑1 (MQ=255)
                        tGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCcgcg                                                                      <  1:1117410/100‑1 (MQ=255)
                        tGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGACGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCcgcg                                                                      <  1:1019177/100‑1 (MQ=255)
                         ggATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGt                                                                     <  1:2556486/100‑1 (MQ=255)
                          gATCGCCAGGGCAGTTTTCCGCTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCAGGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                    <  2:726959/100‑1 (MQ=255)
                          gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                    >  1:1167390/1‑100 (MQ=255)
                          gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                    <  2:1593798/100‑1 (MQ=255)
                          gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                    <  2:1033954/100‑1 (MQ=255)
                           aTCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAg                                                                   <  1:363526/100‑1 (MQ=255)
                            tCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGa                                                                  >  1:78065/1‑100 (MQ=255)
                               ccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcag                                                               <  2:493800/100‑1 (MQ=255)
                               ccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcag                                                               >  1:86346/1‑100 (MQ=255)
                                caggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagc                                                              <  1:1204378/100‑1 (MQ=255)
                                caggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATAAGTACAGGTTGCAGGTTGGCGCGCTCCGCGTAGAcagc                                                              <  2:1585166/100‑1 (MQ=255)
                                 tctccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagca                                                             <  1:1761332/97‑1 (MQ=255)
                                  ggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACagcagc                                                           >  1:2386294/1‑100 (MQ=255)
                                   gccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACagcagc                                                           >  1:1702450/1‑100 (MQ=255)
                                    ccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagcaaga                                                          >  2:1761319/1‑97 (MQ=255)
                                    ccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCg                                                          >  2:1674436/1‑100 (MQ=255)
                                      atttttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCg                                                          <  2:567954/96‑1 (MQ=255)
                                      agTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgt                                                        >  1:2024266/1‑100 (MQ=255)
                                       gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTTTTTTTTTTTATTCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       >  2:2247680/1‑100 (MQ=255)
                                       gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  1:2437692/100‑1 (MQ=255)
                                       gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  1:18108/100‑1 (MQ=255)
                                       gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  2:809256/100‑1 (MQ=255)
                                       gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  1:2520409/100‑1 (MQ=255)
                                          ttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTAt                                                    >  2:447375/1‑100 (MQ=255)
                                           tCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATc                                                   <  2:982529/100‑1 (MQ=255)
                                             cACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATccc                                                 >  1:2014098/1‑100 (MQ=255)
                                               cTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgc                                               >  1:855698/1‑100 (MQ=255)
                                                ttCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgcc                                              >  1:2292979/1‑100 (MQ=255)
                                                 tCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgccg                                             >  1:2500028/1‑100 (MQ=255)
                                                  ccGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCgg                                            <  2:969670/100‑1 (MQ=255)
                                                   cGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCggg                                           <  2:1521675/100‑1 (MQ=255)
                                                   cGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCggg                                           <  1:223144/100‑1 (MQ=255)
                                                   cGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGAGTAGACAGCAGCGGTGTATCCCGCCggg                                           <  1:2503302/100‑1 (MQ=255)
                                                    gtggtggtgGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGc                                          <  2:2497620/100‑1 (MQ=255)
                                                     tggtggtggTCAGTTGGCTGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCTCGCCGGGcc                                         <  2:850718/100‑1 (MQ=255)
                                                        tggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGa                                      >  1:830946/1‑100 (MQ=255)
                                                         ggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGaa                                     >  2:1459805/1‑100 (MQ=255)
                                                         ggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGaa                                     >  1:1371477/1‑100 (MQ=255)
                                                         ggtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGaa                                     >  1:1169254/1‑100 (MQ=255)
                                                          gtggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAc                                    <  1:258443/100‑1 (MQ=255)
                                                           tggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAcc                                   <  1:2196077/100‑1 (MQ=255)
                                                             gtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCa                                 >  2:844615/1‑100 (MQ=255)
                                                              tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                >  1:1176607/1‑100 (MQ=255)
                                                              tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                <  2:2024252/100‑1 (MQ=255)
                                                              tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                <  1:1433676/100‑1 (MQ=255)
                                                                aGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGa                              >  1:1260770/1‑100 (MQ=255)
                                                                 gTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGAt                             >  1:1321470/1‑100 (MQ=255)
                                                                  ttGGCTGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATa                            <  1:890748/100‑1 (MQ=255)
                                                                  ttGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATa                            <  2:1990048/100‑1 (MQ=255)
                                                                     gccgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCAGTCTATCCCGCCGGGCCTGAACCCAGGATAAGc                         <  2:1679939/100‑1 (MQ=255)
                                                                        gccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGt                      <  2:222139/100‑1 (MQ=255)
                                                                           tttttCGCTGCTGGTAATCAGTACGGGTTGCAGGTTGGCGCGTGCCGCGTACACCGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGttt                    <  2:211362/99‑1 (MQ=255)
                                                                                ccATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtgtg              >  2:252977/1‑100 (MQ=255)
                                                                                 cATGCCGCTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTgtgt             >  2:428810/1‑100 (MQ=255)
                                                                                   tGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGttt           <  1:509839/100‑1 (MQ=255)
                                                                                    gCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTg          >  1:76509/1‑100 (MQ=255)
                                                                                    gCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTg          <  2:1300454/100‑1 (MQ=255)
                                                                                    gCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTg          >  1:2112117/1‑100 (MQ=255)
                                                                                      cGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGt        >  1:1524005/1‑100 (MQ=255)
                                                                                       ggTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTc       >  1:245241/1‑100 (MQ=255)
                                                                                            tCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc  <  2:76508/100‑1 (MQ=255)
                                                                                            tCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc  <  2:1089022/100‑1 (MQ=255)
                                                                                            tCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc  >  2:847425/1‑100 (MQ=255)
                                                                                            tCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc  <  1:2351453/100‑1 (MQ=255)
                                                                                                |                                                                                               
TAATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCC  >  NC_000913/931856‑932047

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: