Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,085,667 T→C 5.1% Y226H (TAT→CAT)  phoH → ATP‑binding protein; putative PhoH family P‑loop ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,085,6670TC 5.1% 244.8 / 2.7 78Y226H (TAT→CAT) phoHATP‑binding protein; putative PhoH family P‑loop ATPase
Reads supporting (aligned to +/- strand):  ref base T (35/39);  new base C (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

AAAAACCTGGATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGCAGTCGTCATTCTTGA  >  NC_000913/1085525‑1085815
                                                                                                                                              |                                                                                                                                                    
aaaaaCCTGGATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTc                                                                                                                                                >  3:127763/1‑149 (MQ=255)
aaaaaCCTGGATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTc                                                                                                                                                <  6:304818/149‑1 (MQ=255)
  aaaCCTGGATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGc                                                                                                                                              <  6:355785/149‑1 (MQ=255)
    aCCTGGATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGccc                                                                                                                                            >  8:133183/1‑149 (MQ=255)
       tGGATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGt                                                                                                                                         >  6:335495/1‑149 (MQ=255)
         gATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAACGATGTCGACAGGATTATCGGCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATTTCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTc                                                                                                                                        <  3:281413/148‑1 (MQ=255)
                 cagcaAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTg                                                                                                                               <  2:168816/149‑1 (MQ=255)
                  agcaAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGc                                                                                                                              <  2:189832/149‑1 (MQ=255)
                       aaGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTATGCAATCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTc                                                                                                                         <  5:208239/149‑1 (MQ=255)
                           ggCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGgc                                                                                                                     >  1:362328/1‑149 (MQ=255)
                             cAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGgcgc                                                                                                                   <  7:368089/149‑1 (MQ=255)
                               gagGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCATGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGGCTACGACGTGCTGGTCCGGCGCtt                                                                                                                 <  3:252969/149‑1 (MQ=255)
                                agGCGCTGATTCATAAGGATGTCGACTGGTTTATCGTCTCCCGTCCAGTTCTGCAAGGCGATGAAGATCTTGGCTTCTTTCCTGGAGATATGGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCtt                                                                                                                 <  2:362328/148‑1 (MQ=255)
                                 ggCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAg                                                                                                               <  7:133183/149‑1 (MQ=255)
                                 ggCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAg                                                                                                               <  4:127763/149‑1 (MQ=255)
                                  gCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAgg                                                                                                              <  6:106911/149‑1 (MQ=255)
                                   cccTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGGCTACGACGTGCTGGTCCGGCGCTTAggg                                                                                                             <  6:75455/149‑1 (MQ=255)
                                    ccTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAgggg                                                                                                            >  8:167880/1‑149 (MQ=255)
                                            aTAAGTATGTCGACAGGATTAGCGTCACCCGTCCATTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCATATATTCGCCCGGTCTACGACGTGCTGGTCAGGCGCTTAGTGGCTTCCttt                                                                                                    <  4:311639/149‑1 (MQ=255)
                                             tAAGGATGTCTACAGGATTATCGTCACCCGTCCAGGTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTa                                                                                                   <  6:88547/149‑1 (MQ=255)
                                             tAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAg                                                                                                               <  6:301715/137‑1 (MQ=255)
                                             tAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAg                                                                                                               >  5:301715/1‑137 (MQ=255)
                                             tAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTATCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTa                                                                                                   >  1:135468/1‑149 (MQ=255)
                                             tAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGGTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTa                                                                                                   <  8:137262/149‑1 (MQ=255)
                                                 gATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATg                                                                                                 >  4:297801/1‑147 (MQ=255)
                                                       gACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTACCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGCGGCTTCCTTTATGCAGTActg                                                                                         >  5:257967/1‑149 (MQ=255)
                                                       gACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGACTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTActg                                                                                         >  7:373598/1‑149 (MQ=255)
                                                       gACAGGAGTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGAGCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCGTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTATGGGGCTTCGTTTATGCTGTActg                                                                                         >  6:66493/1‑149 (MQ=255)
                                                        aCAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCgg                                                                                                                      >  2:286249/1‑119 (MQ=255)
                                                        aCAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCgg                                                                                                                      <  1:286249/119‑1 (MQ=255)
                                                        aCAGGATTATCGTCACACGTCCATTTCTTCAAGCCGATGAAGATCTTGGATTATTACCTTGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGGACGGTATACGACTTTCTGGTCCGGCGCTTAGGTGATTCCTTTATGCAGTActt                                                                                         <  2:135468/148‑2 (MQ=255)
                                                            gATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTActgcctgc                                                                                    >  7:371865/1‑149 (MQ=255)
                                                            gATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTActgcctgc                                                                                    >  4:103858/1‑149 (MQ=255)
                                                                      aCCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTAAGCAGTACTGCCTGCGACCGGaaa                                                                           <  7:119342/148‑1 (MQ=255)
                                                                      aCCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCGCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAAt                                                                          <  5:335495/149‑1 (MQ=255)
                                                                          gTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAggg                                                                      <  1:136949/149‑1 (MQ=255)
                                                                                ttCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTg                                                                >  5:357101/1‑149 (MQ=255)
                                                                                   tGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCATAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGa                                                              >  6:26642/1‑148 (MQ=255)
                                                                                    gCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACTACGTGCTGGTCCGTCGCTTAGGGGCTTCCTTTATGAAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGaaa                                                            >  7:141602/1‑149 (MQ=255)
                                                                                     cAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTg                                                                >  6:212417/1‑144 (MQ=255)
                                                                                     cAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTg                                                                <  5:212417/144‑1 (MQ=255)
                                                                                     cAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAAt                                                           >  6:183366/1‑149 (MQ=255)
                                                                                     cAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGGGGa                                                              <  3:103858/146‑1 (MQ=255)
                                                                                       aGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGATTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTg                                                         >  1:196863/1‑149 (MQ=255)
                                                                                                  aTCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGcc                                              >  4:320702/1‑149 (MQ=255)
                                                                                                   tCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCATAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGcc                                              <  5:26642/148‑1 (MQ=255)
                                                                                                     ttGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGc                                                                                                >  5:88878/1‑96 (MQ=255)
                                                                                                     ttGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGc                                                                                                <  6:88878/96‑1 (MQ=255)
                                                                                                       ggCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCta                                            >  6:132011/1‑146 (MQ=255)
                                                                                                       ggCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCta                                            <  5:132011/146‑1 (MQ=255)
                                                                                                       ggCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCt                                             <  5:183366/145‑1 (MQ=255)
                                                                                                        gCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATg                                        >  7:46915/1‑149 (MQ=255)
                                                                                                        gCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCGTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTatat                                         <  4:65611/148‑1 (MQ=255)
                                                                                                         cttcttACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGc                                       >  2:148456/1‑149 (MQ=255)
                                                                                                             ttACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTg                                    >  1:142309/1‑148 (MQ=255)
                                                                                                               aCCTGGAGATATCGAAGAAAAGTTTGCTCCCTATTTTCGCCCGGGCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGTAGGTGGAAATTGCGCCGTTCGCCTATATGCGTgg                                   <  7:167880/147‑1 (MQ=255)
                                                                                                                  tGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTACTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTa                              >  2:107649/1‑149 (MQ=255)
                                                                                                                        tatCGCAGAAAAGTTTGCTCCCCATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCtata                                          <  3:326480/131‑1 (MQ=255)
                                                                                                                        tatCGCAGAAAAGTTTGCTCCCCATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCtata                                          <  3:331809/131‑1 (MQ=255)
                                                                                                                        tatCGCAGAAAAGTTTGCTCCCCATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCtata                                          >  4:326480/1‑131 (MQ=255)
                                                                                                                        tatCGCAGAAAAGTTTGCTCCCCATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCtata                                          >  4:331809/1‑131 (MQ=255)
                                                                                                                          tCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGCCGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGa                       >  5:114683/1‑148 (MQ=255)
                                                                                                                            gCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGTTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGaaaa                    >  6:107684/1‑149 (MQ=255)
                                                                                                                             cAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAAt                   <  3:297801/149‑1 (MQ=255)
                                                                                                                              aGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTg                                                                >  7:317785/1‑103 (MQ=255)
                                                                                                                              aGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTg                                                                <  8:317785/103‑1 (MQ=255)
                                                                                                                               gAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGc                 <  1:4549/149‑1 (MQ=255)
                                                                                                                                aaaaGTTTGCGCCCTATTTTCGCCCGCTCTACGACGTGCTGGGCCGGCGCTTTGGGGCTTCCTTTATGCAGTACTGCCTTCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGCa                <  8:373598/149‑1 (MQ=255)
                                                                                                                                   aGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGCAgtc             <  3:320702/149‑1 (MQ=255)
                                                                                                                                   aGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTACTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGCAgtc             <  1:107649/149‑1 (MQ=255)
                                                                                                                                     tttGCTCCGTATTTTCGCCCTGTCTTCGACGTGCTGGTCAGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGATATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATTTGCGTGGACGTACCTTTGAAAATGCAgtcg            <  8:141602/148‑1 (MQ=255)
                                                                                                                                     tttGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTATGAAAATGCAgtcgt           >  8:322969/1‑149 (MQ=255)
                                                                                                                                       tGCTCCCTATTTTCGCCCGGTCTACGCCGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGGTCGCCTATATGCGTGGACGTACCTTTGAAAATGCAgtcgtc          <  8:289173/148‑1 (MQ=255)
                                                                                                                                         cTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGt                                     >  8:197710/1‑119 (MQ=255)
                                                                                                                                          tCCGTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTActgcct                                                                                      <  1:51337/69‑1 (MQ=255)
                                                                                                                                          tCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTActgcct                                                                                      >  2:51337/1‑69 (MQ=255)
                                                                                                                                           cccTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGCAGTCGTCATTCt     >  4:335376/1‑149 (MQ=255)
                                                                                                                                             cTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAGGTGGTAATTGCGCCGTTCGGCTAAATGCGTGGACGGACCTTTGAAAATTCAGTCGTCATTCttta  >  7:308107/1‑147 (MQ=255)
                                                                                                                                              tATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGCAGTCGTCATTCTTGa  >  4:144326/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
AAAAACCTGGATCAGCGCAGCAAAAGCGGCAGAGGCCCTGATACATAAGGATGTCGACAGGATTATCGTCACCCGTCCAGTTCTGCAAGCCGATGAAGATCTTGGCTTCTTACCTGGAGATATCGCAGAAAAGTTTGCTCCCTATTTTCGCCCGGTCTACGACGTGCTGGTCCGGCGCTTAGGGGCTTCCTTTATGCAGTACTGCCTGCGACCGGAAATAGGGAAGGTGGAAATTGCGCCGTTCGCCTATATGCGTGGACGTACCTTTGAAAATGCAGTCGTCATTCTTGA  >  NC_000913/1085525‑1085815

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