Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,076,800 C→T A695T (GCG→ACG)  putA ← fused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,076,8000CT98.5% 207.6 / NA 67A695T (GCG→ACG) putAfused DNA‑binding transcriptional repressor/proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase PutA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (29/37);  minor base A (1/0);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATACACACCACAGGCCCTAATGGACGGTGGGTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGCGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCGAGATCCAGCCCTGCCGAGTTGTCGCGCCCGTG  >  NC_000913/1076559‑1077043
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
atanccaccacAGGCCCTAATGGACGGTGGGTTTCGTGAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGGTGCTGCATCTGGCTTGCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAAc                                                                                                                                                                                                                                            <  2:593831/246‑1 (MQ=255)
  acacacCACAGGCCCTAATGGACGGTGGGTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCg                                                                                                                                                                                                                                          <  1:196051/251‑1 (MQ=255)
      accacAGGCCCTAATGGACGGTGGGTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACt                                                                                                                                                                                                                                      >  1:482922/1‑251 (MQ=255)
        cacaGGCCCTAATGGACGGTGGGTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACttt                                                                                                                                                                                                                                    >  2:628382/1‑251 (MQ=255)
                cTAATGGACGGTGGGTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGcc                                                                                                                                                                                                                            >  1:64756/1‑251 (MQ=255)
                    tGGACGGTGGCTTACGTTAGCGAAATCATCCCGCACATGGCTGGCGTAGGAGGGGCGAAAATAGACCGCATCGCGCACGTCGGCGATGGAGTTATTGAAGGTTTTTCCGGCCTCACGCACCAGAATACAAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCCAAATCGATGCGTGTTAATCCGGTGGTGTGGCACATCAGATTGTTGTGTTATTAAACGCACTTTCCAggtactgtt                                                                                                                                                                                                                        >  2:42570/1‑243 (MQ=255)
                            gggTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTc                                                                                                                                                                                                                 <  1:494847/250‑1 (MQ=255)
                            gggTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCa                                                                                                                                                                                                                <  1:584480/251‑1 (MQ=255)
                                 tCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACgcgg                                                                                                                                                                                                           >  2:617233/1‑251 (MQ=255)
                                     tAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTg                                                                                                                                                                                                       >  1:347372/1‑251 (MQ=255)
                                          aaaTCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCTCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTc                                                                                                                                                                                                  <  1:145921/251‑1 (MQ=255)
                                              cagcccGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAAGCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACg                                                                                                                                                                                               <  2:68294/247‑1 (MQ=255)
                                               atcCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGc                                                                                                                                                                                              >  2:249387/1‑250 (MQ=255)
                                                   cgcgccTGTCCGGTGTAGTAGTGGAGAAAATCTACCACTTCGCGCACTGCGGCAATGGCGTTACGGAATGTTGTTCCGGCCTCACGAGCCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCTCGGCAGCGCGGTGCAAAATCGCTGCGCTTTCAGCCGGAGTCGTTGTAAATCAGATTGGCGTGTTATTAATCGCACATGCCAGCGCCTGTTCTACTGCGCGCGGTGTTGCTTCACGCACATa                                                                                                                                                                                         <  2:85054/247‑1 (MQ=255)
                                                     cACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGc                                                                                                                                                                                       <  1:116087/251‑1 (MQ=255)
                                                                   agtagtGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGGTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGAAAACCAGATTGGCGTGTTATTACCCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCttt                                                                                                                                                                          <  2:18627/250‑1 (MQ=255)
                                                                    gtagtGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCg                                                                                                                                                                        <  1:628382/251‑1 (MQ=255)
                                                                    ggagtGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCg                                                                                                                                                                        <  2:344598/249‑1 (MQ=255)
                                                                         gtagaAAATCGCCCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTACCCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCGCGCACATAGCCCACAATATCTTTCGGTTcc                                                                                                                                                                   <  2:352318/249‑1 (MQ=255)
                                                                         ggAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTcc                                                                                                                                                                   <  1:526875/251‑1 (MQ=255)
                                                                           agaAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCATGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGc                                                                                                                                                                                       <  1:41362/229‑1 (MQ=255)
                                                                           agaAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCATGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGc                                                                                                                                                                                       >  2:41362/1‑229 (MQ=255)
                                                                               aaTCGACCTCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAggg                                                                                                                                                             <  1:305016/251‑1 (MQ=255)
                                                                                  cGACCGCTTCGCGCACTTCGGCAATTGCGTTACTGAAGGTTTTTCCGGCCGCACGCACCAGAATGCCAATCAGTTGCTGCATCTGGCTTACCATCAGCACGGCAGCGCGGTGCATAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTTGCTTCACGCACATAGCCCTCAATATCTTTCGGTTCCGCAGGGTta                                                                                                                                                          <  2:114564/251‑1 (MQ=255)
                                                                                  cGACCGCTTCGCGCACGTCGGCAATGGCGTTACTGAAGGTTTTGCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTta                                                                                                                                                          <  2:144389/251‑1 (MQ=255)
                                                                                      cGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTtaata                                                                                                                                                       >  2:235294/1‑250 (MQ=255)
                                                                                         ttCGCGCACTTCGGCAATGGCGTCACTGAAGGTTTTCCCAGCCTCAAGCACCAGATTACCAATCAGTTGCTGCATATGGCTTTTCATCAGCACGGAAGTGCGGTGTAAAACAGCTGCGCGTTCAGCCAGAGGAGTAGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCACACAATATCTTTCGGTTCCGCAGGGTTAACAACgg                                                                                                                                                   >  2:385286/1‑251 (MQ=255)
                                                                                                    cagCAATGGCGTTACTGAAGGTTTTTCCTGCCTCACGCACCAGAATACCAATTAGTTTCTGCATCTGGCTTGCCATCAGCACGGCAGCGCGGTGCAAAATCGCTCCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCCCTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCCCGCACATATCCCACAATATCTATCGGTTTCTCAGGGTTAATAACGGGCCACAtctc                                                                                                                                         <  2:1966/248‑1 (MQ=255)
                                                                                                     ggCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCAc                                                                                                                                       <  1:320255/251‑1 (MQ=255)
                                                                                                     ggCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCAc                                                                                                                                       <  1:443739/251‑1 (MQ=255)
                                                                                                      gcgaTGGCGTTCCTGAAGGTTTGTCCGGCCTCACGCACCAGAATGCCAATCAGTTGCTGCATCTTGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCAcc                                                                                                                                      <  2:564508/248‑1 (MQ=255)
                                                                                                        aaTGGCGTTACTGAAGGTTTTACCGGCCTCACGAACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGAAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCAGTTTTACTTCAAGCGGAGTGGCTTAACGCACATAGGCCACAATATCTTTAGGTTCCGCAGGGTTAATAACggatgacaacacacctg                                                                                                                                    >  1:303370/1‑236 (MQ=255)
                                                                                                          tGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTACAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGcc                                                                                                                                  >  2:519269/1‑251 (MQ=255)
                                                                                                           ggCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGAAATCTCACCTGCCg                                                                                                                                 >  1:319116/1‑251 (MQ=255)
                                                                                                                 ccTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAACCCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCCCCTGCCGCTACCg                                                                                                                           <  2:1451/250‑1 (MQ=255)
                                                                                                                    gacggTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTCCCATCAGCCCGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCCCTTTCCAGCGCCTGTTCTACTTCCCGCGGCGTGGCTTCACGCACATAGCCCACAATCTCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGgtt                                                                                                                        <  2:545832/248‑1 (MQ=255)
                                                                                                                     aaGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGgttg                                                                                                                       >  1:177979/1‑251 (MQ=255)
                                                                                                                     aaGGTTTTTCCGGCCTCACGCACCAGAATAACAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGAGCACAATCGCTGCGCGATCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCAAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGCTAATAACGGGCGACATCTCACCTGCCGTTACCGgttg                                                                                                                       >  2:542493/1‑251 (MQ=255)
                                                                                                                       ggTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGgttgtt                                                                                                                     <  1:627713/251‑1 (MQ=255)
                                                                                                                          ttttCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTcc                                                                                                                   <  1:555793/250‑1 (MQ=255)
                                                                                                                                    tCACGCACCAGAATACAAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTCCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACTCGGCGTGGCTTCACGCACATAGACCACACTATCTTTCGGGTCCGCAGGGTTAATAACGGACGACATATCACCTGCCGCGACCGGTTGTTCCAGCATTgcca                                                                                                         >  2:207710/1‑247 (MQ=255)
                                                                                                                                         cACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCt                                                                                                   >  1:413204/1‑251 (MQ=255)
                                                                                                                                         cACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGGTGTTCCAGCATTGGCAAGGCCt                                                                                                   >  1:529693/1‑251 (MQ=255)
                                                                                                                                          aCCAGAATACCGATCAGTTGATGCATCTGGCTTTCCATCCGCACGGCAGCGCGGTGCACAATCGCTGCGCGTTCAGCCGGAGGCCTGACAAACCAGATTAGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATACCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGTCATCTCACCTGCCGCTACCGGTTGTTCCAGCATTCGCCAGGCCTg                                                                                                  >  2:313590/1‑251 (MQ=255)
                                                                                                                                            cAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCAGCAGGGTTAATAACGGGAGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCGAGGTCTga                                                                                                 >  2:309689/1‑249 (MQ=255)
                                                                                                                                                aTACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCAttt                                                                                            >  1:186642/1‑251 (MQ=255)
                                                                                                                                                   ccAATCAGTTGCTGCATCTGGCTTTTCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGTCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGGTCTACTGCACTCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCGGCAGGGTTAATAACGGGCGACATCTCAcc                                                                                                                                      <  2:145862/206‑1 (MQ=255)
                                                                                                                                                   ccAATCAGTTGCTGCATCTGGCTTTCCATAAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGAGTTATTAACCGCACTTTCCAGCGCCTGTTCTAATTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCAcc                                                                                                                                      >  1:145862/1‑206 (MQ=255)
                                                                                                                                                      aTCAGTTGCTGCATCTGGTTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCGGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAg                                                                                      <  2:466654/251‑1 (MQ=255)
                                                                                                                                                         aGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGc                                                                                   >  1:332770/1‑251 (MQ=255)
                                                                                                                                                                      ggCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAgg                                                                      <  1:68681/251‑1 (MQ=255)
                                                                                                                                                                      ggCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAg                                                                       <  1:292078/250‑1 (MQ=255)
                                                                                                                                                                       gCTTTTCATCATCACGGCAGCGCGGTGCAAAATCGCTGCGCTGTTAGCCGGAGGGGTGGCAAACCAGATGGGCGTGTTAATAACCGCACTTTCCAGCGCCTTTTTTCCTTCACCCGGCGTGGCTTCACGCACATAGCCCACAATATAGTGCGGTTCCGCGGGGTTAAGAACGGGCCACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATGGAGCAggg                                                                     <  2:29796/251‑1 (MQ=255)
                                                                                                                                                                                  gCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCCCTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGc                                                          <  1:178329/251‑1 (MQ=255)
                                                                                                                                                                                         aGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCCCTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGccc                                                   >  2:138660/1‑250 (MQ=255)
                                                                                                                                                                                                  cAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCTTGTTATTACCCGCACTTTCCAGCGCGTGTTCGACTGCACGCGGCGAGGCTTCACGCACATAGCCCGCAATTTCTTTCGGTTCCGCAGGGTTAAAAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTATTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGtt                                          <  1:385286/251‑1 (MQ=255)
                                                                                                                                                                                                         gatGCGCGTTCAGCCGGAGGCGTGGGAAACCAGATTGGCGTGTTATTAACCTCACTTTCCAGCGTCTGTTCTGCGTCACGCGGCGTGGTTTCACGCACACACCCCCCAATATCTTTCTGTTCCGCAGGGTTAATAACGGGGCACTGCGCCTCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCTAGGCGGTGTTCGTTAGc                                   <  1:453005/249‑1 (MQ=255)
                                                                                                                                                                                                         gCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATTTCTTTCGGTTCCGCAGGTTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATGTTTGCAGTGCACTATTGAGCAGGGCAGAGGATAGCGAGGCCAGGCGGTGTTCGTTAGc                                   <  2:64756/251‑1 (MQ=255)
                                                                                                                                                                                                            gcgcGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCg                                  >  1:555558/1‑249 (MQ=255)
                                                                                                                                                                                                             cgcgTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCg                                  >  1:319351/1‑248 (MQ=255)
                                                                                                                                                                                                                   tCAGCCGGAGGCGTGGCAACCTCGATGGGCGTGTTATTACCCGCGCTTCCCAGCGCCTTTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTGCGTTAGCGAGATCCa                           <  2:547931/250‑1 (MQ=255)
                                                                                                                                                                                                                     aGCCGGAGGCGTGGCAAACCAGATTGGTGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCCCGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCGAGATCCAGcc                        <  2:332770/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCCCTTTCCAGCGACTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCGAGATCCAGCCCTGCCg                  >  1:377919/1‑251 (MQ=255)
                                                                                                                                                                                                                           aGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCAACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCGAGATCCAGCCCTGCCg                  <  1:155494/251‑1 (MQ=255)
                                                                                                                                                                                                                           aGGCGTGGCAAACCAGATTGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGACCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGTGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCGGGGACCAGCCCTGCCg                  >  2:36417/1‑251 (MQ=255)
                                                                                                                                                                                                                                tGGCAAACCAGATTGGCGTGCTATTAACCGCACTTTCTACAGCCTGGTGAACTTCACGAGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCCCAGGGTTAATAACGGGCGACATCTCCCCTGCCCCCACCGGTTGTTCCAGCATTGGCAAGACCTGCCATTTTTGCAGTGAACTATTAGGCAGGGCAGCGGCGAGCGACGCCCGGCGGTGTTCGTTAGCGCGACCCAGCGCTGCCGAGTTg             >  2:660025/1‑251 (MQ=255)
                                                                                                                                                                                                                                            ttGGCGTGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCGAGATCCAGCCCTGCCGAGTTGTCGCGCCCggg  >  1:586667/1‑248 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
ATACACACCACAGGCCCTAATGGACGGTGGGTTTCGTTAGCGAAATCATCCCGCACCTGTCCGGCGTAGTAGTGGAGAAAATCGACCGCTTCGCGCACTTCGGCAATGGCGTTACTGAAGGTTTTTCCGGCCTCACGCACCAGAATACCAATCAGTTGCTGCATCTGGCTTTCCATCAGCACGGCAGCGCGGTGCAAAATCGCTGCGCGTTCAGCCGGAGGCGTGGCAAACCAGATTGGCGCGTTATTAACCGCACTTTCCAGCGCCTGTTCTACTTCACGCGGCGTGGCTTCACGCACATAGCCCACAATATCTTTCGGTTCCGCAGGGTTAATAACGGGCGACATCTCACCTGCCGCTACCGGTTGTTCCAGCATTGGCAAGGCCTGCCATTTTTGCAGTGCACTATTGAGCAGGGCAGAGGAGAGCGAGGCCAGGCGGTGTTCGTTAGCGAGATCCAGCCCTGCCGAGTTGTCGCGCCCGTG  >  NC_000913/1076559‑1077043

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