Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,160,853 C→T intergenic (‑64/‑270) sthA ← / → fabR pyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,160,8530CT95.7% 190.4 / NA 71intergenic (‑64/‑270)sthA/fabRpyridine nucleotide transhydrogenase, soluble/transcriptional repressor of fabA and fabB
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (31/36);  minor base G (3/0);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAA  >  NC_000913/4160714‑4160992
                                                                                                                                           |                                                                                                                                           
ccAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCaa                                                                                                                                          >  1:324271/1‑143 (MQ=255)
 cAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTCTACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAAt                                                                                                                                         <  1:631215/143‑1 (MQ=255)
  aGGCCCATTGCTGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAAtt                                                                                                                                        <  2:719965/143‑1 (MQ=255)
    gCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTCTACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGGCAAttca                                                                                                                                      >  2:444771/1‑141 (MQ=255)
    gCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGa                                                                                                                                      >  2:169241/1‑143 (MQ=255)
     cccATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGAt                                                                                                                                     <  2:894477/143‑1 (MQ=255)
            cAGCGCCTTCGCCGCCGGGGCCGGGACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTc                                                                                                                              <  2:389646/143‑1 (MQ=255)
            cAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAGGTCCATTGATCCCCGTc                                                                                                                              >  1:1055118/1‑143 (MQ=255)
            cAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATGTTATCGCGGGTAAGTCAATTGATCACAGTc                                                                                                                              >  1:622453/1‑143 (MQ=255)
              gcgcCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCAt                                                                                                                            <  1:1093163/143‑1 (MQ=255)
              gcgcCTTCGCCGCCGGGGCCGGAACCTAATACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCAt                                                                                                                            <  2:339675/143‑1 (MQ=255)
                  ctgcgccgcCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTAt                                                                                                                        <  2:105341/140‑1 (MQ=255)
                  cTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCACTTGATCCGCGTCATTTAt                                                                                                                        >  2:49130/1‑143 (MQ=255)
                   ttCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATc                                                                                                                       <  1:798269/143‑1 (MQ=255)
                      gctgccgGGGCCGGAACCTATTATTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTATCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTgg                                                                                                                    <  2:649085/140‑1 (MQ=255)
                      gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTATGTCAATTGATCCCCGTCATTTATCTgg                                                                                                                    <  2:141304/143‑1 (MQ=255)
                      gccgccGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTgg                                                                                                                    <  2:747334/143‑1 (MQ=255)
                       gcgctgGGGGCGGAGACTCTTAATACGGCATCGTAATCGTAGGAATGTGGCAGGGTAGGGCTTACCTGTTCTGATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGc                                                                                                                   <  1:629911/142‑1 (MQ=255)
                       ccgccgGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATGTTATCGCGGGTAAGTGAATTGATCCCCGTCATTTATCTGGc                                                                                                                   >  2:100786/1‑143 (MQ=255)
                        cgccgGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCt                                                                                                                  <  1:862969/143‑1 (MQ=255)
                         gccgGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCta                                                                                                                 <  1:550156/143‑1 (MQ=255)
                            gggtgcGGAATCTGTTACTATGGCAGCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTatat                                                                                                              <  2:390352/138‑1 (MQ=255)
                               gCCGGAGCCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCt                                                                                                           <  2:686271/143‑1 (MQ=255)
                                           tgctagtgCATCGTAATCGCAGGAATGTGGCATGGTAGGGCTTACCTGTTCTGATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGGAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttg                                                                                               <  2:951745/136‑1 (MQ=255)
                                            aCTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATACCCGTCGTTTATCTGGGTATATGCTGAGCGGCcgttgc                                                                                              >  1:197173/1‑143 (MQ=255)
                                             cTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgct                                                                                             >  1:1074306/1‑143 (MQ=255)
                                              tATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCctttgctt                                                                                            >  1:514599/1‑143 (MQ=255)
                                               gtGGCATCGTAATCGGAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCcgttgcttt                                                                                           <  2:442969/142‑1 (MQ=255)
                                                    atcgtaatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAAGGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCCATACCCTGAGCGACCTCTGCTGTGTCTg                                                                                      >  2:148935/1‑143 (MQ=255)
                                                     tcgtatccctaGTAATGGGGCATGGTCGGGCTTAGCTGTTCTTGCCCCTCCAGCCAGCAGAATGGTAACCTCGTATCGCGGCTCAGTCAATTGCTCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGt                                                                                     <  2:511032/134‑1 (MQ=255)
                                                      cgtaatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGtt                                                                                    >  1:1071173/1‑143 (MQ=255)
                                                      cgtaatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTATTATACATAAAAGCAACAGAATGGTAACATTTTATAGCGGGAAAGTCAATTGATCCGCGTTGGTTGTCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGtt                                                                                    >  1:718731/1‑143 (MQ=255)
                                                       gtaatggtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGttt                                                                                   <  2:363899/143‑1 (MQ=255)
                                                        taatcgtaGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCACGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTc                                                                                  >  2:167761/1‑143 (MQ=255)
                                                           tcgtaGGAATGTGGCATGGTAGTGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTa                                                                               <  2:77903/143‑1 (MQ=255)
                                                              taGGAATGCGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGGCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCGCTTGCTTTGTCGGTTTGTTGCtt                                                                            >  2:315113/1‑143 (MQ=255)
                                                              taGGAACGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCCTTTATCTGGCTCTCTCCTTAGCGGCCTTTGCTTTGTCTCTTTCTTACtt                                                                            >  1:299550/1‑143 (MQ=255)
                                                                ggAATGGGGCATGGTAGGGCTTACCTGTTCCTATACATAAAAGCAACAGAATGGTAACATTTTCTCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACtttt                                                                          <  1:659365/143‑1 (MQ=255)
                                                                 gAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATATATCTGGCCATAACCTGAGCGGCCTTGGCTTTGTCTGTTTCTTAGTTTTg                                                                         >  2:392533/1‑143 (MQ=255)
                                                                  aagggggtaTGGTAGCGCTTACCTGTTCTTATGCATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGTTCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGc                                                                        <  2:764761/135‑1 (MQ=255)
                                                                  aaTGTGGCATGGGAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGc                                                                        <  2:766586/143‑1 (MQ=255)
                                                                      tgGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGACTGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTg                                                                    <  2:888415/143‑1 (MQ=255)
                                                                      tgGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTACGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTg                                                                    >  1:1153930/1‑143 (MQ=255)
                                                                      tgGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATTTCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTg                                                                    >  1:784272/1‑143 (MQ=255)
                                                                                gggCTTACCTGTTCTTATACATAACAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCGCCGTCTTGTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTAtt                                                          <  2:1003033/143‑1 (MQ=255)
                                                                                gggCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTTTATCCTGAGCGGCCTTTGCATTGTCTGTTGCTTACTTTTGCCCTCAGGTTTTAta                                                          >  2:786677/1‑142 (MQ=255)
                                                                                 ggCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATATGGCTATATCCTGAGCGGCCTCTGCTTTGTCTGTTTCTTACTTTTGCCCTGCCGCCTTATTg                                                         >  2:228541/1‑143 (MQ=255)
                                                                                   cTTACCTGTCCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTTACGTTTTATTGGa                                                       <  1:75355/143‑1 (MQ=255)
                                                                                    ttACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAt                                                      >  2:628682/1‑143 (MQ=255)
                                                                                      aCCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCACTTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAttt                                                    <  1:897593/143‑1 (MQ=255)
                                                                                        cTGTTCTTATACATAAAAGCAACAGAATGGTAACAGTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGAttttt                                                  <  1:55975/143‑1 (MQ=255)
                                                                                                ataCATAAAAGCAACAGAATGGTACCATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGTCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATcgacga                                          <  1:820522/143‑1 (MQ=255)
                                                                                                ataCATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATcgacga                                          <  1:1150109/143‑1 (MQ=255)
                                                                                                    aTAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATAct                                      >  1:566220/1‑143 (MQ=255)
                                                                                                    aTAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCATCATATATCTGGCTATATCCTGAGCGGCCTTTGCATTGTCTGTTTCTTACTTTTGGCATGCCGTTTTATTGGATTTTTATCTAGTATAct                                      >  2:260093/1‑143 (MQ=255)
                                                                                                        aagtagcAGAATGGTAACATTTTATCGCGGGTAGGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTcc                                  <  2:657485/137‑1 (MQ=255)
                                                                                                         aGCAACAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTTTCGACGATattgtgcg                                 >  2:971949/1‑136 (MQ=255)
                                                                                                            caccacggggTGACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTAGATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACCTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCCttt                              >  2:199818/9‑143 (MQ=255)
                                                                                                             aCAGAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTCTTGGATCTTTATCGACGATCCTCTCCGTTTa                             >  1:8788/1‑143 (MQ=255)
                                                                                                              cAGAATGGTAACATTAGATGGCGGTGAGGTCAATTGATCCCCGTCATTTGTCTGGCTATATCCTGAGCGGCGTGTGCTTTGTCTGTTTCTTACTTTTGCGCTGACGTTTTCTTGGATTGTTGTGGACGATACTCTCGGGTTaa                            >  2:1010177/1‑143 (MQ=255)
                                                                                                                gAATGGTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGc                          >  2:224933/1‑143 (MQ=255)
                                                                                                                    ggTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGCTTAACCGGCa                      >  2:921122/1‑143 (MQ=255)
                                                                                                                     gtaaacttttATCGCGGGTAAGTCAATTGCTCCCCGTCCGTTCTTTGTTGTTATCCGGAGCGGCCTTTGCCTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAg                     <  1:315113/137‑1 (MQ=255)
                                                                                                                     gTAACATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTCTATCCTGAACGCCCTTCGCTTTGTCTGTTTCTTACTTTTGCCCCGACGCTTTATTGCCTTTTTAGTGCCGATACTCTCCGTTTCAGCgaccg                     >  2:953200/1‑139 (MQ=255)
                                                                                                                       aaCATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGt                   <  1:188115/143‑1 (MQ=255)
                                                                                                                        aCATTTTATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGCCTGTTTCTTACTTTTGCCCTGACGCTTTACTGGATTTCTGTCGACGATACTCTCCGTTTaacgcgcagctt                  >  2:25989/1‑133 (MQ=255)
                                                                                                                         cATTTTACAGGGGGTAAGGCAATTGATCCCCGTCATTTATATGGCTATATACTGAGCAGCGTTTTCTTTGTCTGTTTCTTAGTTTTGCTCTAACGTTTTTTTGTTTTTTTATCGCCTATCCGCCGCTTTTATGCAGCAGGttt                 >  2:199551/1‑143 (MQ=255)
                                                                                                                             ttATCGCGGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGc             <  2:139122/143‑1 (MQ=255)
                                                                                                                                   cgGGTAAGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGCTTAAGCGGCAGGTCTCCGCTGTCCg       >  1:114411/1‑143 (MQ=255)
                                                                                                                                      gTAAGCTAATTCCTCCCCTTCATTTATCTGGCTATATGCTGAGCGGCTTTTGCTTTGTCTCTTTCTTACTTCTGCCCGCACGTTTGAGTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaa    <  2:1062883/143‑1 (MQ=255)
                                                                                                                                        aaGTCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTaaaa  <  1:169241/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CCAGGCCCATTGCAGCGCCTTCGCCGCCGGGGCCGGAACCTATTACTATGGCATCGTAATCGTAGGAATGTGGCATGGTAGGGCTTACCTGTTCTTATACATAAAAGCAACAGAATGGTAACATTTTATCGCGGGTAAGCCAATTGATCCCCGTCATTTATCTGGCTATATCCTGAGCGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAA  >  NC_000913/4160714‑4160992

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