Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,123,930 A→T V176E (GTA→GAA)  cytR ← DNA‑binding transcriptional repressor CytR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,123,9300AT98.7% 266.0 / ‑4.7 75coding (527/1026 nt)cytRDNA‑binding transcriptional repressor CytR
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (26/48);  minor base . (1/0);  total (27/48)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGA  >  NC_000913/4123683‑4124172
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
gAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGCAGGTGAAGTCGCCACGGGCGATGCATTGCGGATCACCCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCCTCTCTTCGGGGCCGGCTATACCGCCAATCCGTTTTTGCCCTTTCTCATATAAATAATTTTCTg                                                                                                                                                                                                                                                  >  1:540239/1‑251 (MQ=255)
   cgcCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCAGCGCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCAt                                                                                                                                                                                                                                               <  1:504085/251‑1 (MQ=255)
         cacaTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGGGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCAGCAAATg                                                                                                                                                                                                                                         <  2:453640/251‑1 (MQ=255)
          acaTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgc                                                                                                                                                                                                                                        <  1:237680/251‑1 (MQ=255)
          acaTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgc                                                                                                                                                                                                                                        <  2:318007/251‑1 (MQ=255)
           caTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgcg                                                                                                                                                                                                                                       <  1:85028/251‑1 (MQ=255)
            aTCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAACGGCCTCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgcgg                                                                                                                                                                                                                                      >  1:94461/1‑251 (MQ=255)
             tCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATAAAGCCAATCCGTTTTTGCCCTTGCTCCTATAAATAATTTTCTTCATCAAATtcggc                                                                                                                                                                                                                                     >  1:474071/1‑251 (MQ=255)
              cGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGACCAAGCAGCTGCGGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCGTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATgcggc                                                                                                                                                                                                                                     <  2:249644/250‑1 (MQ=255)
               gCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGCCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGcggcgg                                                                                                                                                                                                                                   <  2:364867/251‑1 (MQ=255)
                  aTGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCGTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCATCCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTGTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGt                                                                                                                                                                                                                                  <  2:515516/249‑1 (MQ=255)
                      cAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGact                                                                                                                                                                                                                            >  1:493571/1‑249 (MQ=255)
                      cAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCGTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGAtt                                                                                                                                                                                                                            <  1:523074/251‑1 (MQ=255)
                      cAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCGTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCAGAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCGCTTCGTGACCGGCTAGACAGCCAATCCGGGTATGCCCTTGCGCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGAtt                                                                                                                                                                                                                            <  2:522665/251‑1 (MQ=255)
                      cAGAAGACAGCAGTAGGCGGTTGTGGAAGAGCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGAtt                                                                                                                                                                                                                            <  2:130670/251‑1 (MQ=255)
                          agacaggcgTAGGCAGGTGTGGAAGATCAGGCCGCTGCTGGCTCGCTTTGCTTCCGGCTGCGAAGGTGAAGTGGCCACGGGCGACGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCCGTATTGACCCAGCGGAATCTCTTCCGGACCGACGATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATTATTTTCTACATCAAATGCGGCGGGCAGATTGTCg                                                                                                                                                                                                                        <  2:521908/243‑1 (MQ=255)
                            aCAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCCTTGCTTCCCGCTTCGAAGGCGAAGTCGCCACGGGCGATGTATTGCGGATAAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGGGACACAGCGGCATCTCTTCGGGACCGGCTATACAACCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGcggcgggccgaatgttga                                                                                                                                                                                                                       >  1:258417/1‑238 (MQ=255)
                             caggtgTAGGCGGTGGGGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTCCCGGCTGCGAAGGGGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGTCTTGCAGGCGGTATTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATAGAAATAATTTTCTGCATCAAATGCGGCTTTCATATCGTCGata                                                                                                                                                                                                                    <  2:381858/246‑1 (MQ=255)
                              agcagTAGGCAGTTGAGGAAGATCAAGCAACTGCTGCCTCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCAAAGGGCGCTGTATTGCGGATCAACCATAATGCCGCAGCGCCGAAGCGGCTGAACATAGCCTCGCAGGCGGTTGTGTTACAGCTGCATTTCTTAGGGATCGGCTCTACAGCCAATCCCCTTATGACCTTGCCCATAAACAAAAATTTCTGCATCACATGcggcgggccgatttttggta                                                                                                                                                                                                                     >  2:240368/1‑236 (MQ=255)
                                    aGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGGTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGGCAGGATGTCGGTATggaatg                                                                                                                                                                                                             >  1:231740/1‑246 (MQ=255)
                                         gTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGc                                                                                                                                                                                                        <  1:351176/251‑1 (MQ=255)
                                         gTTGTGGAAGAGCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCTATGGATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATAATAACTGTAGGc                                                                                                                                                                                                        <  2:368263/251‑1 (MQ=255)
                                               gAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGcagctcc                                                                                                                                                                                                  >  1:122090/1‑251 (MQ=255)
                                                         cagcTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACCCAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCCTCaa                                                                                                                                                                                                                                            >  2:426328/1‑200 (MQ=255)
                                                         cagcTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCaa                                                                                                                                                                                                                                            <  1:426328/200‑1 (MQ=255)
                                                                     cGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCATCTCCAGCTCCGGTGCAAAttcgttcg                                                                                                                                                                            <  1:445907/251‑1 (MQ=255)
                                                                                  gCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGGATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATcgg                                                                                                                                                               <  1:124879/251‑1 (MQ=255)
                                                                                    ttCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATCATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCCTCTCTTCGGGACCGGCTATACCGCCAATCCGTTTATGACCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGGTTTTTGATATGAACTGTAAGGAGCTCCCGCTCCGGGGCGAATTCGGTCGCCATCAACATcgggg                                                                                                                                                             >  1:226432/1‑249 (MQ=255)
                                                                                     tCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATcggc                                                                                                                                                              <  2:493571/249‑1 (MQ=255)
                                                                                     tCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATcggc                                                                                                                                                              <  1:122846/249‑1 (MQ=255)
                                                                                       gAAGGTGAAGACGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCCGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGAAAAAGCGGCATCTCATCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGGCAGATTGTTGATTATAACAGTAGGCAGCTCCAGCCCCGGTGCAAAtttgttttcccacaccctccgtgggg                                                                                                                                                          >  2:519020/1‑231 (MQ=255)
                                                                                        aaGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAg                                                                                                                                                         <  1:19080/251‑1 (MQ=255)
                                                                                          ggTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGAt                                                                                                                                                       <  1:337284/251‑1 (MQ=255)
                                                                                             gAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCGGAACATAGCCTTGCAGGCGGTAGGGACACAGCGGCATCTCTGCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCAGCAAATGCGGCGGTCAGATTGGCGATATGAACTGGAGGCAGCGCCTGCTCCGGTGCAAATTCGTTCGCCAGCTCCATCGGCGGCAGATTAc                                                                                                                                                    <  2:542235/251‑1 (MQ=255)
                                                                                                gTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACCGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCCAATGCGGCGGTCAGATTGTTGATATGAACTGTAGGCCGCTCCCGGTCCGGGGCAAATTTGTTCGCCATCACCATCGGCGGCAGAATACGcc                                                                                                                                                 >  1:296843/1‑250 (MQ=255)
                                                                                                          gCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAAt                                                                                                                                       <  2:231740/251‑1 (MQ=255)
                                                                                                             aTGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTCGTGACACAGCGGCCTCTCTTCGGGACCGGCCATACCGCCAATCCGTTTATGCCCTTGCTCCTATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGGTATTAAATGTAGGCTGCTCCCGCCCCCGGGCGAATTTGTTCGCCCTCACCATccgcgggcggattccgtgttcccgaactgg                                                                                                                                    >  2:74237/1‑227 (MQ=255)
                                                                                                             aTGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTcc                                                                                                                                            <  1:74237/243‑1 (MQ=255)
                                                                                                              tGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGGGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGAATGGCGATATTAACTGTAGGCAGCTCCAGGTCCGGGGCAAATTTGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTTCTCAATgacg                                                                                                                                   >  1:118677/1‑248 (MQ=255)
                                                                                                                    gCGGATCAACCATAATGCCGCAGCGACGCCGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCCAATGCGGCGGTTAAATTGTCGATATGAAATGTAGGCCGCTCCAGCTCCGGTGCCAATTCCTTCGCCCTCACCCTCGGCGGCCGATTTCCCTGTTTCTCCATGCTGGcctcc                                                                                                                             >  1:163670/1‑247 (MQ=255)
                                                                                                                     cGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCaa                                                                                                                            <  2:122090/251‑1 (MQ=255)
                                                                                                                     cGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCaa                                                                                                                            <  2:18246/251‑1 (MQ=255)
                                                                                                                          cAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGc                                                                                                                       <  2:163670/251‑1 (MQ=255)
                                                                                                                             ccATAATGCCGCTGCGACGCGGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGc                                                                                                                    <  2:491957/251‑1 (MQ=255)
                                                                                                                              cATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTGCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGcc                                                                                                                   <  2:189959/251‑1 (MQ=255)
                                                                                                                               aTAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCt                                                                                                                                                 <  1:273922/220‑1 (MQ=255)
                                                                                                                                 aaTGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCGCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCtt                                                                                                                 <  1:341002/250‑1 (MQ=255)
                                                                                                                                                  gcgcCTGAACATAGCCTTGCAGGCGGTAGAGACACCGCGGCATCTTTTCGAGACCGGCTATACAGCCAATCCGTTTATGCCCCTGCTCCTATAAATAATTTTCTGCATCAAATGCGGCGGTCAGGTTGTCGATTTGAAATCTGGGCCGCTCCAGCTCCGGTGCCAATTCGGTCGCCCTCACCAACCGCAGCACATTCCGCCGTTCCTCAAAGCAGGCCTCAACAGGCAGCCCTGAACCCCGACACAACATg                                                                                               >  1:531481/1‑251 (MQ=255)
                                                                                                                                                   cgcCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCAGCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCTTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGc                                                                                              <  2:540239/251‑1 (MQ=255)
                                                                                                                                                    gcCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGcc                                                                                             <  2:127759/251‑1 (MQ=255)
                                                                                                                                                    gcCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATATCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGGCCGAATTTAGATATGAACTGTTAGCAGCTCCAACCCCGGGGGAAAATCGTTCGCCATCACCAACGGCGGCAGATTTAGCTGTTTCTCAATGCCGGCCTTAAAAGGCCGACCTGAAACCAGCAacaacaatcc                                                                                             >  1:98694/1‑247 (MQ=255)
                                                                                                                                                     ccTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAAttc                                                                                                                                                                                 <  1:138706/166‑1 (MQ=255)
                                                                                                                                                        gAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCa                                                                                         <  1:240773/251‑1 (MQ=255)
                                                                                                                                                        gAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCAGCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCGCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCGCAAGGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCa                                                                                         <  1:100041/251‑1 (MQ=255)
                                                                                                                                                                    gCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGt                                                                              <  1:197613/250‑1 (MQ=255)
                                                                                                                                                                       ggCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCACTCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCGAATGCGGCGGTCAGATTGTTCATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGACTTGCACCCAGCAACAAAATGCCATCAATTTGCTTGGtgttg                                                                          >  1:243095/1‑248 (MQ=255)
                                                                                                                                                                             aGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCaa                                                                     <  1:551738/250‑1 (MQ=255)
                                                                                                                                                                                  cacaGCGGCATCTCCTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGAt                                                               <  2:327057/251‑1 (MQ=255)
                                                                                                                                                                                      gCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAg                                                           <  2:94461/251‑1 (MQ=255)
                                                                                                                                                                                               ctTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGtttttt                                                    <  2:441203/249‑1 (MQ=255)
                                                                                                                                                                                                   gggACCGGCTATACAGCCAATCCGTTTATTCCCTTGCTCAAATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGAGATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCctgct                                              <  2:448252/251‑1 (MQ=255)
                                                                                                                                                                                                    gtaCCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCctgctg                                             <  2:118677/249‑1 (MQ=255)
                                                                                                                                                                                                      aCCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCCGCCTTGAACCCAGCAACAAAATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGAt                                           >  2:23307/1‑251 (MQ=255)
                                                                                                                                                                                                          gCTATACAGCCAATCCGTTTATGCCCTTGCTCATATACATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAAATGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGCCAGACTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGGTTTTTCCTGCTGAttttt                                       >  1:390265/1‑250 (MQ=255)
                                                                                                                                                                                                             attcAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCGGTTCCGCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGAt                                     <  2:437867/247‑1 (MQ=255)
                                                                                                                                                                                                                 aGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATgcgc                                 >  1:322824/1‑250 (MQ=255)
                                                                                                                                                                                                                          gTTTATGCCCTTGCTCATATAAATAATTTTCTGCATCAACTGCGGCGCTCAAATTGTCGATATGACCTGTTGGCAGCTCCAGCTCCGGTGAAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCCTCAATTTGCTTGGAGATGATCAAATCGATAAAGGTTTTTTCCTGCTGGTTTTGATGCGCACAGTCGcc                        >  1:500366/1‑250 (MQ=255)
                                                                                                                                                                                                                               tGCCCTTGCTCATATAACTAATTTTCTGCATCAAATGAGGCGGTCAGACTTTCGATATGAACTGTAGGCAGCTCCACCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCCGATTAACCTGTTCCTCCATGCTGGCATCAAACGGCAGCCTTGAAACCAGCCACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGGTTTTTTCTGCTGCTTTTTATTCGCCCCGTTGCCAATccg                  >  2:542073/1‑249 (MQ=255)
                                                                                                                                                                                                                                  ccTTGCTCATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACCACCTGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCAc               >  1:10985/1‑251 (MQ=255)
                                                                                                                                                                                                                                     tGCTCATATAAATAATTTTCTGCATCAAAGGCGGCGGTCAGATTGTCGATATGAACCGTAGGCAGCTCCAGCTCCGGTGCAAATGCGTTCGCCATCACCATCGGCGGCGGATTACGCGGTTCCTCAATGCTGGCATCAAACGTCTGCCTTGAACCCAGCAACAACATGCCATGAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTGTGCGTGCTGATGTTGATGCTGACAGTCCCCAATCAGCACCa             <  2:41801/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cATATAAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGatat        <  1:68738/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tataAATAATTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATccgaacccgctatcc      >  1:322770/1‑237 (MQ=255)
                                                                                                                                                                                                                                              aaataaTTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTTCAGCTCCGGTGCAAATTAGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATTAAACGGCAACCTTGAACCCAGAAAAAACATGCCCTCAATTTGCTTGGTGATGATCAAATCGTAAAAGGTTTTTTCCTGCTGATTTTTATGCGCACAGTCGCCAATCAGCACCATATATCCGTg   >  1:322576/1‑251 (MQ=255)
                                                                                                                                                                                                                                              aaataaTTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCCCCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCCACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGGTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTg   >  1:507986/1‑251 (MQ=255)
                                                                                                                                                                                                                                                aataatTTTCTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGGTTTTTCCTGCTGATTTTGATGTTCACAGTCGCCAATCAGCACCAGATATCCGTGa  <  2:560504/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
GAGCGCCATCACATCGCTATGGCAGAAGACAGCAGTAGGCGGTTGTGGAAGATCAAGCAGCTGCTGCATCGCTTTGCTTCCGGCTTCGAAGGTGAAGTCGCCACGGGCGATGTATTGCGGATCAACCATAATGCCGCAGCGACGCAGCGCCTGAACATAGCCTTGCAGGCGGTAGTGACACAGCGGCATCTCTTCGGGACCGGCTATACAGCCAATCCGTTTATGCCCTTGCTCATATAAATAATTTACTGCATCAAATGCGGCGGTCAGATTGTCGATATGAACTGTAGGCAGCTCCAGCTCCGGTGCAAATTCGTTCGCCATCACCATCGGCGGCAGATTACGCTGTTCCTCAATGCTGGCATCAAACGGCAGCCTTGAACCCAGCAACAACATGCCATCAATTTGCTTGGTGATGATCAAATCGATAAAGGTTTTTTCCTGCTGATTTTGATGCGCACAGTCGCCAATCAGCACCAGATATCCGTGA  >  NC_000913/4123683‑4124172

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