Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 944,073 C→T 5.9% T14M (ACG→ATG)  dmsC → dimethyl sulfoxide reductase, anaerobic, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913944,0730CT 5.9% 217.5 / 2.9 68T14M (ACG→ATG) dmsCdimethyl sulfoxide reductase, anaerobic, subunit C
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base T (2/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTAT  >  NC_000913/943927‑944220
                                                                                                                                                  |                                                                                                                                                   
ccGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACg                                                                                                                                                    <  2:323640/148‑1 (MQ=255)
     aGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTc                                                                                                                                              <  4:37736/149‑1 (MQ=255)
                     ccGTATATTGTGATCAAACCGAATGCCAATAGCCGCCCGACCGTGTATACCACTGGCTATCTGGCAAAGCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAg                                                                                                                              <  2:216056/149‑1 (MQ=255)
                      cGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAgg                                                                                                                             <  5:396353/149‑1 (MQ=255)
                          tatTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAggtggt                                                                                                                         >  8:65113/1‑149 (MQ=255)
                          tatTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAggtggt                                                                                                                         >  6:222756/1‑149 (MQ=255)
                             tgtgATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGttt                                                                                                                       >  8:161669/1‑148 (MQ=255)
                               tgATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTa                                                                                                                     <  2:405965/148‑1 (MQ=255)
                                 aTCAAACCCAATGCCAATAGGCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATc                                                                                                                   <  7:76471/148‑1 (MQ=255)
                                    aaaCCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTc                                                                                                               >  1:361003/1‑149 (MQ=255)
                                        ccAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTAt                                                                                                                    >  1:413286/1‑140 (MQ=255)
                                        ccAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTAt                                                                                                                    <  2:413286/140‑1 (MQ=255)
                                             gCCAATAGCCGCCCGCCCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCCCGGTCTTCGGGCAATGTGTAGCAGGGGGTTTTATCGTTCTGGCTTTgg                                                                                                      >  6:329105/1‑149 (MQ=255)
                                                aaTAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATgtgt                                                                                                                                   <  6:209570/117‑1 (MQ=255)
                                                aaTAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATgtgt                                                                                                                                   >  5:209570/1‑117 (MQ=255)
                                                       gccCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCa                                                                                             >  2:141529/1‑148 (MQ=255)
                                                          cGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAgg                                                                                         >  4:44141/1‑149 (MQ=255)
                                                                gggATACAACTGGCTATCTGGCAAACCCTAAGGAGGAGTGAGATGGGAAGTGGATGGGATGAATGGCCGCTGATGATCTTCACTGTCTTCGGGCAATGAGTATCAGGTGGTTTTATCGTGATGGCTTTGGCGCTGCTCAAAGGCGACCt                                                                                   >  2:354414/1‑149 (MQ=255)
                                                                  gATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGc                                                                                 <  1:141529/149‑1 (MQ=255)
                                                                     aCCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGgcgc                                                                                                   >  8:248810/1‑128 (MQ=255)
                                                                     aCCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGgcgc                                                                                                   <  7:248810/128‑1 (MQ=255)
                                                                     aCCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAg                                                                              >  5:126408/1‑149 (MQ=255)
                                                                      cgaCTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTgg                                                                                                            <  6:60485/116‑1 (MQ=255)
                                                                        aCTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAg                                                                              >  3:261327/1‑146 (MQ=255)
                                                                             cTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      <  6:126408/149‑1 (MQ=255)
                                                                              tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAAtgt                                                                                                                                     <  6:146658/85‑1 (MQ=255)
                                                                              tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAAtgt                                                                                                                                     >  5:146658/1‑85 (MQ=255)
                                                                              tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      >  4:314981/1‑148 (MQ=255)
                                                                              tATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      >  2:188497/1‑148 (MQ=255)
                                                                               aTCTGGCAAACCCGAAGGAGGTGTTAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGAGCTGCTCAAAGGCTCCCTGCGAGAAGAAGcc                                                                      >  6:73172/1‑147 (MQ=255)
                                                                               aTCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGcc                                                                      >  5:158467/1‑147 (MQ=255)
                                                                                tCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCca                                                                    <  7:65113/148‑1 (MQ=255)
                                                                                tCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCATGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCca                                                                    <  5:222756/148‑1 (MQ=255)
                                                                                    gCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCagcagc                                                               >  3:155080/1‑149 (MQ=255)
                                                                                     cAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCg                                                              <  2:361003/149‑1 (MQ=255)
                                                                                      aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAggtg                                                                                                                           <  7:236188/87‑1 (MQ=255)
                                                                                      aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAggtg                                                                                                                           >  8:236188/1‑87 (MQ=255)
                                                                                      aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCgt                                                             >  4:324695/1‑149 (MQ=255)
                                                                                      aaaCCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCgt                                                             >  4:251716/1‑149 (MQ=255)
                                                                                         cccGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCgtgt                                                           >  1:39192/1‑148 (MQ=255)
                                                                                           cGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGtt                                                                                                                        >  6:308924/1‑85 (MQ=255)
                                                                                           cGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGtt                                                                                                                        <  5:308924/85‑1 (MQ=255)
                                                                                              aggaggTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGc                                                     <  3:251716/149‑1 (MQ=255)
                                                                                                      tgAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCcag                                                                   >  1:42724/1‑127 (MQ=255)
                                                                                                      tgAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCcag                                                                   <  2:42724/127‑1 (MQ=255)
                                                                                                      tgAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGt                                             >  1:36342/1‑149 (MQ=255)
                                                                                                           tGGGATGTGGATGGCATGAATGGCCGCTGATGATCTTCACTGTATTCTGGCAATGTGTAGCAGGTGGGTTTATCGTTCTGGCTGTGGGGCTGCTCAAAGGAGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTgg                                         <  2:260287/148‑1 (MQ=255)
                                                                                                            gggAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGtt                                       >  2:348945/1‑149 (MQ=255)
                                                                                                             ggAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGtt                                       <  1:348945/148‑1 (MQ=255)
                                                                                                              gAAGTGGCTGGCATGAATGGCCGCTGATGTTCTTCTCGGGCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGcta                                                                               <  3:218883/107‑3 (MQ=255)
                                                                                                                 gTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATgg                                  <  3:44141/149‑1 (MQ=255)
                                                                                                                    gATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTAAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGt                                >  5:110427/1‑148 (MQ=255)
                                                                                                                         cATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGAt                          >  3:49173/1‑149 (MQ=255)
                                                                                                                           ttaaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCaaat                                                                                          >  3:36675/3‑82 (MQ=255)
                                                                                                                           ttaaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAg                                                                                          <  4:36675/81‑1 (MQ=255)
                                                                                                                           tGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATgg                        <  2:39192/149‑1 (MQ=255)
                                                                                                                             aaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCa                                                <  8:381183/123‑1 (MQ=255)
                                                                                                                             aaTGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCa                                                >  7:381183/1‑123 (MQ=255)
                                                                                                                                ggCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCAtt                   >  6:335871/1‑149 (MQ=255)
                                                                                                                                  ccGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTgg                 <  4:49173/149‑1 (MQ=255)
                                                                                                                                        atgatCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGt                                        >  6:254123/1‑120 (MQ=255)
                                                                                                                                        atgatCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGt                                        <  5:254123/120‑1 (MQ=255)
                                                                                                                                        atgatCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTAt           <  5:335871/149‑1 (MQ=255)
                                                                                                                                         tgatCTTCACGGTATGCGCTCAATTTGTAGCAGGTGGTTTTATCGTTCTTGCTTTGGCGCTGCTCAAAGGCGACATGCGAGCAGAAGCACAGCATCGTGTTATCGCCTTCATGTTTGTTTTATGGGTGCTGATGGGCATTGGCTTTATc          <  5:73172/149‑1 (MQ=255)
                                                                                                                                            tCTTCACGGTCTTCGGGCATTGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGACGGCATGTTTGGTTTATGGGTGCTGATGGGCATTGTCTTTATCGcc       >  6:95832/1‑149 (MQ=255)
                                                                                                                                                cACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGc        <  3:324695/144‑1 (MQ=255)
                                                                                                                                                 aCGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTa   <  4:155080/148‑1 (MQ=255)
                                                                                                                                                 aCGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTAt  >  3:343958/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CCGCGAGCTCACTTTACCAAACCGAATATTGTGATCAAACCCAATGCCAATAGCCGCCCGACCGGGGATACCACTGGCTATCTGGCAAACCCGAAGGAGGTGTGAGATGGGAAGTGGATGGCATGAATGGCCGCTGATGATCTTCACGGTCTTCGGGCAATGTGTAGCAGGTGGTTTTATCGTTCTGGCTTTGGCGCTGCTCAAAGGCGACCTGCGAGCAGAAGCCCAGCAGCGTGTTATCGCCTGCATGTTTGGTTTATGGGTGCTGATGGGCATTGGCTTTATCGCCTCTAT  >  NC_000913/943927‑944220

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