Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,406,129 T→C 100% intergenic (+144/‑61) purA → / → nsrR adenylosuccinate synthetase/nitric oxide‑sensitive repressor for NO regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,406,1290TC99.7% 1137.9 / ‑5.0 296intergenic (+144/‑61)purA/nsrRadenylosuccinate synthetase/nitric oxide‑sensitive repressor for NO regulon
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (169/125);  total (171/125)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 296 total aligned reads displayed.

GTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTTC  >  NC_000913/4405982‑4406275
                                                                                                                                                   |                                                                                                                                                  
gTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCtg                                                                                                                                                  >  1:273806/1‑150 (MQ=255)
gTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCtg                                                                                                                                                  >  1:3011600/1‑150 (MQ=255)
 tAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCtgt                                                                                                                                                 >  2:2019041/1‑150 (MQ=255)
    ttCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGc                                                                                                                                              >  2:271504/1‑150 (MQ=255)
    ttCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGc                                                                                                                                              >  1:2497390/1‑150 (MQ=255)
    ttCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGc                                                                                                                                              >  1:423402/1‑150 (MQ=255)
          tACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGttt                                                                                                                                        >  1:372872/1‑150 (MQ=255)
          tACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTa                                                                                                                                       <  2:777762/150‑1 (MQ=255)
          tACGCCTGGAATATTTTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTGGGGAACGCCATTAGGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCTAACTTAAATAAATTAGCGACACAGCCTGTGGCTGGttt                                                                                                                                        <  1:3477959/150‑1 (MQ=255)
              ccTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTTAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatca                                                                                                                                    <  1:1030611/150‑1 (MQ=255)
              ccTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTatca                                                                                                                                    >  2:1058809/1‑150 (MQ=255)
                  gCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCaa                                                                                                                                >  2:2068916/1‑150 (MQ=255)
                      atatTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAAtata                                                                                                                            <  1:3443095/150‑1 (MQ=255)
                       tatTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATaa                                                                                                                           <  1:580116/150‑1 (MQ=255)
                         ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                         <  2:1334533/150‑1 (MQ=255)
                         ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                         >  1:1435547/1‑150 (MQ=255)
                         ttttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAAt                                                                                                                         >  2:281237/1‑150 (MQ=255)
                          tttGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATg                                                                                                                        >  1:1927483/1‑150 (MQ=255)
                              cctgccGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTAtt                                                                                                                    >  1:2335317/1‑150 (MQ=255)
                               ctgccGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTAttt                                                                                                                   >  1:131790/1‑150 (MQ=255)
                                    gggCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTcc                                                                                                              <  1:684811/150‑1 (MQ=255)
                                     ggCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTccc                                                                                                             >  2:1472214/1‑150 (MQ=255)
                                     ggCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTccc                                                                                                             <  2:118507/150‑1 (MQ=255)
                                         aaCAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGAtt                                                                                                         <  2:1837392/150‑1 (MQ=255)
                                          aCAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGAttt                                                                                                        <  1:987628/150‑1 (MQ=255)
                                           cAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTc                                                                                                       <  1:87211/150‑1 (MQ=255)
                                            aGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTcc                                                                                                      >  2:2943926/1‑150 (MQ=255)
                                            aGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTcc                                                                                                      <  1:3863850/150‑1 (MQ=255)
                                            aGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGATATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTcc                                                                                                      >  2:3866775/1‑150 (MQ=255)
                                              tgtgATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCt                                                                                                    >  2:630014/1‑150 (MQ=255)
                                              tgtgATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCt                                                                                                    <  1:3306850/150‑1 (MQ=255)
                                                tgATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCttt                                                                                                  >  1:1960827/1‑150 (MQ=255)
                                                 gATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCtttt                                                                                                 <  1:2629490/150‑1 (MQ=255)
                                                 gATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCtttt                                                                                                 <  1:3528773/150‑1 (MQ=255)
                                                   tACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGa                                                                                               <  2:2031537/150‑1 (MQ=255)
                                                   tACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCGTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGa                                                                                               <  1:2684772/150‑1 (MQ=255)
                                                    aCATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAg                                                                                              <  1:1793263/150‑1 (MQ=255)
                                                       ttGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGtt                                                                                           >  2:3253344/1‑150 (MQ=255)
                                                       ttGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGtt                                                                                           >  1:334498/1‑150 (MQ=255)
                                                       ttGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGtt                                                                                           >  1:3150584/1‑150 (MQ=255)
                                                       gtGCTGGGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGtt                                                                                           <  1:771924/149‑1 (MQ=255)
                                                        tGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTg                                                                                          <  2:585774/150‑1 (MQ=255)
                                                          cTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTTTTAAATCCTCGCATTTTTCCTACCCCGAACTGAAATAAATTAGCGAAACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGAt                                                                                        <  2:727728/150‑1 (MQ=255)
                                                          cTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGAt                                                                                        >  2:3018720/1‑150 (MQ=255)
                                                          cTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGAt                                                                                        >  2:828184/1‑150 (MQ=255)
                                                               tCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCa                                                                                   <  2:3788670/150‑1 (MQ=255)
                                                                    tAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGTTGTGCAGTTaa                                                                              >  1:166632/1‑150 (MQ=255)
                                                                    tAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTaa                                                                              <  1:3277645/150‑1 (MQ=255)
                                                                     aaGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAAc                                                                             >  2:1663908/1‑150 (MQ=255)
                                                                      aGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACg                                                                            >  2:3679904/1‑150 (MQ=255)
                                                                       gCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGa                                                                           >  2:2469212/1‑150 (MQ=255)
                                                                       gCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGa                                                                           <  1:3807990/150‑1 (MQ=255)
                                                                           ttACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGttt                                                                       >  2:1622474/1‑150 (MQ=255)
                                                                           ttACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGttt                                                                       <  2:2653648/150‑1 (MQ=255)
                                                                           ttACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGttt                                                                       >  2:1622217/1‑150 (MQ=255)
                                                                           ttACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGttt                                                                       >  2:2997916/1‑150 (MQ=255)
                                                                           ttACGCTATCCGACACAGTGTTAAAAAATCGATTTTTAACTTCCCCCAACTGAAATAAATTAGCGACACAGAATGTGAATGGAATATAATCAATATAAAAGTATTTTTTCACGATTTCCCTAATGATGTTGATGTGCAGTTAAAGAGttt                                                                       >  2:2356602/1‑150 (MQ=255)
                                                                              cGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCAc                                                                    <  1:746276/150‑1 (MQ=255)
                                                                                 tATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTga                                                                 >  2:1289026/1‑150 (MQ=255)
                                                                                   tCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCGGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGGTTCACTgatt                                                               >  1:3109106/1‑150 (MQ=255)
                                                                                    ccGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAgg                                                                                             >  2:541003/1‑119 (MQ=255)
                                                                                     cGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattac                                                             <  2:1917503/150‑1 (MQ=255)
                                                                                      gACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGGTGGTTTATCATCAATTTAAATGTTTTTTTTGCGGATTTACCTTTTGTGGTTGTTGTGCGTTTTTCGAGTTTCACTgatttcg                                                            >  1:1483915/1‑150 (MQ=255)
                                                                                      gACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacg                                                            >  1:1381113/1‑150 (MQ=255)
                                                                                       acacAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacgg                                                           >  2:3355666/1‑150 (MQ=255)
                                                                                         acaGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacggat                                                         >  2:3178485/1‑150 (MQ=255)
                                                                                           aGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTACCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacggatta                                                       >  1:3626309/1‑150 (MQ=255)
                                                                                            gtgtTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTgattacggattac                                                      <  1:1263488/150‑1 (MQ=255)
                                                                                               ttAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTg                                                   >  2:1213312/1‑150 (MQ=255)
                                                                                                tAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTgc                                                  <  1:1187211/150‑1 (MQ=255)
                                                                                                   aTCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCt                                               <  2:2151385/150‑1 (MQ=255)
                                                                                                     ccTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGa                                             <  1:1876768/150‑1 (MQ=255)
                                                                                                       tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATc                                           >  1:1291058/1‑150 (MQ=255)
                                                                                                       tCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATc                                           >  2:4089848/1‑150 (MQ=255)
                                                                                                        cGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCt                                          >  2:2421014/1‑150 (MQ=255)
                                                                                                         gCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTa                                         >  1:1597784/1‑150 (MQ=255)
                                                                                                           ttttttCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACa                                       >  2:2090139/1‑150 (MQ=255)
                                                                                                            tttttCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACAt                                      >  1:2211140/1‑150 (MQ=255)
                                                                                                             ttttCCTTCCCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATg                                     >  2:3163414/1‑150 (MQ=255)
                                                                                                                 ccttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGt                                 >  2:4036368/1‑150 (MQ=255)
                                                                                                                   ttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCa                               >  1:4032751/1‑150 (MQ=255)
                                                                                                                   ttccCCGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCa                               >  2:3439496/1‑150 (MQ=255)
                                                                                                                       ccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGc                           >  2:2510146/1‑150 (MQ=255)
                                                                                                                       ccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGc                           <  1:488830/150‑1 (MQ=255)
                                                                                                                       ccGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTGTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGc                           >  1:2472009/1‑150 (MQ=255)
                                                                                                                        cGAACTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGcc                          <  2:710622/150‑1 (MQ=255)
                                                                                                                            cTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGaa                      <  2:1562081/150‑1 (MQ=255)
                                                                                                                            cTGAAATAAATTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGaa                      <  1:2068916/150‑1 (MQ=255)
                                                                                                                               aaataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAggg                   >  2:2595949/1‑150 (MQ=255)
                                                                                                                               aaataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAggg                   <  1:1003132/150‑1 (MQ=255)
                                                                                                                               aaataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCGTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAggg                   >  2:669806/1‑150 (MQ=255)
                                                                                                                                aataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGc                  <  1:3764953/150‑1 (MQ=255)
                                                                                                                                 ataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGCTTTCCCATTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCg                 <  2:1591621/150‑1 (MQ=255)
                                                                                                                                 ataaatTAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCg                 <  1:529148/150‑1 (MQ=255)
                                                                                                                                      ttAGCGACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGa            >  2:2876153/1‑150 (MQ=255)
                                                                                                                                           gACACAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGt       >  1:3905692/1‑150 (MQ=255)
                                                                                                                                            acacAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTa      >  2:225817/1‑150 (MQ=255)
                                                                                                                                            acacAGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTa      >  1:2893905/1‑150 (MQ=255)
                                                                                                                                               caGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTAttt   >  2:3637161/1‑150 (MQ=255)
                                                                                                                                                aGCCTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTAAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTTc  >  2:3567401/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTAATTCTGGTACGCCTGGCAGATATTTTGCCTGCCGGGCGAACAGTGTGATACATTGCTGTGTCGGGTAAGCCATTACGCTATCCGACACAGTGTTAAATCCTCGCTTTTTTCCTTCCCCGAACTGAAATAAATTAGCGACACAGCTTGTGGCTGGTTTATCATCAATATAAATGTATTTTTTCCCGATTTCCCTTTTGAGGTTGATGTGCAGTTAACGAGTTTCACTGATTACGGATTACGTGCGCTGATCTACATGGCGTCATTGCCAGAAGGGCGGATGACCAGTATTTC  >  NC_000913/4405982‑4406275

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: