Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,335,684 G→A 7.0% intergenic (+226/‑147) yciX → / → acnA uncharacterized protein/aconitate hydratase 1; aconitase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,335,6840GA 7.0% 279.9 / 8.5 86intergenic (+226/‑147)yciX/acnAuncharacterized protein/aconitate hydratase 1; aconitase A
Reads supporting (aligned to +/- strand):  ref base G (39/41);  new base A (4/2);  total (43/43)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTAT  >  NC_000913/1335536‑1335832
                                                                                                                                                    |                                                                                                                                                    
ctcAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcg                                                                                                                                                      >  2:561018/1‑149 (MQ=255)
   aaCACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTgcggcg                                                                                                                                                   >  1:366262/1‑149 (MQ=255)
     cacaAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCggg                                                                                                                                                 <  2:198150/149‑1 (MQ=255)
     cacaAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCggg                                                                                                                                                 <  3:666314/149‑1 (MQ=255)
     cacaAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCggg                                                                                                                                                 <  2:81925/149‑1 (MQ=255)
       caAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGt                                                                                                                                               >  3:323688/1‑149 (MQ=255)
        aaaaTAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGtt                                                                                                                                              >  2:292106/1‑149 (MQ=255)
          aataatAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGt                                                                                                                                            >  4:53570/1‑149 (MQ=255)
             aataaCCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAt                                                                                                                                         <  3:51833/149‑1 (MQ=255)
                    gTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGTTTGTCAtttgcttt                                                                                                                                  <  3:54891/149‑1 (MQ=255)
                    gTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgcttt                                                                                                                                  >  3:132206/1‑149 (MQ=255)
                    gTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCCGCAGCGGGGtcaa                                                                                                                                            >  4:422318/1‑136 (MQ=255)
                    gTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCCGCAGCGGGGtcaa                                                                                                                                            <  3:422318/139‑4 (MQ=255)
                      tAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCAtttgctttgc                                                                                                                                <  4:546168/149‑1 (MQ=255)
                            tgtAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAg                                                                                                                          <  4:595264/149‑1 (MQ=255)
                                   tATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTtctc                                                                                                                   >  4:144181/1‑149 (MQ=255)
                                   tATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTtctc                                                                                                                   <  2:366262/149‑1 (MQ=255)
                                   tATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTtctc                                                                                                                   >  4:527235/1‑149 (MQ=255)
                                   tATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTtctc                                                                                                                   >  4:133782/1‑149 (MQ=255)
                                    aTGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTctcc                                                                                                                  <  3:736927/149‑1 (MQ=255)
                                       aTCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCctct                                                                                                               >  1:581346/1‑149 (MQ=255)
                                             acaAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATc                                                                                                         >  3:265276/1‑149 (MQ=255)
                                               aaaTATGACATATTGGTCCTGGATGGGCGCGTTTTCTCTGACGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCaa                                                                                                       >  2:346635/1‑149 (MQ=255)
                                               aaaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCaa                                                                                                       >  1:715589/1‑149 (MQ=255)
                                                aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTAt                                                                                                          >  3:83246/1‑145 (MQ=255)
                                                aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTAt                                                                                                          <  4:83246/145‑1 (MQ=255)
                                                aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAt                                                                                                      <  4:415986/149‑1 (MQ=255)
                                                aaTATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAt                                                                                                      >  3:415986/1‑149 (MQ=255)
                                                 atatGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAAtt                                                                                                     <  4:269874/149‑1 (MQ=255)
                                                    tGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTgg                                                                                                  <  3:712535/149‑1 (MQ=255)
                                                    tGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTgg                                                                                                  <  3:253386/149‑1 (MQ=255)
                                                       aaTATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgtt                                                                                               <  3:449881/149‑1 (MQ=255)
                                                          atTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGgttgtt                                                                                            <  2:178117/149‑1 (MQ=255)
                                                            tGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTAt                                                                                          <  2:715589/149‑1 (MQ=255)
                                                              gTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCa                                                                                        >  4:448960/1‑149 (MQ=255)
                                                                 cTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAAt                                                                                     >  1:346790/1‑149 (MQ=255)
                                                                 cTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAAt                                                                                     >  4:256588/1‑149 (MQ=255)
                                                                  tGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATc                                                                                    <  3:368576/149‑1 (MQ=255)
                                                                   ggATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCg                                                                                   >  3:558738/1‑149 (MQ=255)
                                                                       gggCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTAc                                                                               <  3:133782/149‑1 (MQ=255)
                                                                            cgTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGAt                                                                          <  2:346790/149‑1 (MQ=255)
                                                                                ttCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGttt                                                                      <  4:134194/149‑1 (MQ=255)
                                                                                  ctctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgt                                                                    <  1:292106/149‑1 (MQ=255)
                                                                                  ctctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgt                                                                    <  4:558738/149‑1 (MQ=255)
                                                                                  ctctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgt                                                                    <  2:183635/149‑1 (MQ=255)
                                                                                  cactGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgt                                                                    <  3:617820/147‑1 (MQ=255)
                                                                                   tctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCTATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgtg                                                                   >  1:258370/1‑149 (MQ=255)
                                                                                   tctGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTtgtg                                                                   <  4:323688/149‑1 (MQ=255)
                                                                                         gATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATc                                                                                         <  4:87523/121‑1 (MQ=255)
                                                                                         gATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATc                                                                                         >  3:87523/1‑121 (MQ=255)
                                                                                             ctctTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTaa                                                         <  2:376291/149‑1 (MQ=255)
                                                                                             ctctTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTaa                                                         <  4:265276/149‑1 (MQ=255)
                                                                                             ctctTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTaa                                                         <  3:387419/149‑1 (MQ=255)
                                                                                              tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTTTGTTATCTTTAat                                                        >  3:286390/1‑149 (MQ=255)
                                                                                              tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAat                                                        >  1:159695/1‑149 (MQ=255)
                                                                                              tcttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCCTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAat                                                        >  3:162666/1‑149 (MQ=255)
                                                                                                ttctGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAatat                                                      >  3:31727/1‑149 (MQ=255)
                                                                                                   tGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCa                                                   <  3:483920/149‑1 (MQ=255)
                                                                                                      gCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCAccc                                                >  3:721937/1‑149 (MQ=255)
                                                                                                      gCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGCTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCAccc                                                <  1:701630/149‑1 (MQ=255)
                                                                                                        cGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTg                                              >  3:54324/1‑149 (MQ=255)
                                                                                                          aaCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTaa                                                         >  3:415124/1‑136 (MQ=255)
                                                                                                          aaCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTaa                                                         <  4:415124/136‑1 (MQ=255)
                                                                                                          aaCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGaa                                            >  3:373591/1‑149 (MQ=255)
                                                                                                           aCTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAag                                           <  4:668053/149‑1 (MQ=255)
                                                                                                             tcctcTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAagag                                         >  1:150390/1‑149 (MQ=255)
                                                                                                                tctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAAt                                                                                     <  3:41135/102‑1 (MQ=255)
                                                                                                                tctGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAAt                                                                                     >  4:41135/1‑102 (MQ=255)
                                                                                                                  tGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCa                                    <  2:150390/149‑1 (MQ=255)
                                                                                                                   ggATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAg                                   <  1:565055/149‑1 (MQ=255)
                                                                                                                   ggATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAg                                   <  3:219685/149‑1 (MQ=255)
                                                                                                                       cATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTc                                                    <  3:452478/128‑1 (MQ=255)
                                                                                                                       cATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTc                                                    >  4:452478/1‑128 (MQ=255)
                                                                                                                        aTTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCtt                              <  3:53570/149‑1 (MQ=255)
                                                                                                                         ttACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTc                             >  2:254615/1‑149 (MQ=255)
                                                                                                                          tACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCg                            >  1:773194/1‑149 (MQ=255)
                                                                                                                               aCTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAAcc                       >  3:529222/1‑149 (MQ=255)
                                                                                                                               aCTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAAcc                       >  1:155629/1‑149 (MQ=255)
                                                                                                                                 tGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCt                     <  1:507091/149‑1 (MQ=255)
                                                                                                                                     tgacaaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCa                 >  1:113521/6‑149 (MQ=255)
                                                                                                                                      cTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCAt                >  2:668167/1‑149 (MQ=255)
                                                                                                                                      cTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCAt                >  4:512194/1‑149 (MQ=255)
                                                                                                                                      cTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCAt                >  3:484782/1‑149 (MQ=255)
                                                                                                                                         aaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCttt                                                           <  4:344005/102‑1 (MQ=255)
                                                                                                                                         aaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCttt                                                           >  3:344005/2‑103 (MQ=255)
                                                                                                                                         aaCCCCGCCGCAGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTaa             >  4:177208/2‑149 (MQ=255)
                                                                                                                                            cccGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAagga          >  3:527967/1‑149 (MQ=255)
                                                                                                                                             ccGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAaggag         >  1:375578/1‑149 (MQ=255)
                                                                                                                                                    ggcggGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTAt  >  2:8006/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCAACACAAAATAATAACCGTTAAGGGTGTAGCCTATGATCAACACAAATATGAAATATTGGTCCTGGATGGGCGCGTTTTCTCTGTCGATGCTCTTCTGGGCCGAACTCCTCTGGATCATTACTCACTGATCCTTGACCCCGCTGCGGCGGGGTTGTCATTTGCTTTGCCACAAGGTTTCTCCTCTTTTATCAATTTGGGTTGTTATCAAATCGTTACGCGATGTTTGTGTTATCTTTAATATTCACCCTGAAGAGAATCAGGGCTTCGCAACCCTGTCATTAAGGAGGAGCTAT  >  NC_000913/1335536‑1335832

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