Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,087,892 C→G 6.9% intergenic (+32/‑392) metK → / → galP S‑adenosylmethionine synthetase/D‑galactose transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,087,8920CG 6.9% 331.9 / 11.7 100intergenic (+32/‑392)metK/galPS‑adenosylmethionine synthetase/D‑galactose transporter
Reads supporting (aligned to +/- strand):  ref base C (44/48);  new base G (2/5);  total (46/54)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 101 total aligned reads displayed.

GCTGCACCCGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAG  >  NC_000913/3087746‑3088037
                                                                                                                                                  |                                                                                                                                                 
gctgcACCCGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTcgc                                                                                                                                                 >  2:131453/1‑149 (MQ=255)
    cACCCGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTg                                                                                                                                             <  5:14349/149‑1 (MQ=255)
     aCCCGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTgg                                                                                                                                            <  7:167542/149‑1 (MQ=255)
      cccGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGc                                                                                                                                           <  3:356504/149‑1 (MQ=255)
       ccGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGcc                                                                                                                                          >  3:194678/1‑149 (MQ=255)
       ccGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGcc                                                                                                                                          >  2:118576/1‑149 (MQ=255)
          ttgtggttggtgtaatcggtttcagttgtgatgataatctgcggtgaaacgtgatcaggaagaatgatgttattgtgcgggcggtttttttgccggatgcggcgtgaacgcctatccaaaagaaAAAGGCCAGCCTCGCGCTGGCCttt                                                                                                                                       >  8:93108‑M1/125‑149 (MQ=255)
              aCAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCttt                                                                                                                                       >  5:456536/1‑145 (MQ=255)
               cAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCGCGAGGCTGGCCtttgaa                                                                                                                                    <  4:259420/147‑4 (MQ=255)
               cAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCtttt                                                                                                                                      <  1:301438/145‑1 (MQ=255)
               cAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCttt                                                                                                                                       <  1:379771/144‑1 (MQ=255)
                aaaGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCttt                                                                                                                                 <  1:245216/149‑1 (MQ=255)
                 aaGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCtttt                                                                                                                                >  4:372315/1‑149 (MQ=255)
                 aaGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTACCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGGCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCtttt                                                                                                                                <  5:193801/149‑1 (MQ=255)
                  aGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTg                                                                                                                               <  2:259170/149‑1 (MQ=255)
                  aGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTg                                                                                                                               >  4:137810/1‑149 (MQ=255)
                   gAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTgg                                                                                                                              <  1:97283/149‑1 (MQ=255)
                      cccGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATa                                                                                                                           <  2:96349/148‑1 (MQ=255)
                      aCCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGa                                                                                                                             >  1:352134/1‑147 (MQ=255)
                      aCCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGa                                                                                                                             >  1:246766/1‑147 (MQ=255)
                      aCCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATa                                                                                                                           <  2:246766/149‑1 (MQ=255)
                            gcaTACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCggt                                                                                                                     <  1:163488/149‑3 (MQ=255)
                                 cGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCAcgc                                                                                                                >  4:185496/1‑149 (MQ=255)
                                 cGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCAcgc                                                                                                                <  8:4483/149‑1 (MQ=255)
                                 cGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCAcgc                                                                                                                >  6:13971/1‑149 (MQ=255)
                                  ggTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGTCGTTCAcgcc                                                                                                               >  3:140626/1‑149 (MQ=255)
                                  ggTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCAcgcc                                                                                                               <  3:137810/149‑1 (MQ=255)
                                   gTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCGCGAGGCTGGCCtttgaac                                                                                                                                   <  2:215164/128‑5 (MQ=255)
                                   gTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCGCGAGGCTGGCCtttgaac                                                                                                                                   >  1:215164/1‑124 (MQ=255)
                                       cTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTg                                                                                                                               <  7:84101/128‑1 (MQ=255)
                                       cTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTg                                                                                                                               >  8:84101/1‑128 (MQ=255)
                                            gTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGc                                                                                                      <  1:118576/148‑1 (MQ=255)
                                                tGAACATTTCCCGTGGTAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGTCaaaa                                                                                                  >  4:304183/1‑148 (MQ=255)
                                                  aaCATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCtt                                                                                                                                        >  3:3300/1‑108 (MQ=255)
                                                  aaCATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCtt                                                                                                                                        <  4:3300/108‑1 (MQ=255)
                                                   aCATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAAc                                                                                              <  4:130117/149‑1 (MQ=255)
                                                     aTTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAAccg                                                                                            <  8:293396/149‑1 (MQ=255)
                                                          ccGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAAccgcccg                                                                                        >  4:372915/1‑148 (MQ=255)
                                                               ggAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACaa                                                                                   <  7:253707/148‑1 (MQ=255)
                                                               ggAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACaa                                                                                   >  5:447383/1‑148 (MQ=255)
                                                                 aaaaaaCCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATaa                                                                                >  1:84720/1‑149 (MQ=255)
                                                                     aaCCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAAcatc                                                                            <  4:194678/149‑1 (MQ=255)
                                                                         gACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGca                                                                                      <  2:285878/135‑1 (MQ=255)
                                                                         gACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGca                                                                                      >  1:285878/1‑135 (MQ=255)
                                                                           cAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAttctt                                                                      <  5:13971/149‑1 (MQ=255)
                                                                            aaaGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAttcttc                                                                     <  3:372315/149‑1 (MQ=255)
                                                                             aaGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTcc                                                                    <  6:447383/149‑1 (MQ=255)
                                                                                cgcAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGa                                                                 >  4:315675/1‑149 (MQ=255)
                                                                                   aGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCa                                                              <  1:131453/149‑1 (MQ=255)
                                                                                       gctgcGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGtt                                                          >  8:452672/1‑149 (MQ=255)
                                                                                         tgcGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGttt                                                         >  7:33657/1‑148 (MQ=255)
                                                                                          gcgcGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCa                                                       >  7:79067/1‑149 (MQ=255)
                                                                                            gcgATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCAcc                                                     >  6:377703/1‑149 (MQ=255)
                                                                                                  ctgcCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAAcatca                                                                           <  8:239502/121‑1 (MQ=255)
                                                                                                  ctgcCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAAcatca                                                                           >  7:239502/1‑121 (MQ=255)
                                                                                                   tgcCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCCCCGCAGAt                                               >  4:419992/1‑148 (MQ=255)
                                                                                                      cGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTat                                            >  6:219361/1‑148 (MQ=255)
                                                                                                        gTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatcat                                         >  7:99431/1‑149 (MQ=255)
                                                                                                           tGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAAca                                                                              >  6:332791/1‑109 (MQ=255)
                                                                                                           tGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAAca                                                                              <  5:332791/109‑1 (MQ=255)
                                                                                                           tGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTAtcatca                                       >  4:368969/1‑148 (MQ=255)
                                                                                                           tGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTAtcatca                                       >  2:152138/1‑148 (MQ=255)
                                                                                                            gAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAAcc                                                                                             >  7:462439/1‑93 (MQ=255)
                                                                                                            gAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAAcc                                                                                             <  8:462439/93‑1 (MQ=255)
                                                                                                             aaGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGc                                                   <  3:257806/134‑1 (MQ=255)
                                                                                                             aaGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGc                                                   >  4:257806/1‑134 (MQ=255)
                                                                                                              ggcgtgaacgcctatccaaaagaaAAAGGCCAGCGCGAGGCTGGCCTTTTTCTTTTGGATAGGCCTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAAAAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACaa                                    <  6:226986‑M1/124‑1 (MQ=255)
                                                                                                              aGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCAc                                                      <  5:375097/130‑1 (MQ=255)
                                                                                                              aGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCAc                                                      >  6:375097/1‑130 (MQ=255)
                                                                                                               gTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACt                                  >  4:53959/1‑149 (MQ=255)
                                                                                                               gTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACt                                  >  3:196318/1‑149 (MQ=255)
                                                                                                                tAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc                                           <  8:174630/139‑1 (MQ=255)
                                                                                                                tAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTatc                                           >  7:174630/1‑139 (MQ=255)
                                                                                                                  aTCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg                                                 >  8:120955/1‑131 (MQ=255)
                                                                                                                  aTCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAg                                                 <  7:120955/131‑1 (MQ=255)
                                                                                                                   tctttcttCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGaaa                              >  8:353069/1‑149 (MQ=255)
                                                                                                                        cttcACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGa                          <  8:33657/148‑1 (MQ=255)
                                                                                                                         ttcACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGa                          <  2:44716/147‑1 (MQ=255)
                                                                                                                         ttcACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGa                          >  1:44716/1‑147 (MQ=255)
                                                                                                                         ctatccaaaagaaAAAGGCCAGCGCGAGGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGAt                         <  2:381327‑M1/135‑1 (MQ=255)
                                                                                                                           cACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGAtt                        <  1:152138/147‑1 (MQ=255)
                                                                                                                           cACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTa                       <  4:208867/148‑1 (MQ=255)
                                                                                                                             ccTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTacac                    >  8:153305/1‑149 (MQ=255)
                                                                                                                              cTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAttct                                                                       <  1:316776/97‑1 (MQ=255)
                                                                                                                              cTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAttct                                                                       >  2:316776/1‑97 (MQ=255)
                                                                                                                              cTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACacc                   <  4:196318/149‑1 (MQ=255)
                                                                                                                               aaaagaaAAAGGCCAGCGCGAGGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAttcttc                                                                     <  7:139566‑M1/91‑1 (MQ=255)
                                                                                                                               aaaagaaAAAGGCCAGCGCGAGGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCAttcttc                                                                     >  8:139566‑M1/8‑98 (MQ=255)
                                                                                                                               aaaagaaAAAGGCCAGCGCAAGGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACacc                   <  1:195456‑M1/141‑1 (MQ=255)
                                                                                                                                 cGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACaccaac                <  3:372915/149‑1 (MQ=255)
                                                                                                                                  gTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATAACGTTTCACCGCAGATTATCATCACAACGGAAACCGAATACaccaac                >  2:342773/1‑148 (MQ=255)
                                                                                                                                   ttCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACaccaacca              <  3:87413/149‑1 (MQ=255)
                                                                                                                                      AAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAAc                             >  2:180754‑M1/1‑131 (MQ=255)
                                                                                                                                      AAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAAc                             <  1:180754‑M1/131‑1 (MQ=255)
                                                                                                                                      AAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCacaa           <  6:65404‑M1/149‑1 (MQ=255)
                                                                                                                                       AAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCacaa           <  1:298690‑M1/148‑1 (MQ=255)
                                                                                                                                       AAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCacaa           <  2:385438‑M1/148‑1 (MQ=255)
                                                                                                                                       AAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCacaa           <  2:276453‑M1/148‑1 (MQ=255)
                                                                                                                                            CAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGAc      <  8:454454‑M1/148‑1 (MQ=255)
                                                                                                                                               CCTCGCGCTGGCCTTTTTCTTTTTGATAGGCGTTCACGCAGCATGCGGCAAAAAATCCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAg  >  4:140626‑M1/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTGCACCCGATCTACAAAGAAACCGCAGCATACGGTCACTTTGGTCGTGAACATTTCCCGTGGGAAAAAACCGACAAAGCGCAGCTGCTGCGCGATGCTGCCGGTCTGAAGTAATCTTTCTTCACCTGCGTTCAAAGGCCAGCCTCGCGCTGGCCTTTTTCTTTTGGATAGGCGTTCACGCCGCATCCGGCAAAAAAACCGCCCGCACAATAACATCATTCTTCCTGATCACGTTTCACCGCAGATTATCATCACAACTGAAACCGATTACACCAACCACAACAGACAAAG  >  NC_000913/3087746‑3088037

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