Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,424,457 (A)8→9 intergenic (+71/+59) ulaF → / ← yjfY L‑ribulose‑5‑phosphate 4‑epimerase UlaF/DUF1471 domain‑containing protein YjfY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,424,4491.A98.6% 884.1 / 3.8 215intergenic (+63/+67)ulaF/yjfYL‑ribulose‑5‑phosphate 4‑epimerase UlaF/DUF1471 domain‑containing protein YjfY
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base A (122/90);  total (123/92)
Fisher's exact test for biased strand distribution p-value = 5.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 266 total aligned reads displayed.

GTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCA  >  NC_000913/4424349‑4424544
                                                                                                      |                                                                                                
gTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTc                                                                                                   <  1:5849766/101‑1 (MQ=255)
ctAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTc                                                                                                   >  1:5850207/2‑101 (MQ=255)
 tAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCa                                                                                                 <  1:5849667/101‑1 (MQ=255)
 tAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCa                                                                                                 >  2:5850211/1‑101 (MQ=255)
  aaaCACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaa                                                                                                >  2:5850213/1‑101 (MQ=255)
  aaaCACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaa                                                                                                >  1:5850214/1‑101 (MQ=255)
   aaCACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaa                                                                                              >  2:5850217/1‑101 (MQ=255)
      acGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaa                                                                                           >  2:5850226/1‑101 (MQ=255)
      acGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaa                                                                                           >  1:5850225/1‑101 (MQ=255)
      acGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaa                                                                                           >  1:5850222/1‑101 (MQ=255)
       cGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaa                                                                                          <  1:5849771/101‑1 (MQ=255)
       cGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaa                                                                                          >  2:5850228/1‑101 (MQ=255)
        ggTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaa                                                                                         >  1:5850232/1‑101 (MQ=255)
          tCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaaag                                                                                       >  2:5850235/1‑99 (MQ=255)
           ccTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaaagc                                                                                      <  1:5849847/101‑4 (MQ=255)
            cTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaaagcc                                                                                     <  2:5850012/101‑5 (MQ=255)
            cTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaaagcc                                                                                     >  1:5850237/1‑97 (MQ=255)
             tAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaaagccc                                                                                    <  1:5850072/101‑6 (MQ=255)
             tAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaaagccc                                                                                    <  1:5849858/101‑6 (MQ=255)
               aCGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGccccc                                                                                   >  2:5850242/1‑101 (MQ=255)
                cGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCt                                                                                  <  1:5849870/101‑1 (MQ=255)
                cGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCt                                                                                  <  1:5849901/101‑1 (MQ=255)
                cGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCt                                                                                  <  1:5850033/101‑1 (MQ=255)
                cGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCt                                                                                  >  1:5850243/1‑101 (MQ=255)
                 gCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTg                                                                                 >  2:5850245/1‑101 (MQ=255)
                   ttattaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCg                                                                               >  2:5850249/1‑101 (MQ=255)
                      ttaCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAg                                                                            >  1:5850252/1‑101 (MQ=255)
                       taCGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAgg                                                                           >  1:5850254/1‑101 (MQ=255)
                         cGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAgggg                                                                         >  2:5850258/1‑101 (MQ=255)
                          gggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCATCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAggggg                                                                        <  2:5849848/101‑1 (MQ=255)
                           ggCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGgc                                                                       >  1:5850260/1‑101 (MQ=255)
                             cAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGgcaa                                                                     >  2:5850263/1‑101 (MQ=255)
                               gAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                                                   >  1:5850266/1‑101 (MQ=255)
                               gAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                                                   <  2:5849960/101‑1 (MQ=255)
                                aaGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGgcaaagc                                                                  >  2:5850268/1‑101 (MQ=255)
                                   tAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGgcaaagcaaa                                                               >  1:5850272/1‑101 (MQ=255)
                                    aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAc                                                              >  2:5850275/1‑101 (MQ=255)
                                    aGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAAc                                                              >  1:5850276/1‑101 (MQ=255)
                                      aaCACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCTCTGCGAAGGGGGCAAAGCAAACTa                                                            <  2:5849779/101‑1 (MQ=255)
                                        cacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAgg                                                                           >  2:5850281/1‑84 (MQ=255)
                                        cacGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAGAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATg                                                          >  1:5850280/1‑101 (MQ=255)
                                             gcTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAAt                                                     >  1:5850284/1‑101 (MQ=255)
                                               tGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGt                                                   >  2:5850288/1‑101 (MQ=255)
                                               tGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCaaaaaaaaagc                                                                                      >  1:5850289/1‑62 (MQ=255)
                                                  ggAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTc                                                <  1:5849812/101‑1 (MQ=255)
                                                   gAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                                               <  2:5849873/101‑1 (MQ=255)
                                                   gAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCg                                               >  2:5850294/1‑101 (MQ=255)
                                                    aaaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGt                                              >  2:5850296/1‑101 (MQ=255)
                                                     aaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGtt                                             <  2:5850025/101‑1 (MQ=255)
                                                     aaTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGtt                                             <  2:5849939/101‑1 (MQ=255)
                                                      aTTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTg                                            >  2:5850299/1‑101 (MQ=255)
                                                       tttCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTgg                                           >  1:5850300/1‑101 (MQ=255)
                                                          ccttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTta                                        >  1:5850304/1‑101 (MQ=255)
                                                           cttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat                                       >  2:5850305/1‑101 (MQ=255)
                                                           cttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat                                       <  2:5850021/101‑1 (MQ=255)
                                                           cttcctCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTtat                                       <  2:5849948/101‑1 (MQ=255)
                                                                tCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGt                                          >  1:5850310/1‑93 (MQ=255)
                                                                tCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTg                                  >  2:5850309/1‑101 (MQ=255)
                                                                  gggAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                >  2:5850314/1‑101 (MQ=255)
                                                                  gggAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGt                                <  2:5849878/101‑1 (MQ=255)
                                                                     agaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCt                             >  2:5850318/1‑101 (MQ=255)
                                                                      gaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTa                            >  2:5850321/1‑101 (MQ=255)
                                                                      gaTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTa                            >  1:5850320/1‑101 (MQ=255)
                                                                       aTAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAg                           <  1:5849935/101‑1 (MQ=255)
                                                                            tGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATa                      <  2:5849826/101‑1 (MQ=255)
                                                                            tGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATa                      >  1:5850325/1‑101 (MQ=255)
                                                                             ggTCTAATTCCGCATCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATaa                     <  2:5849973/101‑1 (MQ=255)
                                                                              gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATaaa                    >  2:5850328/1‑101 (MQ=255)
                                                                              gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATaaa                    <  1:5849951/101‑1 (MQ=255)
                                                                              gTCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATaaa                    <  1:5849996/101‑1 (MQ=255)
                                                                               tCTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAat                   >  1:5850329/1‑101 (MQ=255)
                                                                                cTAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAata                  <  1:5849822/101‑1 (MQ=255)
                                                                                 tAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATGCCTGGTGCTAGCGATAAatat                 <  1:2605759/101‑1 (MQ=255)
                                                                                 tAATTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAatat                 >  2:5850333/1‑101 (MQ=255)
                                                                                   aTTCCGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATcc               >  2:5850335/1‑101 (MQ=255)
                                                                                      ccGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCcgcg            >  1:5850340/1‑101 (MQ=255)
                                                                                      ccGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCcgcg            >  2:5850339/1‑101 (MQ=255)
                                                                                       cGCAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCgcgc           >  1:5850342/1‑101 (MQ=255)
                                                                                         ccgCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGttt                                                 <  2:5849934/59‑1 (MQ=255)
                                                                                         catccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGc                          <  2:5849662/81‑1 (MQ=255)
                                                                                         cAGCCTTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGgcaaag                                                                   <  2:5849765/43‑1 (MQ=37)
                                                                                         cAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTg         >  1:5850346/1‑101 (MQ=255)
                                                                                         cAGCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTg         >  2:5850348/1‑101 (MQ=255)
                                                                                           gCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTgg                                 <  2:5849754/75‑1 (MQ=255)
                                                                                           gCCGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGc       >  1:5850352/1‑101 (MQ=255)
                                                                                            ccGTTTTTAAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTAt                                                           <  2:5849737/48‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGtt                                                  <  2:5849688/57‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTg                                            <  1:5849799/63‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGtt                                         <  1:5849851/66‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTgg                                 <  2:5849797/74‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCg                         <  2:5849788/82‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGa                        <  1:5849892/83‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTg         <  2:5850081/98‑1 (MQ=255)
                                                                                            ccGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCg      >  2:5850354/1‑101 (MQ=255)
                                                                                             cGTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGt     >  1:5850355/1‑101 (MQ=255)
                                                                                              gTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGaaa                                                      >  2:5850362/1‑51 (MQ=255)
                                                                                              gTTTTTCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTg    >  1:5850360/1‑101 (MQ=255)
                                                                                               tttttCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCt          <  1:5849920/94‑1 (MQ=255)
                                                                                               tttttCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGc   <  2:5849900/101‑1 (MQ=255)
                                                                                                ttttCAAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTg                               >  1:5850366/1‑72 (MQ=255)
                                                                                                      |                                                                                                
GTAAACACGGTCCTAACGCTTATTACGGGCAGAAGTAGAACACGCGCTGCGGAAATTTCCTTCCTCGGGAGATAACTGGTCTAATTCCGCAGCCGTTTTTCAAAAAAAAGCCCCCTGCGAAGGGGGCAAAGCAAACTATGGCAATGTTTCGTTGGTTATACCTGGTGCTAGCGATAAATATCCGCGCTGGCGTGCA  >  NC_000913/4424349‑4424544

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