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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,289,0940TA 6.3% 420.9 / 10.5 134I444N (ATT→AAT) nrfAcytochrome c552 nitrite reductase
Reads supporting (aligned to +/- strand):  ref base T (64/61);  new base A (4/5);  total (68/66)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

CCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTG  >  NC_000913/4288946‑4289241
                                                                                                                                                    |                                                                                                                                                   
 cGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGCCCATTg                                                                                                                                                   >  1:28963/1‑150 (MQ=255)
  ggAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTgg                                                                                                                                                  >  2:98616/1‑150 (MQ=255)
  ggAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTgg                                                                                                                                                  >  2:796917/1‑150 (MQ=255)
  ggAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTgg                                                                                                                                                  >  2:547332/1‑150 (MQ=255)
  ggAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTgg                                                                                                                                                  <  2:455598/150‑1 (MQ=255)
   gaagaaGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGc                                                                                                                                                 >  2:479340/1‑150 (MQ=255)
    aagaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                                                                                <  2:187073/150‑1 (MQ=255)
     agaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCt                                                                                                                                               >  2:309411/1‑150 (MQ=255)
      gaagGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTg                                                                                                                                              >  2:862351/1‑150 (MQ=255)
         ggTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAAc                                                                                                                                           <  2:544183/150‑1 (MQ=255)
          gTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctgg                                                                                                                                          <  1:65235/150‑5 (MQ=255)
          gTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctgg                                                                                                                                          <  1:1408143/150‑5 (MQ=255)
           tttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctggg                                                                                                                                         <  2:1041488/150‑6 (MQ=255)
            ttACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCAATGGCCTgctgggc                                                                                                                                        <  1:205763/150‑7 (MQ=255)
               cGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGaa                                                                                                                                     >  2:42113/1‑150 (MQ=255)
                ggATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCAGCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAc                                                                                                                                    <  1:541513/150‑1 (MQ=255)
                ggATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAc                                                                                                                                    <  1:1027772/150‑1 (MQ=255)
                ggATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAc                                                                                                                                    <  2:573246/150‑1 (MQ=255)
                ggATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAAc                                                                                                                                    <  1:1291998/150‑1 (MQ=255)
                 gATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACa                                                                                                                                   >  2:455996/1‑150 (MQ=255)
                    gCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAAt                                                                                                                                >  2:481803/1‑150 (MQ=255)
                    gCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAAt                                                                                                                                >  1:1354101/1‑150 (MQ=255)
                     cTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATc                                                                                                                               >  1:249296/1‑150 (MQ=255)
                     cTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATc                                                                                                                               >  1:844288/1‑150 (MQ=255)
                       cGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCaa                                                                                                                             <  2:229515/150‑1 (MQ=255)
                        ggTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAg                                                                                                                            >  2:439443/1‑150 (MQ=255)
                         gTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAgg                                                                                                                           <  2:434550/150‑1 (MQ=255)
                            cGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCg                                                                                                                        <  1:769303/150‑1 (MQ=255)
                             ggCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCga                                                                                                                       >  1:31958/1‑150 (MQ=255)
                                  tGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGc                                                                                                                  <  2:1124090/150‑1 (MQ=255)
                                    gATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAg                                                                                                                >  1:195515/1‑150 (MQ=255)
                                      tAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCAAGAAGCAGGa                                                                                                              >  1:485419/1‑150 (MQ=255)
                                        aaGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACt                                                                                                            >  2:982784/1‑150 (MQ=255)
                                             gcggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcatc                                                                                                       <  2:283343/150‑1 (MQ=255)
                                              cggATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcatca                                                                                                      <  1:451951/150‑1 (MQ=255)
                                                    cacgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACgg                                                                                                >  1:1069030/1‑150 (MQ=255)
                                                    cacgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACgg                                                                                                >  1:1407964/1‑150 (MQ=255)
                                                     acgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGt                                                                                               <  1:439443/150‑1 (MQ=255)
                                                      cgcacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCCACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTg                                                                                              >  1:746999/1‑150 (MQ=255)
                                                        cacCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGAt                                                                                            >  1:798076/1‑150 (MQ=255)
                                                          ccAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATcc                                                                                          >  2:1188043/1‑150 (MQ=255)
                                                           cAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGAt                                                                                            >  1:640692/1‑147 (MQ=255)
                                                           cAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGAt                                                                                            <  2:640692/147‑1 (MQ=255)
                                                            aaaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCg                                                                                        >  1:952597/1‑150 (MQ=255)
                                                             aaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGc                                                                                       <  2:79368/150‑1 (MQ=255)
                                                             aaCTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGc                                                                                       <  1:1348977/150‑1 (MQ=255)
                                                                 ggCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcatca                                                                                                      >  1:587212/1‑131 (MQ=255)
                                                                 ggCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTtcatca                                                                                                      <  2:587212/131‑1 (MQ=255)
                                                                  gcgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                                                                                  <  1:454212/150‑1 (MQ=255)
                                                                  gcgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                                                                                  <  2:622366/150‑1 (MQ=255)
                                                                  gcgcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTgg                                                                                  >  2:901944/1‑150 (MQ=255)
                                                                    gcgcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGga                                                                                <  2:923517/150‑1 (MQ=255)
                                                                      gcCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACAAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGGCGAGAAGCAGGACTTCATCATAACGGTTATCCCGCAGTGGgaag                                                                              <  2:315703/150‑1 (MQ=255)
                                                                           ctgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGcc                                                                                                                                                >  2:132274/1‑79 (MQ=255)
                                                                            tgGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAgg                                                                        >  1:1250807/1‑150 (MQ=255)
                                                                                  ccAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTa                                                                  <  1:1362337/150‑1 (MQ=255)
                                                                                    aaaGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTaaa                                                                >  2:126985/1‑150 (MQ=255)
                                                                                         cATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACgg                                                           >  1:869523/1‑150 (MQ=255)
                                                                                         cATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACgg                                                           <  2:860188/150‑1 (MQ=255)
                                                                                          aTCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGt                                                          >  2:502268/1‑150 (MQ=255)
                                                                                             aCCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTg                                                       >  2:1339403/1‑150 (MQ=255)
                                                                                              cccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGt                                                      >  2:515757/1‑150 (MQ=255)
                                                                                              cacATGAAAACCAAAACCCTCATATATCAACCAAAGAGAAACCACAACATTAAAATACCCTGAAAATGGAACAAATAAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGt                                                      >  2:1079085/3‑150 (MQ=255)
                                                                                               ccATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGtt                                                     <  1:409883/150‑1 (MQ=255)
                                                                                                cATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTa                                                    >  1:957615/1‑150 (MQ=255)
                                                                                                 aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                   >  1:66187/1‑150 (MQ=255)
                                                                                                 aTGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTaa                                                   <  1:862351/150‑1 (MQ=255)
                                                                                                    aaaTCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGcc                                                <  1:807771/150‑1 (MQ=255)
                                                                                                    aaaTCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGcc                                                <  1:548179/150‑1 (MQ=255)
                                                                                                       tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAAt                                             >  1:173333/1‑150 (MQ=255)
                                                                                                       tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAAt                                             <  2:605952/150‑1 (MQ=255)
                                                                                                       tCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAAt                                             >  2:167837/1‑150 (MQ=255)
                                                                                                        ccAGATCCCGGATATCTCATCCAAAGAGAAAGCCCAGGAGTCGATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCATGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATa                                            <  2:541585/150‑1 (MQ=255)
                                                                                                         cAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGc                                                                       <  1:127217/122‑1 (MQ=255)
                                                                                                            aTCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAAccc                                        >  1:603363/1‑150 (MQ=255)
                                                                                                               ccGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAGCAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGt                                     >  1:113098/1‑150 (MQ=255)
                                                                                                               ccGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGt                                     >  1:181085/1‑150 (MQ=255)
                                                                                                                      tctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCt                              >  1:1156770/1‑150 (MQ=255)
                                                                                                                       ctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTc                             >  2:82846/1‑150 (MQ=255)
                                                                                                                       ctcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTc                             >  1:285361/1‑150 (MQ=255)
                                                                                                                        tcAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCg                            >  2:1173599/1‑150 (MQ=255)
                                                                                                                         cAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGc                           >  2:880167/1‑150 (MQ=255)
                                                                                                                              aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAggg                      >  1:286838/1‑150 (MQ=255)
                                                                                                                              aaaGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAggg                      <  2:181085/150‑1 (MQ=255)
                                                                                                                                  agaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAGCAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCgg                  <  2:113098/150‑1 (MQ=255)
                                                                                                                                  agaAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCgg                  <  2:1069030/150‑1 (MQ=255)
                                                                                                                                       gCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAac             >  2:1105667/1‑150 (MQ=255)
                                                                                                                                        cccAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAaca            <  2:66187/150‑1 (MQ=255)
                                                                                                                                        atgttcagGCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAaca            >  1:153063/6‑150 (MQ=255)
                                                                                                                                         ccAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAacac           <  1:1173599/150‑1 (MQ=255)
                                                                                                                                            gcagGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAAt        <  2:249296/150‑1 (MQ=255)
                                                                                                                                              agGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATgg      <  1:982784/150‑1 (MQ=255)
                                                                                                                                              agGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATgg      <  1:509207/150‑1 (MQ=255)
                                                                                                                                              agGCCATTGGCCTGAACATGGAACAAATCAAGGCCAAGAAGCAGGACTTCATCAAAACGGTGATCCCGGAGTTGTAAGAACAGGCACGTAAAAACTGTCTGTTAAGCCAATAACCCCGTTCCTCCTCTCAAGGGGCGGAAAACTCAATgg      <  2:485419/150‑1 (MQ=255)
                                                                                                                                               ggCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGa     >  1:962869/1‑150 (MQ=255)
                                                                                                                                               ggCCAATGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGCAAGAACAGCCACGTAAAAACGGACTGTTAAGCCAAAAACCCCGTTCCGCCTCGCAAGGGCCGCAAAACACAATgca     >  1:406883/1‑148 (MQ=255)
                                                                                                                                                gCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAg    <  1:1095644/150‑1 (MQ=255)
                                                                                                                                                 ccATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGt   >  2:521700/1‑150 (MQ=255)
                                                                                                                                                 ccATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGt   <  2:31958/150‑1 (MQ=255)
                                                                                                                                                  cATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTg  <  2:869523/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCGGAAGAAGGTTTACGGATGCTCGGTACGGCGATGGATAAAGCGGCGGATGCACGCACCAAACTGGCGCGCCTGCTGGCGACCAAAGGCATCACCCATGAAATCCAGATCCCGGATATCTCAACCAAAGAGAAAGCCCAGCAGGCCATTGGCCTGAACATGGAACAAATCAAGGCCGAGAAGCAGGACTTCATCAAAACGGTGATCCCGCAGTGGGAAGAACAGGCACGTAAAAACGGTCTGTTAAGCCAATAACCCCGTTCCGCCTCGCAAGGGGCGGAAAACACAATGGAGTG  >  NC_000913/4288946‑4289241

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