Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,512,561 Δ1 bp 68.6% intergenic (‑217/+73) yhfZ ← / ← trpS putative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,512,5600T.68.6% 224.0 / 246.7 172intergenic (‑216/+74)yhfZ/trpSputative DNA‑binding transcriptional regulator/tryptophanyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base . (56/62);  total (85/87)
Fisher's exact test for biased strand distribution p-value = 5.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 179 total aligned reads displayed.

ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAA  >  NC_000913/3512461‑3512656
                                                                                                   |                                                                                                
aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt                                                                                                  >  1:2964412/1‑100 (MQ=255)
aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt                                                                                                  >  2:3275777/1‑100 (MQ=255)
aTAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAt                                                                                                  >  2:1519242/1‑100 (MQ=255)
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                               >  2:2739997/1‑97 (MQ=255)
   aaaaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCAtact                                                                                               >  2:2530660/1‑97 (MQ=255)
     aaaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATACGTGCATTACtt                                                                                             >  1:1384944/1‑100 (MQ=255)
      aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTg                                                                                           >  2:828304/1‑100 (MQ=255)
      aaaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTg                                                                                           <  1:3195424/100‑1 (MQ=255)
       cttccgatctCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGt                                                                                          <  2:173992/91‑1 (MQ=25)
       aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGt                                                                                          <  2:3160049/100‑1 (MQ=255)
       aaCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGt                                                                                          <  1:2693735/100‑1 (MQ=255)
        aCCGTGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTa                                                                                         <  2:3911663/100‑1 (MQ=255)
        aCCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTCATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTa                                                                                         <  1:4085318/100‑1 (MQ=255)
         ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTgaga                                                                                        >  2:1288208/1‑97 (MQ=255)
         ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAg                                                                                        >  1:1952057/1‑100 (MQ=255)
         ccGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAg                                                                                        <  1:3333715/100‑1 (MQ=255)
         ccGGGTTTCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAg                                                                                        >  1:697561/1‑99 (MQ=255)
            ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGcc                                                                                     <  1:1150260/100‑1 (MQ=255)
            ggTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGcc                                                                                     <  2:3620492/100‑1 (MQ=255)
             gTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGcc                                                                                     >  2:2925301/1‑100 (MQ=255)
                 ccccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCACTATGTGCATTACTTTGTAAGCCGGa                                                                                  >  2:3678788/1‑99 (MQ=255)
                 ccccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGAt                                                                                 >  2:3733/1‑100 (MQ=255)
                 ccccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAgatcggaag                                                                                >  1:173992/1‑92 (MQ=37)
                 ccccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATa                                                                                >  1:2263623/1‑100 (MQ=255)
                  cccGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATaa                                                                               >  1:3565985/1‑100 (MQ=255)
                    cGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAgg                                                                             <  2:616633/100‑1 (MQ=255)
                    cGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAgg                                                                             <  1:1277486/100‑1 (MQ=255)
                     ggTTTCAGAGTGATGATAAAAGCAAAATTTCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGc                                                                            <  1:920989/100‑1 (MQ=255)
                      gTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGc                                                                            >  1:3633283/1‑100 (MQ=255)
                         tCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTc                                                                        >  1:3508456/1‑100 (MQ=255)
                             agTGATGATAAAAGCTAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgc                                                                    <  2:4386409/100‑1 (MQ=255)
                             agTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgc                                                                    <  1:1446683/100‑1 (MQ=255)
                              gTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                   <  1:4138767/100‑1 (MQ=255)
                              gTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAACATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                   <  2:2268461/100‑1 (MQ=255)
                               tgatgaTAAAAGCAACATTGCCTGATGCGCTACGCGTATCAGGCGTACTTTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAAGCGTTCAcgcc                                                                   <  1:2542770/100‑1 (MQ=37)
                               tgatgaTAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCAcgcc                                                                   <  2:1329173/100‑1 (MQ=255)
                               tgatgaTAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                  <  2:1391070/100‑1 (MQ=255)
                                gatgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                  <  2:2149979/100‑1 (MQ=255)
                                gatgatAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCAcgccg                                                                  >  1:261952/1‑100 (MQ=255)
                                       aatgCAAAATTGCCTGATGCGCTACGCTTATCAGGTCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATGCGGCGTGGACGCCTTATCCGGc                                                          <  1:60901/97‑1 (MQ=21)
                                       aaaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                          <  2:1410645/100‑1 (MQ=255)
                                       aaaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                          >  2:3738399/1‑100 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                          <  1:2335514/100‑1 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                         >  1:2746960/1‑100 (MQ=255)
                                        aaGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                         <  1:214313/100‑1 (MQ=255)
                                          gCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                                                       <  1:2662964/100‑1 (MQ=255)
                                           cAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                      >  1:2423127/1‑100 (MQ=255)
                                           cAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                                                      >  2:2061284/1‑100 (MQ=255)
                                                 taccTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGGCGCATCCGGCATGAACAAAg                                                <  1:2684446/98‑1 (MQ=255)
                                                 tGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                                <  2:52658/100‑1 (MQ=255)
                                                 tGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                                <  2:3489049/100‑1 (MQ=255)
                                                 tGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAg                                                <  2:1125686/100‑1 (MQ=255)
                                                  gCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCCATATGTGCATACTTTGTAGGCCGGACAAAGCGTTCACGCCGCATCCGACATGAACACAgc                                               >  1:3783222/1‑100 (MQ=25)
                                                  gCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAgc                                               <  2:166187/100‑1 (MQ=255)
                                                   ccTGATGCGCTACGCTTATCTGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTACGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACaaa                                                 >  1:994006/1‑98 (MQ=255)
                                                    cTGATGCGCTACGCTTATCAAGCCTACATCTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACgcgc                                             >  1:873297/1‑100 (MQ=21)
                                                         gcgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAttt                                        <  2:4275986/100‑1 (MQ=255)
                                                          cgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAttt                                        >  2:4178763/1‑100 (MQ=255)
                                                          cgcTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAAttt                                        >  2:4178865/1‑100 (MQ=255)
                                                            cTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAGGCGCAATTTGc                                      <  1:2316717/99‑1 (MQ=255)
                                                             tACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCa                                    >  1:3308480/1‑100 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                   >  2:3490735/1‑100 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAg                                   <  1:4178732/100‑1 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                  <  2:3536081/100‑1 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                  <  2:382395/100‑1 (MQ=255)
                                                               cGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                  >  2:237470/1‑100 (MQ=255)
                                                               cGCTTATAAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGTATAAGGCGTTCACGCCGCATCCTGCATGAACAAAGCGCAATTTGc                                      >  1:1045295/1‑97 (MQ=255)
                                                               agcgcaTCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                  <  1:224128/95‑1 (MQ=255)
                                                                gCTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGc                                  <  1:3881090/99‑1 (MQ=255)
                                                                 cTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCa                                 <  1:3211122/100‑1 (MQ=255)
                                                                 cTTATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCaa                                >  1:566906/1‑100 (MQ=255)
                                                                   tATCAGGCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATa                              >  1:2968922/1‑100 (MQ=255)
                                                                     tCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAg                             <  1:310742/100‑1 (MQ=255)
                                                                        gtccTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGa                          <  2:2275514/98‑1 (MQ=255)
                                                                         gCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAt                        >  1:1819238/1‑100 (MQ=255)
                                                                         gCCTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAAt                        <  2:1360386/100‑1 (MQ=255)
                                                                          ccTACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCAGCATTCGGCATGAAAAAAGCGCAATTTGCCAGCAATAGTGAAtt                       >  1:1875829/1‑100 (MQ=25)
                                                                            tACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTAc                     >  2:2460880/1‑100 (MQ=255)
                                                                            tACATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTAc                     >  2:677964/1‑100 (MQ=255)
                                                                             aCATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACg                    <  2:2218390/100‑1 (MQ=255)
                                                                              cATTTCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACgg                   <  1:1519236/100‑1 (MQ=255)
                                                                                tttCCTTGCAATATGTGCATACTTTGTAGGCCCGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACg                    <  1:4045229/97‑1 (MQ=255)
                                                                                  tCCTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTc               >  1:2834324/1‑100 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAATGCGTTCACGCCTCATACGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc             <  1:288143/100‑1 (MQ=25)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc             >  1:1700801/1‑100 (MQ=255)
                                                                                    cTTGCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc             <  1:632527/100‑1 (MQ=255)
                                                                                     ttGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGc             >  1:2543309/1‑100 (MQ=255)
                                                                                      tGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGcc            >  2:644044/1‑100 (MQ=255)
                                                                                       gCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCca           <  2:2607488/100‑1 (MQ=255)
                                                                                       gCAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGcccc          <  2:2968939/100‑3 (MQ=255)
                                                                                        cAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcac          <  2:728633/100‑1 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca         <  2:3976961/100‑1 (MQ=255)
                                                                                        cAATATGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca         <  2:2556387/100‑1 (MQ=255)
                                                                                         aaTATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCcaca         <  2:2385400/100‑1 (MQ=255)
                                                                                          atatGTGCATACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCTCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACaaa       <  1:3527136/100‑1 (MQ=255)
                                                                                            atGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACaaaa      >  2:3055328/1‑100 (MQ=255)
                                                                                             tgtgCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAAc     <  1:226484/100‑1 (MQ=255)
                                                                                               tgCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCa   <  1:2072613/100‑1 (MQ=255)
                                                                                                gCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCaa  >  2:1443166/1‑100 (MQ=255)
                                                                                                gCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCaa  >  1:443192/1‑100 (MQ=255)
                                                                                                   |                                                                                                
ATAAAAAAACCGGGTTTCCCCGGTTTCAGAGTGATGATAAAAGCAAAATTGCCTGATGCGCTACGCTTATCAGGCCTACATTTCCTTGCAATATGTGCATTACTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCGCAATTTGCCAGCAATAGTGAATTACGGCTTCGCCACAAAACCAA  >  NC_000913/3512461‑3512656

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: