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,9520CT95.1% 292.7 / 3.6 100V33V (GTG→GTAtrxBthioredoxin reductase, FAD/NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base T (46/49);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

AATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCC  >  NC_000913/931857‑932047
                                                                                               |                                                                                               
aataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcag                                                                                             <  2:905456/100‑1 (MQ=255)
 ataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcagg                                                                                            >  1:1570792/1‑100 (MQ=255)
 ataaCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcagg                                                                                            >  2:1963354/1‑100 (MQ=255)
    aCGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttg                                                                                         >  2:3010269/1‑100 (MQ=255)
     cGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgc                                                                                        <  1:2112212/100‑1 (MQ=255)
       gACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgca                                                                                       >  2:2556514/1‑99 (MQ=255)
         ccGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggt                                                                                    <  2:2317438/100‑1 (MQ=255)
         ccGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttgcaggt                                                                                    >  2:596600/1‑100 (MQ=255)
         ccGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCAcaggttgcaggt                                                                                    >  1:2245969/1‑100 (MQ=255)
         ccGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGACTTTTTCCATGCCGGTAATCAGTAcaggttgcaggt                                                                                    >  1:1003831/1‑100 (MQ=255)
           ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTTGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggttg                                                                                  >  1:755680/1‑100 (MQ=255)
           ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTAcaggttgcaggttg                                                                                  >  2:1309595/1‑100 (MQ=255)
           ggTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTACATGCCGGTAATCAGTAcaggttgaaggttg                                                                                  >  1:1264087/1‑100 (MQ=255)
            gTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTgg                                                                                 <  2:551333/100‑1 (MQ=255)
            gTCAGATCGTTTGGATAGAAACGCCAGTTTTCCACTTCAGTTGTGGTGATCAATTGTCCTCCTTTTTCCATTCCTTTAATCGGTACAGGTTGTAGGTTgg                                                                                 >  2:1763257/1‑100 (MQ=255)
             tCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTACGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGAAATCAGTACAGGTTGCAGGTTGgc                                                                                >  1:3022835/1‑100 (MQ=255)
                 aTCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGgcgcgc                                                                            >  1:1763344/1‑100 (MQ=255)
                    gTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgcc                                                                         >  1:881618/1‑100 (MQ=255)
                    gTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgcc                                                                         <  2:2298541/100‑1 (MQ=255)
                    gTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgcc                                                                         >  2:1096345/1‑100 (MQ=255)
                    gTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGcgcc                                                                         >  1:1461653/1‑100 (MQ=255)
                        ggATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGt                                                                     <  2:1554299/100‑1 (MQ=255)
                         gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTTGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGtn                                                                    <  1:164674/100‑2 (MQ=255)
                         gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                    <  1:2317604/100‑1 (MQ=255)
                         gATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTa                                                                    <  1:2473681/100‑1 (MQ=255)
                          aTCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCTTTTTTCCATGCCGGTAATCATTACAGGTTTCAGGTTGGCGTGCGGCGCTTAg                                                                   <  2:1763324/100‑1 (MQ=255)
                          aTCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGTTAATCAGTACAGGTTGCAGTTTGGCGCGCGCCGCGTAg                                                                   >  1:399575/1‑100 (MQ=255)
                           tCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGa                                                                  <  1:1247385/100‑1 (MQ=255)
                              ccaggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGACGCGCGCCGCGTAGAcag                                                               >  1:2723235/1‑100 (MQ=255)
                               caggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagc                                                              <  1:2253481/100‑1 (MQ=255)
                               caggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagc                                                              <  1:179270/100‑1 (MQ=255)
                                 ggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagcag                                                            >  2:1151036/1‑100 (MQ=255)
                                 ggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagcag                                                            >  2:1729065/1‑100 (MQ=255)
                                 ggccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGAcagcag                                                            >  1:2636912/1‑100 (MQ=255)
                                  gccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTTTTTTTTTTTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACagcagc                                                           >  2:2389046/1‑100 (MQ=255)
                                  gccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACagcagc                                                           <  2:2037910/100‑1 (MQ=255)
                                  gccagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACagcagc                                                           <  2:725919/100‑1 (MQ=255)
                                    cagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCgg                                                         <  1:92950/100‑1 (MQ=255)
                                    cagTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCgg                                                         <  2:1819370/100‑1 (MQ=255)
                                      gTTTTCCACTTCCTTGGTGGTGTTTATTTGGCCGCATTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  1:1792071/100‑1 (MQ=255)
                                      gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTTCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  1:811253/100‑1 (MQ=255)
                                      gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATTCCGGTAATCAGTACAGGTTGCAAGTTAGAGCGCGCCGCGTAGACAGCAGCGgtg                                                       >  1:2507146/1‑100 (MQ=255)
                                      gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  1:2211772/100‑1 (MQ=255)
                                      gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  1:1377931/100‑1 (MQ=255)
                                      gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       >  2:2908675/1‑100 (MQ=255)
                                      gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  2:2628305/100‑1 (MQ=255)
                                      gTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGgtg                                                       <  2:2697299/100‑1 (MQ=255)
                                        tttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTa                                                     >  2:918094/1‑100 (MQ=255)
                                         ttCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTTTATTTTTTTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTAt                                                    >  2:697258/1‑100 (MQ=255)
                                         ttCCACTACCGCGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTAt                                                    >  2:1582251/1‑100 (MQ=255)
                                           ccACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCTATGTTGTTCATCAGTACAGGTTGCAGGTTGGCGCGCGCCTCGTAGACAGCAGCGGTGTa                                                     >  2:1478319/1‑97 (MQ=255)
                                           ccACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTACTCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATcc                                                  >  1:1087747/1‑100 (MQ=255)
                                            cACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATccc                                                 >  2:1691323/1‑100 (MQ=255)
                                             aCTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccg                                                >  1:3006661/1‑100 (MQ=255)
                                              cTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATccc                                                 <  2:755676/98‑1 (MQ=255)
                                               ttCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgcc                                              <  2:900065/100‑1 (MQ=255)
                                                tccggtggggtggtCAGTTGGCCGCCGGTTTCCATGCCGGTAATCAGTACATGTTGCAGGTTGGCTCGCGCCTCGTAGACAGCAGCGGTGTACAccgc                                               <  2:2413384/90‑1 (MQ=255)
                                                tCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGACGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCccgccg                                             <  1:2111832/100‑1 (MQ=255)
                                                 ccGTGGTGGTGGTCAGTTGGTCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCgg                                            <  1:2913668/100‑1 (MQ=255)
                                                  cGTGGTGTTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCggg                                           <  1:1317107/100‑1 (MQ=255)
                                                      gtggtggtCAGTTGGCCGCCTTTTTCCATGCCGGAAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTg                                       <  2:3022802/100‑1 (MQ=255)
                                                          tggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAc                                    >  2:484283/1‑99 (MQ=255)
                                                           ggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAAccc                                  >  2:1605970/1‑100 (MQ=255)
                                                           ggtCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGATGCTGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCAGGGCCTGAAccc                                  >  2:864648/1‑100 (MQ=255)
                                                             tCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAg                                <  1:1729085/100‑1 (MQ=255)
                                                              cAGTTGGCCGCCTTTTTCCATGCCTGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAgg                               <  1:640349/100‑1 (MQ=255)
                                                              cAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAgg                               <  2:399572/100‑1 (MQ=255)
                                                              cAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAgg                               >  2:1918598/1‑100 (MQ=255)
                                                              cAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAgg                               <  2:2245946/100‑1 (MQ=255)
                                                               aGTTGGCCGCCTTTTTCCTTGCCGGTAATCAGTACAGGTTGCAGGTTGGCGTGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGa                              <  1:204755/100‑1 (MQ=255)
                                                                  tGGCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGATGTATCCCGCCGGGCCTGAACCCAGGATaa                           >  1:1539971/1‑100 (MQ=255)
                                                                   ggCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAg                          <  1:353341/100‑1 (MQ=255)
                                                                   ggCCGCCTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAg                          >  2:2574044/1‑100 (MQ=255)
                                                                   ggCCGCCTTTTTACATACCGGTAATCAGTACAGGTTGCAGGTTAGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAg                          >  2:1819223/1‑100 (MQ=255)
                                                                    gccgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCACGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGc                         >  2:1094224/1‑100 (MQ=255)
                                                                    gccgccTTTTTCCATGCCGGTAATCAGTACAGGGTGCAGGTTGGCGCGCGCCGCGTAGAAAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAg                          >  2:1564113/1‑99 (MQ=255)
                                                                      cgccTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAg                       <  1:1720263/100‑1 (MQ=255)
                                                                      cgccTTTTTCCATGCCGGTAAACAGTACAGGTTTGAGGTTGGCGCGCGCCGCGTGGAGAGCAGCGGTGTATCCCGACGGGCCTGAACTGAGGATaagaag                       >  2:600208/1‑100 (MQ=255)
                                                                        ccTTTTTCCATGCCGGTAATCAGTACAGGTTGCTGGTTGGCGCGCGCTGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGtt                     <  1:2307900/100‑1 (MQ=255)
                                                                         cTTTTTCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGc                         <  2:147621/95‑1 (MQ=255)
                                                                         cTTTTTCCATGCCGGTAATCAGTACAGATTGCAGGTTGGCGCGCACCGCGTAGACAGCAGCGGTGTATCCCGCCGAGCCTGAACTCAGGATAAGCAGttt                    >  2:1628274/1‑100 (MQ=255)
                                                                           ttttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTAc                  <  1:2908707/100‑1 (MQ=255)
                                                                            tttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACt                 <  1:280549/100‑1 (MQ=255)
                                                                             ttCCATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACtg                <  2:1052346/100‑1 (MQ=255)
                                                                                cctGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGACGGGCCTGAACCCAGGATAAGCAGTTTACTgtgt             >  2:334143/3‑100 (MQ=255)
                                                                                cATGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTgtgt             >  2:2168276/1‑100 (MQ=255)
                                                                                 aTGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCTCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGtt            <  1:1811834/100‑1 (MQ=255)
                                                                                  tGCCGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGttt           <  2:1464082/100‑1 (MQ=255)
                                                                                    ccGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTgg         <  1:1414054/100‑1 (MQ=255)
                                                                                    ccGGTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTg          <  2:2474213/99‑1 (MQ=255)
                                                                                      ggTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTc       >  2:2556712/1‑100 (MQ=255)
                                                                                      ggTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTc       <  1:402007/100‑1 (MQ=255)
                                                                                       gTAATCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCg      >  2:510373/1‑100 (MQ=255)
                                                                                         aaTCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGACTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTg    >  1:2759952/1‑100 (MQ=255)
                                                                                          aTCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGc   >  1:2231081/1‑100 (MQ=255)
                                                                                          aTCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGc   <  2:1073471/100‑1 (MQ=255)
                                                                                          aTCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGc   <  2:1095491/100‑1 (MQ=255)
                                                                                           tCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGTGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc  <  1:1272670/100‑1 (MQ=255)
                                                                                           tCAGTACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc  >  1:877525/1‑100 (MQ=255)
                                                                                           tCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGcc  <  1:596604/100‑1 (MQ=255)
                                                                                               |                                                                                               
AATAACGGACCGGTCAGATCGTTTGGATCGCCAGGCCAGTTTTCCACTTCCGTGGTGGTGGTCAGTTGGCCGCCTTTTTCCATGCCGGTAATCAGCACAGGTTGCAGGTTGGCGCGCGCCGCGTAGACAGCAGCGGTGTATCCCGCCGGGCCTGAACCCAGGATAAGCAGTTTACTGTGTTTGGTCGTGCC  >  NC_000913/931857‑932047

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