Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,289,094 T→A 7.6% I444N (ATT→AAT)  nrfA → nitrite reductase, formate‑dependent, cytochrome

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,289,0940TA 7.6% 398.4 / 17.8 120I444N (ATT→AAT) nrfAnitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (60/50);  new base A (4/5);  total (64/56)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 122 total aligned reads displayed.

CGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTGA  >  NC_000913/4288947‑4289242
                                                                                                                                                   |                                                                                                                                                    
  gaagaaGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTgg                                                                                                                                                   >  3:672331/1‑149 (MQ=255)
    agaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                                                                                 <  4:476503/149‑1 (MQ=255)
     gaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCt                                                                                                                                                <  3:2625/149‑1 (MQ=255)
     gaagGTTTACGGATGCTCGGTACGGCGATGGAGAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATACAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGCCCATTGGCCt                                                                                                                                                >  4:756408/1‑149 (MQ=255)
         gTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctg                                                                                                                                            <  1:793479/149‑4 (MQ=255)
         gTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctg                                                                                                                                            <  2:600101/149‑4 (MQ=255)
          tttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctgg                                                                                                                                           <  4:193055/149‑5 (MQ=255)
             aCGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATgg                                                                                                                                        <  2:409279/149‑1 (MQ=255)
               ggATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGTCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGaa                                                                                                                                      <  2:737329/149‑1 (MQ=255)
                    cTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAAt                                                                                                                                 <  1:319200/149‑1 (MQ=255)
                        gTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAg                                                                                                                             <  3:550808/149‑1 (MQ=255)
                               gATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCgaga                                                                                                                      <  4:533139/149‑1 (MQ=255)
                                aTGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGaa                                                                                                                     >  3:172738/1‑149 (MQ=255)
                                   gATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCa                                                                                                                  <  1:125186/149‑1 (MQ=255)
                                    aTAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAg                                                                                                                 <  1:472209/149‑1 (MQ=255)
                                           ggcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtca                                                                                                          <  2:485660/149‑1 (MQ=255)
                                                     cgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATTTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                <  4:499268/149‑1 (MQ=255)
                                                     cgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACggc                                                                                                <  1:368699/149‑2 (MQ=255)
                                                     cgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                <  4:116620/149‑1 (MQ=255)
                                                     cgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                <  2:440950/149‑1 (MQ=255)
                                                       cacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGa                                                                                              <  4:844465/149‑1 (MQ=255)
                                                         ccAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATc                                                                                            <  4:402235/149‑1 (MQ=255)
                                                         ccAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATc                                                                                            >  1:336650/1‑149 (MQ=255)
                                                          cAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATcc                                                                                           <  4:673963/149‑1 (MQ=255)
                                                          cAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATAGCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAGAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATcc                                                                                           >  1:731410/1‑149 (MQ=255)
                                                            aaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCg                                                                                         >  3:171244/1‑149 (MQ=255)
                                                             aCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGc                                                                                        <  1:363548/149‑1 (MQ=255)
                                                                ggCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGt                                                                                     <  1:797980/149‑1 (MQ=255)
                                                                  cgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                                                                                   <  1:175745/149‑1 (MQ=255)
                                                                  cgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                                                                                   <  4:199234/149‑1 (MQ=255)
                                                                  cgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                                                                                   >  4:604356/1‑149 (MQ=255)
                                                                   gcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTggg                                                                                  <  2:446160/149‑1 (MQ=255)
                                                                   gcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTggg                                                                                  <  4:584730/149‑1 (MQ=255)
                                                                     gcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACGTCATCAAAACGGTGATCCCGCAGTgggca                                                                                <  3:756408/149‑3 (MQ=255)
                                                                        tgctgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACCAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGgaagaa                                                                             <  1:813665/149‑1 (MQ=255)
                                                                           tgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAg                                                                          >  4:775726/1‑149 (MQ=255)
                                                                                  cAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGCAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTa                                                                   >  1:174980/1‑149 (MQ=255)
                                                                                   aaaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaa                                                                  >  1:895198/1‑149 (MQ=255)
                                                                                    aaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaa                                                                 >  4:161344/1‑149 (MQ=255)
                                                                                    aaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaa                                                                 >  4:841222/1‑149 (MQ=255)
                                                                                      ggCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaaaa                                                               >  1:398789/1‑149 (MQ=255)
                                                                                         aTCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGTTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACgg                                                            >  3:872088/1‑149 (MQ=255)
                                                                                         aTCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACgg                                                            >  2:95522/1‑149 (MQ=255)
                                                                                         aTCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACgg                                                            >  2:575149/1‑149 (MQ=255)
                                                                                          tCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGt                                                           >  3:211488/1‑149 (MQ=255)
                                                                                              ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCa                                                                                                                  >  1:13449/1‑90 (MQ=255)
                                                                                              ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAgg                                                                         <  1:809365/131‑1 (MQ=255)
                                                                                              ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGt                                                       >  2:565468/1‑149 (MQ=255)
                                                                                              ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGt                                                       >  1:490172/1‑149 (MQ=255)
                                                                                              ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGACCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAgg                                                                         >  2:809365/1‑131 (MQ=255)
                                                                                               cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                      >  3:701726/1‑149 (MQ=255)
                                                                                               cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                      >  2:838717/1‑149 (MQ=255)
                                                                                               cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                      >  2:283505/1‑149 (MQ=255)
                                                                                                aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                     >  4:855710/1‑149 (MQ=255)
                                                                                                 tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                >  1:153452/1‑105 (MQ=255)
                                                                                                 tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                                <  2:153452/105‑1 (MQ=255)
                                                                                                 tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                    >  3:737243/1‑149 (MQ=255)
                                                                                                 tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                    <  4:212551/149‑1 (MQ=255)
                                                                                                 tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                    <  4:737243/149‑1 (MQ=255)
                                                                                                 tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                    >  3:125458/1‑149 (MQ=255)
                                                                                                 tGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                    >  2:75452/1‑149 (MQ=255)
                                                                                                      tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCaa                                               <  1:283505/149‑1 (MQ=255)
                                                                                                      tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCaa                                               >  2:466425/1‑149 (MQ=255)
                                                                                                       ccAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAAt                                              <  4:332777/149‑1 (MQ=255)
                                                                                                         aGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGGAAAAACGGTCTGTTAAGCCAATaa                                            >  3:177971/1‑149 (MQ=255)
                                                                                                                ggATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGtt                                     >  2:363045/1‑149 (MQ=255)
                                                                                                                 gATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTc                                    <  1:723497/149‑1 (MQ=255)
                                                                                                                      ctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCt                               >  3:770147/1‑149 (MQ=255)
                                                                                                                       tcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGgaa                                                                                <  3:856313/99‑1 (MQ=255)
                                                                                                                       tcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGgaa                                                                                >  4:856313/1‑99 (MQ=255)
                                                                                                                       tcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTc                              >  1:732393/1‑149 (MQ=255)
                                                                                                                         aaCTAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGc                            >  3:229689/1‑149 (MQ=255)
                                                                                                                            cAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAg                         >  3:741358/1‑149 (MQ=255)
                                                                                                                            cAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAg                         >  1:438473/1‑149 (MQ=255)
                                                                                                                             aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGt                                      >  3:3315/1‑135 (MQ=255)
                                                                                                                             aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGt                                      <  4:3315/135‑1 (MQ=255)
                                                                                                                             aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAgg                        >  4:108182/1‑149 (MQ=255)
                                                                                                                             aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAgg                        >  1:693735/1‑149 (MQ=255)
                                                                                                                              aaGAGAAAGCCCAGCAGGCCATTGGGCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTCATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGTCAATAACCCCGTTCCGCCACGCAAggg                       >  2:406518/1‑149 (MQ=255)
                                                                                                                              aaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAggg                       <  2:168683/149‑1 (MQ=255)
                                                                                                                              aaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAggg                       >  4:51146/1‑149 (MQ=255)
                                                                                                                                gagaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGc                     >  1:197983/1‑149 (MQ=255)
                                                                                                                                gagaAAGCCCAGCAGGCCACTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGGGAGCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAAGAACCCCGTTCCGCCTCGCAAGGGGc                     <  2:829082/149‑1 (MQ=255)
                                                                                                                                 agaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCg                    >  3:116525/1‑149 (MQ=255)
                                                                                                                                   aaaGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGa                  >  3:665572/1‑149 (MQ=255)
                                                                                                                                   aaaGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGa                  >  4:365089/1‑149 (MQ=255)
                                                                                                                                        tgttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAaca             >  2:712960/5‑149 (MQ=255)
                                                                                                                                          ttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGt                                                           >  2:619803/3‑101 (MQ=255)
                                                                                                                                          ttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGt                                                           <  1:619803/99‑1 (MQ=255)
                                                                                                                                          ttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGt                                                           >  1:269193/3‑101 (MQ=255)
                                                                                                                                          agcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTc                                                          >  2:344814/1‑102 (MQ=255)
                                                                                                                                          agcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTc                                                          <  1:344814/102‑1 (MQ=255)
                                                                                                                                          agcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAAcaca           <  1:557250/149‑1 (MQ=255)
                                                                                                                                           gcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACaa          <  4:125458/149‑1 (MQ=255)
                                                                                                                                           gcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACaa          <  3:266388/149‑1 (MQ=255)
                                                                                                                                             agGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATg        >  2:204850/1‑149 (MQ=255)
                                                                                                                                             agGCCAATGGCCTGCACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATg        >  3:587335/1‑149 (MQ=255)
                                                                                                                                                 cATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGt    >  2:668517/1‑149 (MQ=255)
                                                                                                                                                 cATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGt    <  3:51146/149‑1 (MQ=255)
                                                                                                                                                   ttGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTGa  >  2:26150/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTGA  >  NC_000913/4288947‑4289242

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