Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,287,009 A→T 7.3% intergenic (+483/+57) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,0090AT 7.3% 329.6 / 7.6 116intergenic (+483/+57)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base A (56/48);  new base T (2/8);  total (59/57)
Fisher's exact test for biased strand distribution p-value = 5.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCG  >  NC_000913/1286860‑1287154
                                                                                                                                                     |                                                                                                                                                 
ttttgatttCCGGGGCTTTTGACAGTGTTGGTGGTGGGGTAAGGATTACTCCGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                                                                                   <  1:1838472/143‑1 (MQ=14)
 tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                                                                                  <  2:966282/150‑1 (MQ=255)
     aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCCGCAACGCACACACCCTGACGCACCACCCGAACCTTAGTCGAACCTTCTCAACCGTACCCCCAAGGGc                                                                                                                                              >  2:2876899/1‑150 (MQ=14)
     aaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCCGCAACGCACGCTCGCTTACCCACCACTCGGACCTTAGTACAAGCATCACAACCATACCCCCAAGGGc                                                                                                                                              >  2:2877067/1‑150 (MQ=2)
              gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGTAGCTTCTCATCCTTCCCCGCATGGGCAGAATAttt                                                                                                                                     <  1:2395314/150‑1 (MQ=25)
                ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                                                                                   <  1:660983/150‑1 (MQ=38)
                 tttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTCCGCCCATCGCCCAACGGCTCGATGCCTTCCGCAACGCACGCACCCAGACGCCCCAACCGAACCTTACTCGAAGCTACACATCAAACCCCCCATGGGCACACTAATTGAt                                                                                                                                  >  2:2269309/1‑150 (MQ=2)
                     gaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                                                              <  1:652249/149‑1 (MQ=39)
                     gaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                                                              <  1:2521080/149‑1 (MQ=39)
                     aaaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCg                                                                                                                              >  2:921930/1‑150 (MQ=39)
                      aaGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCgg                                                                                                                             >  2:1991381/1‑150 (MQ=21)
                       aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                                                            <  2:203453/150‑1 (MQ=39)
                       aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                                                            <  2:1280814/150‑1 (MQ=39)
                       aGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGa                                                                                                                            <  2:2827083/150‑1 (MQ=39)
                                     gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                                                                              <  1:74595/150‑1 (MQ=255)
                                     gggAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCttccttaaa                                                                                                              <  2:39241/150‑8 (MQ=255)
                                       gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                                                                                            <  1:2895801/150‑1 (MQ=255)
                                        caggaggaCTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCGCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCCTCCTTCCCCGCATGGGCAGAATATTGGATTGCGGATTCACTTGAGAGTTCAg                                                                                                           <  2:1804679/143‑1 (MQ=14)
                                        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                                                                                           <  2:1408480/150‑1 (MQ=255)
                                        aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCACTTGAGAGTTCAg                                                                                                           <  2:2445377/150‑1 (MQ=255)
                                         aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                                          <  1:921930/150‑1 (MQ=255)
                                         aGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAgg                                                                                                          <  1:2934863/150‑1 (MQ=255)
                                            aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCGTAGTCGAAGCTTCTCAACCGTCCCCGCAGGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGAc                                                                                                       >  2:2430836/1‑150 (MQ=11)
                                                cTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACtttt                                                                                                   <  1:204855/150‑1 (MQ=255)
                                                  cAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCTTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTACCCGCATGGGCAGAATATTTGATTGCGGATTCGGGTGTGAGTTCAGGGGCTTTTGa                                                                                                 >  1:2581725/1‑150 (MQ=14)
                                                    gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                                                                               <  1:1755570/150‑1 (MQ=255)
                                                     cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                                                              <  1:2447624/150‑1 (MQ=255)
                                                     cgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAg                                                                                              <  1:381777/150‑1 (MQ=255)
                                                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                                                             <  2:1093348/150‑1 (MQ=255)
                                                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGTATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                                                             <  2:2804471/150‑1 (MQ=255)
                                                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                                                             <  2:669993/150‑1 (MQ=255)
                                                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                                                             <  2:2683693/150‑1 (MQ=255)
                                                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                                                             <  1:1310088/150‑1 (MQ=255)
                                                      gcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGt                                                                                             <  1:1844997/150‑1 (MQ=255)
                                                        tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCCAGGCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGA                                                                                           <  2:434012‑M1/150‑1 (MQ=255)
                                                          cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACGTTTGAAAGTGATG                                                                                         >  2:954068‑M1/1‑150 (MQ=255)
                                                          cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATG                                                                                         <  2:924569‑M1/150‑1 (MQ=255)
                                                                  gCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGG                                                                                 <  2:2882484‑M1/150‑1 (MQ=255)
                                                                   cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCAGGGGCAGAATAGTTGATTGCGGATTCGCTTGTGCGGTCAGGGACGTTTGAAAGTGAGGGGGGGGGGG                                                                                >  2:786186‑M1/1‑150 (MQ=255)
                                                                   cccTTCGGGCCGTTGCCCGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGG                                                                                <  1:3059601‑M1/150‑1 (MQ=255)
                                                                       tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                                                            <  1:3107364‑M1/150‑1 (MQ=255)
                                                                       tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                                                            <  1:353352‑M1/150‑1 (MQ=255)
                                                                       tCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAG                                                                            <  1:1408839‑M1/150‑1 (MQ=255)
                                                                            tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa                                                                       <  2:1946636/150‑1 (MQ=255)
                                                                              gTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt                                                                     <  1:2738383/150‑1 (MQ=255)
                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTccac                                                                    >  2:3099761‑M1/1‑146 (MQ=255)
                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGATTGTACAGGGACGGTGGAGAGTGGTGGTGGTGGGGGAAGGATTACTc                                                                    >  2:366704/1‑150 (MQ=17)
                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                    >  2:1335094/1‑150 (MQ=255)
                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                    >  2:3099799/1‑150 (MQ=255)
                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                    >  2:565340/1‑150 (MQ=255)
                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTGTGGAAAGTGATGGTGGTGGGGGAAGGATTACTc                                                                    >  2:366719/1‑150 (MQ=255)
                                                                               ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTACTGGGAGTGTAGAGAGTGATCGTGGTGGGGGAAGGATTACTc                                                                    >  2:3115643/1‑150 (MQ=17)
                                                                                tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGACGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa                                                                   <  1:1602239/150‑1 (MQ=255)
                                                                                   cTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCATAGTCGAAGCTTCTCAACCATCCCCCCATGGGCCGAAAATTTGATTTCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTaaacacag                                                                >  1:2262688/1‑143 (MQ=17)
                                                                                     gCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCt                                                              >  2:552586/1‑150 (MQ=255)
                                                                                              ctctctCGCTGTCGCTCTAGCCGAACCTTCGTCGAAGCGTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                                                     <  1:2033284/150‑1 (MQ=14)
                                                                                              ctctctCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                                                     <  1:2836715/150‑1 (MQ=255)
                                                                                              ctctctCGCTGACGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTAGATTGCGGATTCGCCTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                                                     <  1:2033288/150‑1 (MQ=255)
                                                                                              ctctctCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                                                     >  2:215812/1‑150 (MQ=255)
                                                                                              ctctctCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTc                                                     >  2:215865/1‑150 (MQ=255)
                                                                                                     gCTGTCGCTCGAGACGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                              <  2:607040/150‑1 (MQ=255)
                                                                                                     gCTGCCGCTCGAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                              <  2:2981451/150‑1 (MQ=255)
                                                                                                     gCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTc                                              >  1:2892706/1‑150 (MQ=255)
                                                                                                        gACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTCATTGCGGATTCGCTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTagacccatcgaagtcgatgacggaacaga                                           >  2:2321633‑M1/1‑121 (MQ=255)
                                                                                                            cTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGt                                       >  1:985496/1‑150 (MQ=255)
                                                                                                            cTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGt                                       >  1:2339413/1‑150 (MQ=255)
                                                                                                             tccagtcGAACCTGGGTCGAAGCTTCTCCTCCTTCCCCGCTTGGGCAGAATAGTTGATTGCGGATTCGCGTGAGAGTTCAGGGCCTTTTGAAAGTGATGGGGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                      <  1:1091491/144‑1 (MQ=1)
                                                                                                             tCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                      >  1:968967/1‑150 (MQ=255)
                                                                                                             tCGAGCCGAACCGTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtt                                      <  1:1091563/150‑1 (MQ=255)
                                                                                                             ncGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCCGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGGGGTGGGGGACGGATTACGCAGCGCTGCGCGCTTCGCCCTTCGGGTCGtg                                      >  1:8986/2‑149 (MQ=255)
                                                                                                               aagtcGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGc                                    >  1:2623403/5‑150 (MQ=9)
                                                                                                                  tcGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                 >  1:332040/2‑150 (MQ=18)
                                                                                                                   cGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGTTCGtgtccggc                                >  2:2620153/1‑143 (MQ=255)
                                                                                                                     aaCCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTTCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                              <  2:2101091/150‑1 (MQ=255)
                                                                                                                        cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCaa                           >  1:3116881/1‑150 (MQ=255)
                                                                                                                        cTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATAGTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTTCGCGCATCGCCCTTCGGGTCGTTGCCAGCGGCaa                           >  1:1308980/1‑150 (MQ=255)
                                                                                                                         ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaacgtctgagacgatgtcgggaggtgtacattctgctcccaat                          >  2:2770971‑M1/1‑104 (MQ=255)
                                                                                                                         ttAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCACCGCTGCGCGCATCGCCCTTCGGGTCGTGGCCGGAGGCAAc                          >  2:535665/1‑150 (MQ=14)
                                                                                                                         ttAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTcgaaccttcgaagtcgatgacggcagatttacagtctgctcccttt                          <  1:684830‑M1/150‑47 (MQ=255)
                                                                                                                             tCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTACGGGTCGTTGCCTGCGGCAACGctc                      >  2:1795850/1‑150 (MQ=255)
                                                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGGTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTAGTAAGCGCTTAGCTCATCACCCTTCGGGTCGTTGCCTGGGGCAAAGttctc                    >  2:2722919/1‑150 (MQ=14)
                                                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctctc                    <  1:565340/150‑1 (MQ=255)
                                                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGctctc                    >  2:1807976/1‑150 (MQ=255)
                                                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTGTCCTGGGGCAACGctctc                    >  2:2241539/1‑150 (MQ=255)
                                                                                                                               gAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTGGGGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCCCTTCGCCCGTCGGGTCGTGGCCGCCGGCAACGccctc                    >  2:2561965/1‑150 (MQ=14)
                                                                                                                                    ttCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt               <  2:332040/150‑1 (MQ=37)
                                                                                                                                    gtctcCTCCGTCCCCGCCTGGGCAGAATATTTGCTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGACAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCGGCGCGCGTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCt               <  2:152968/149‑1 (MQ=14)
                                                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGATTCCGCCGGCAAAGCACTCTCGCtc              >  1:3009775/1‑149 (MQ=14)
                                                                                                                                     tctcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCATCGCCCTTCGGGTCGTTGCCTGCGGCAACGCACTCTCGCTg              >  1:3008749/1‑150 (MQ=255)
                                                                                                                                       tcATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCGCTCGCGGgc            >  2:661951/1‑150 (MQ=255)
                                                                                                                                          tCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCGCTCGCTGgcgcg         >  1:906184/1‑149 (MQ=255)
                                                                                                                                           ccTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGATTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTc        <  2:2099636/150‑1 (MQ=255)
                                                                                                                                             ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGa      >  1:2420946/1‑150 (MQ=255)
                                                                                                                                             ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCCACCCTCACTCGCGGGCGCTCGa      >  1:2420948/1‑150 (MQ=14)
                                                                                                                                             ttCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCACTCGCTGGCGCTCGa      >  2:1459179/1‑150 (MQ=255)
                                                                                                                                               ccccGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGt    <  2:2953414/150‑1 (MQ=12)
                                                                                                                                                 ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGACGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCGCACTCGCAGGCGCACGAGTCg  >  2:424194/1‑150 (MQ=14)
                                                                                                                                                 ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCg  >  2:2525981/1‑150 (MQ=255)
                                                                                                                                                 ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCg  >  2:2756665/1‑150 (MQ=255)
                                                                                                                                                 ccGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCACGAGGCg  >  2:2172272/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCG  >  NC_000913/1286860‑1287154

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