Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,550,578 C→G 5.3% intergenic (‑33/+63) yfeX ← / ← yfeY porphyrinogen oxidase, cytoplasmic/RpoE‑regulated lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,550,5780CG 5.3% 260.1 / 5.0 75intergenic (‑33/+63)yfeX/yfeYporphyrinogen oxidase, cytoplasmic/RpoE‑regulated lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (38/33);  new base G (2/2);  total (40/35)
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.63e-01

GCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTC  >  NC_000913/2550433‑2550724
                                                                                                                                                 |                                                                                                                                                  
gCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGt                                                                                                                                                 >  7:26772/1‑149 (MQ=255)
  aaaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCCATCCGTgg                                                                                                                                               >  3:70422/1‑149 (MQ=255)
  aaaTGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTgg                                                                                                                                               >  4:164268/1‑149 (MQ=255)
     tGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAt                                                                                                                                            <  4:201236/149‑1 (MQ=255)
       ttttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCAttg                                                                                                                                          >  8:52161/1‑149 (MQ=255)
        tttACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGGATTGCGTGGCAtta                                                                                                                                          <  1:150153/148‑2 (MQ=255)
            cTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCATATCCCCGCCCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGATACATTGTTCCTCCTGAAAATAATAATTCCACGCAATCCGTGTCATTGTTGa                                                                                                                                      <  6:25417/148‑1 (MQ=255)
                 cgcacgcaGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTa                                                                                                                                >  4:138124/1‑149 (MQ=255)
                 cgcacgcaGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTa                                                                                                                                <  3:118863/149‑1 (MQ=255)
                    acgcaGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTAttt                                                                                                                             >  4:195619/1‑149 (MQ=255)
                       caGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTAc                                                                                                                          <  3:53899/149‑1 (MQ=255)
                       caGGGCGTCAACTTCCCCTTTCACGTGGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGAAAAATGCCCCTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCCCGCAATCCGTGGCATTGTTGAGCGCTATTTTAc                                                                                                                          <  4:202255/149‑1 (MQ=255)
                            cGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGtt                                                                                                                     <  3:262477/149‑1 (MQ=255)
                               cAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGttttt                                                                                                                  >  8:336817/1‑149 (MQ=255)
                                 aCTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTc                                                                                                                <  1:297890/149‑1 (MQ=255)
                                   ttCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTc                                                                                                                >  3:118469/1‑147 (MQ=255)
                                   ttCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTc                                                                                                                <  4:118469/147‑1 (MQ=255)
                                        cTTTCACGTCGGCTTCGATCCAAATCGCCGCTCGGCAATGTTCTGGCAAACTGCCACTCTGAACCTGACACATTGATCCTCCTGAAAATATTAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATATTACGTGTTTTTTCgctgctg                                                                                                         <  2:77338/149‑3 (MQ=255)
                                           tCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGcggcggcgg                                                                                                      >  1:225721/1‑149 (MQ=255)
                                             aCGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                                                                     <  5:326009/148‑1 (MQ=255)
                                             aCGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGt                                                                                                     <  3:138124/148‑1 (MQ=255)
                                               gTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTg                                                                                                  >  4:148024/1‑149 (MQ=255)
                                               gTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTg                                                                                                  >  2:230686/1‑149 (MQ=255)
                                                 tGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACGTTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGtt                                                                                                <  3:184109/149‑1 (MQ=255)
                                                   gCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGtt                                                                                                >  1:13887/1‑147 (MQ=255)
                                                   gCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTc                                                                                               <  2:178762/148‑1 (MQ=255)
                                                          tCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCa                                                                                        >  1:343456/1‑148 (MQ=255)
                                                             aaaTCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAt                                                                                     >  6:132638/1‑148 (MQ=255)
                                                               aTCGCCGCGCGCCAATGATCTGGCAAACTGCCACTCAGCACCTGCGACATTGTTCCTCCTGAAAATAATACTCCCCCCCAAGCCGGCGCTTTGTTGAGCGCTATGTACCCGGCCTTTCCGCCGCGGCGGTCGGTTCCCACGCCAATTaa                                                                                  <  4:70384/149‑1 (MQ=255)
                                                                         ggCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAAcgccgc                                                                            >  6:335137/1‑145 (MQ=255)
                                                                         ggCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGa                                                                        <  3:164268/149‑1 (MQ=255)
                                                                            aaTGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAATTAATAATGCCACGAAAGCCGTGGCATTGTTGAGCGGTATTTTACGTGTTTTTTCCCGCCGCCGTTTTGTACAGAGGCAAATTACCGCCCCAAGaaaa                                                                     >  4:135449/1‑149 (MQ=255)
                                                                            aaTGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATaa                                                                     >  2:219204/1‑149 (MQ=255)
                                                                                   tGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCGCGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAAGTAACGCCGCCAGATAATTTTACt                                                              >  4:296642/1‑149 (MQ=255)
                                                                                     gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGtt                                                                                                >  6:51879/1‑113 (MQ=255)
                                                                                     gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGtt                                                                                                <  5:51879/113‑1 (MQ=255)
                                                                                     gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACt                                                              <  8:26772/147‑1 (MQ=255)
                                                                                     gCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGa                                                            <  1:230686/149‑1 (MQ=255)
                                                                                             ccACTCTGCACCTGAGCCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTATACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCa                                                     <  4:207694/148‑1 (MQ=255)
                                                                                                 tctGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGtttt                                                >  8:146660/1‑149 (MQ=255)
                                                                                                   tGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGtttt                                                <  7:89084/147‑1 (MQ=255)
                                                                                                     gaCCTGAGTAATTGTTACTGCTGTAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCa                                             <  3:135449/147‑1 (MQ=255)
                                                                                                     aaCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGttttt                                               <  7:146660/146‑1 (MQ=255)
                                                                                                     aaCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCa                                             <  5:132638/148‑1 (MQ=255)
                                                                                                     aaCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAg                                            >  8:96055/1‑149 (MQ=255)
                                                                                                     aaCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAg                                            >  5:314783/1‑149 (MQ=255)
                                                                                                      aCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAgg                                           <  3:195619/149‑1 (MQ=255)
                                                                                                      aCCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAgg                                           <  5:335137/149‑1 (MQ=255)
                                                                                                         tGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGtt                                                  >  4:153853/1‑139 (MQ=255)
                                                                                                         tGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGtt                                                  <  3:153853/139‑1 (MQ=255)
                                                                                                         tGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTa                                        >  6:105563/1‑149 (MQ=255)
                                                                                                          gagaCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTAt                                       >  4:222653/1‑149 (MQ=255)
                                                                                                           agaCATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATc                                      <  5:105563/149‑1 (MQ=255)
                                                                                                               aTTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTTACGCCGCCAGATAATTTTACTGACTATCCAGTTTTTCAGGTTATCGTcc                                  >  6:304685/1‑149 (MQ=255)
                                                                                                                      ctcctGAAAATAATAATGCCACGCAATCCTTTGCATTTTTTATCGCTATTTTACGTGTTTTTTCGCGGCGTCGTTTTGTTCATACGCAAATTAACGCCTCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGc                            >  2:223869/1‑148 (MQ=255)
                                                                                                                      ctcctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGc                            <  7:96055/148‑1 (MQ=255)
                                                                                                                       tcctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCAt                          >  8:230535/1‑149 (MQ=255)
                                                                                                                         ctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGTGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAgg                                           >  5:179349/1‑130 (MQ=255)
                                                                                                                         ctGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGTGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAgg                                           <  6:179349/130‑1 (MQ=255)
                                                                                                                           gAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTagc                      >  3:38094/1‑147 (MQ=255)
                                                                                                                            aaaaTAATAATGCCACGCAATCCGTGTCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAc                      <  2:22410/148‑1 (MQ=255)
                                                                                                                            aaaaTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAc                      <  4:233594/148‑1 (MQ=255)
                                                                                                                            aaaaTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAc                      <  3:222653/148‑1 (MQ=255)
                                                                                                                             aaaTAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAAccc                    >  4:133202/1‑149 (MQ=255)
                                                                                                                               tcaacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCtt                  >  7:296612/6‑149 (MQ=255)
                                                                                                                               tcaacaatGCCACGGATTGCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCtt                  >  3:190573/6‑149 (MQ=255)
                                                                                                                                 aataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCAtt                  >  2:292893/1‑147 (MQ=255)
                                                                                                                                  ataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAg               <  1:219204/149‑1 (MQ=255)
                                                                                                                                  ataatGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAg               >  7:328698/1‑149 (MQ=255)
                                                                                                                                      tGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTa                                                                <  2:194540/96‑1 (MQ=255)
                                                                                                                                      tGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTa                                                                >  1:194540/1‑96 (MQ=255)
                                                                                                                                        ccACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCAc         >  6:293414/1‑149 (MQ=255)
                                                                                                                                         cACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACg        >  7:170599/1‑149 (MQ=255)
                                                                                                                                              aaTCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCAtt   >  6:284189/1‑149 (MQ=255)
                                                                                                                                               aTCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTc  >  4:312430/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCAAATGTTTTACTGGCCGCACGCAGGGCGTCAACTTCCCCTTTCACGTTGGCTTCGATCCAAATCGCCGCGCGGCAATGTTCTGGCAAAATGCCACTCTGAACCTGAGACATTGTTCCTCCTGAAAATAATAATGCCACGCAATCCGTGGCATTGTTGAGCGCTATTTTACGTGTTTTTTCGCGGCGGCGGTTTGTTCAGACGCAAATTAACGCCGCCAGATAATTTTACTGACTTTCCAGTTTTTCAGGGTATCGTCCGAAGGCATTAACCCTTCAGGACCACGCCATTC  >  NC_000913/2550433‑2550724

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