Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,032,708 C→G 5.5% H448D (CAT→GAT)  uacT → uric acid permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,032,7080CG 5.5% 391.1 / 9.3 108H448D (CAT→GAT) uacTuric acid permease
Reads supporting (aligned to +/- strand):  ref base C (51/51);  new base G (3/3);  total (54/54)
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.99e-01

Only 100 of 111 total aligned reads displayed.

GTCTCGGCGTTGGCATGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTA  >  NC_000913/3032562‑3032828
                                                                                                                                                  |                                                                                                                        
  ctcGGCGTTGGCATGACTCCGACGCTCTCACACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGGACGTCTTCTTTAATGGCTATCAGCATCATg                                                                                                                       >  1:897762/1‑148 (MQ=255)
    cGGCGTTGGCATGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCt                                                                                                                     >  3:314161/1‑148 (MQ=255)
    cGGCGTTGGCATGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCt                                                                                                                     >  4:202375/1‑148 (MQ=255)
     ggCGTTGGCATGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTg                                                                                                                    <  3:202375/148‑1 (MQ=255)
        gTTGGCATGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTGGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGAcc                                                                                                                 <  2:897762/148‑1 (MQ=255)
             cATGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTg                                                                                                            >  3:569035/1‑148 (MQ=255)
               tGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGCCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGaa                                                                                                          >  3:794350/1‑148 (MQ=255)
                    ccGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGa                                                                                                           >  3:336078/1‑142 (MQ=255)
                    ccGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGa                                                                                                           <  4:336078/142‑1 (MQ=255)
                    ccGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAAt                                                                                                     >  3:212149/1‑148 (MQ=255)
                    ccGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAAt                                                                                                     >  2:76764/1‑148 (MQ=255)
                       aCGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCg                                                                                                  <  1:76764/148‑1 (MQ=255)
                        cGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGt                                                                                                 <  4:794350/148‑1 (MQ=255)
                         gCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGtc                                                                                                <  4:569035/148‑1 (MQ=255)
                            ctctcACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATGAGCATCATGCTGACCTGGTGa                                                                                                           >  4:741516/1‑134 (MQ=255)
                            ctctcACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATGATGCTGatagc                                                                                                               <  2:583420/130‑5 (MQ=255)
                            ctctcACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATGATGCTGatagc                                                                                                               >  1:583420/1‑126 (MQ=255)
                            ctctcACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAAAGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGa                                                                                                           <  3:741516/134‑1 (MQ=255)
                                cACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTg                                                                                                               <  1:240636/126‑1 (MQ=255)
                                cACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTg                                                                                                               >  2:240636/1‑126 (MQ=255)
                                   gATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAg                                                                                      <  4:212149/148‑1 (MQ=255)
                                     tttcttttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGaa                                                                                                          >  4:665823/1‑126 (MQ=255)
                                     tttcttttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGaa                                                                                                          <  3:665823/126‑1 (MQ=255)
                                       tcttttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGtc                                                                                                <  4:623955/134‑1 (MQ=255)
                                       tcttttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGtc                                                                                                >  3:623955/1‑134 (MQ=255)
                                       tcttttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGAttt                                                                                  <  4:314161/148‑1 (MQ=255)
                                          tttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGAcc                                                                                                                 >  1:456923/1‑114 (MQ=255)
                                          tttctAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGAcc                                                                                                                 <  2:456923/114‑1 (MQ=255)
                                               aaGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCt                                                                                                                <  1:741545/110‑1 (MQ=255)
                                               aaGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCt                                                                                                                >  2:741545/1‑110 (MQ=255)
                                                      cGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTgg                                                                                                              <  4:428248/105‑1 (MQ=255)
                                                      cGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTgg                                                                                                              >  3:428248/1‑105 (MQ=255)
                                                        gCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGa                                                                                                           >  4:325008/1‑106 (MQ=255)
                                                        gCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGa                                                                                                           <  3:325008/106‑1 (MQ=255)
                                                              cTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGtc                                                                                                <  2:687047/111‑1 (MQ=255)
                                                                   aCCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAgg                                                                                                        >  3:523731/1‑98 (MQ=255)
                                                                   aCCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAgg                                                                                                        <  4:523731/98‑1 (MQ=255)
                                                                   aCCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGGGAAGGAATCCGTCTCTGATaa                                                                                        >  1:233194/1‑114 (MQ=255)
                                                                    ccGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTc                                                     >  4:1117466/1‑148 (MQ=255)
                                                                        tgctACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTaaa                                                                               <  2:549345/118‑1 (MQ=255)
                                                                        tgctACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTaaa                                                                               >  1:549345/1‑118 (MQ=255)
                                                                         gctACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGtct                                                                                               >  1:553368/1‑101 (MQ=255)
                                                                           tACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTgg                                                                                                              <  3:410440/84‑1 (MQ=255)
                                                                           tACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTgg                                                                                                              >  4:410440/1‑84 (MQ=255)
                                                                             cATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATGATGCTGata                                                                                                                 <  1:761139/79‑3 (MQ=255)
                                                                             cATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATGATGCTGata                                                                                                                 >  2:761139/1‑77 (MQ=255)
                                                                             cATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAAt                                                                                                     <  1:982746/91‑1 (MQ=255)
                                                                             cATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAAt                                                                                                     >  2:982746/1‑91 (MQ=255)
                                                                               tAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGt                                                                   >  4:1105337/1‑123 (MQ=255)
                                                                                aGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTc                                                                           <  2:405677/114‑1 (MQ=255)
                                                                                aGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTc                                                                           >  1:405677/1‑114 (MQ=255)
                                                                                  cGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTaaa                                                                               >  2:1079561/1‑108 (MQ=255)
                                                                                  cGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTaaa                                                                               <  1:1079561/108‑1 (MQ=255)
                                                                                  cGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTa                                                              <  1:951270/125‑1 (MQ=255)
                                                                                  cGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTa                                                              >  2:951270/1‑125 (MQ=255)
                                                                                     cATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAg                                                                             <  1:832197/107‑1 (MQ=255)
                                                                                     cATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAg                                                                             >  2:832197/1‑107 (MQ=255)
                                                                                     cATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTAtgt                                                           <  1:572895/125‑1 (MQ=255)
                                                                                     cATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTAtgt                                                           >  2:572895/1‑125 (MQ=255)
                                                                                      aTTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATgcg                                              <  1:980813/137‑1 (MQ=255)
                                                                                      aTTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATgcg                                              >  2:980813/1‑137 (MQ=255)
                                                                                       ttATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCt                                                                                                                <  1:102873/70‑1 (MQ=255)
                                                                                       ttATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCt                                                                                                                >  2:102873/1‑70 (MQ=255)
                                                                                       ttATGCTCACAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTc                                                                           <  1:171197/107‑1 (MQ=255)
                                                                                       ttATGCTCACAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTc                                                                           >  2:171197/1‑107 (MQ=255)
                                                                                            cTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATg                             <  3:1117466/148‑1 (MQ=255)
                                                                                               gCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGAcc                                                                                                                 <  1:683616/61‑1 (MQ=255)
                                                                                               gCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGAcc                                                                                                                 >  2:683616/1‑61 (MQ=255)
                                                                                                             cGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAg                                          <  1:77947/118‑1 (MQ=255)
                                                                                                             cGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTAGCTTCTGATGCGCAAGCTGAAGaaaaa                               <  1:701809/129‑1 (MQ=255)
                                                                                                             cGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTAGCTTCTGATGCGCAAGCTGAAGaaaaa                               >  2:701809/1‑129 (MQ=255)
                                                                                                                   gCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGt                                                                   >  1:328166/1‑87 (MQ=255)
                                                                                                                   gCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGt                                                                   <  2:328166/87‑1 (MQ=255)
                                                                                                                    cTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGaaaa                                >  3:377038/1‑121 (MQ=255)
                                                                                                                       aaCGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAgg                                                                        >  2:271477/1‑78 (MQ=255)
                                                                                                                       aaCGTCTTCTTTAATGGCTATCAACATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAgg                                                                        <  1:271477/78‑1 (MQ=255)
                                                                                                                        aCGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTa                                                                                 <  3:1043339/68‑1 (MQ=255)
                                                                                                                        aCGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTa                                                                                 >  4:1034274/1‑68 (MQ=255)
                                                                                                                        aCGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTa                                                                                 <  3:1034274/68‑1 (MQ=255)
                                                                                                                        aCGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGAt                                                 >  2:336196/1‑100 (MQ=255)
                                                                                                                           tcttctTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGAt                                                 <  4:385894/97‑1 (MQ=255)
                                                                                                                           tcttctTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTga                                      >  2:1082013/1‑108 (MQ=255)
                                                                                                                           tcttctTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTga                                      <  1:1082013/108‑1 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTgg                                                         <  2:716891/86‑1 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTgg                                                         >  1:716891/1‑86 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATgcg                                              >  2:161508/1‑97 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATgcg                                              <  1:161508/97‑1 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATg                                                <  1:98374/95‑1 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATg                                                >  2:98374/1‑95 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGa          <  2:53002/133‑1 (MQ=255)
                                                                                                                              tcttTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGa          >  1:53002/1‑133 (MQ=255)
                                                                                                                                 ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAg                                          <  2:970446/98‑1 (MQ=255)
                                                                                                                                 ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGa                            <  4:200728/112‑1 (MQ=255)
                                                                                                                                 ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGa                            >  3:200728/1‑112 (MQ=255)
                                                                                                                                 ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTa  <  2:1014097/138‑1 (MQ=255)
                                                                                                                                 ttAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTa  >  1:1014097/1‑138 (MQ=255)
                                                                                                                                              gCATGATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCtt                                                      >  1:416122/1‑73 (MQ=255)
                                                                                                                                              gCATGATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCtt                                                      <  2:416122/73‑1 (MQ=255)
                                                                                                                                                  catGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGaaaaa                               >  3:149914/1‑92 (MQ=255)
                                                                                                                                                  catGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGaaaaa                               <  4:149914/92‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                        
GTCTCGGCGTTGGCATGACTCCGACGCTCTCTCACGATTTCTTTTCTAAGTTACCGGCCGTACTGCAACCGCTGCTACATAGCGGCATTATGCTCGCAACCCTTAGCGCCGTTGTGCTGAACGTCTTCTTTAATGGCTATCAGCATCATGCTGACCTGGTGAAGGAATCCGTCTCTGATAAAGATTTAAAAGTCAGGACAGTACGTATGTGGCTTCTGATGCGCAAGCTGAAGAAAAATGAGCATGGAGAATAATATGAATTTTTTA  >  NC_000913/3032562‑3032828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: