Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,737,707 C→T 5.4% intergenic (+61/‑38) pheL → / → pheA pheA gene leader peptide/chorismate mutase and prephenate dehydratase, P‑protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,737,7070CT 5.4% 306.3 / 4.3 92intergenic (+61/‑38)pheL/pheApheA gene leader peptide/chorismate mutase and prephenate dehydratase, P‑protein
Reads supporting (aligned to +/- strand):  ref base C (49/38);  new base T (3/2);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAA  >  NC_000913/2737562‑2737853
                                                                                                                                                 |                                                                                                                                                  
tACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCggg                                                                                                                                                 <  3:537273/149‑1 (MQ=255)
tACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCggg                                                                                                                                                 <  4:108931/149‑1 (MQ=255)
tACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCggg                                                                                                                                                 <  4:20330/149‑1 (MQ=255)
tACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgg                                                                                                                                                  <  2:376114/148‑1 (MQ=255)
tACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgg                                                                                                                                                  <  2:738348/148‑1 (MQ=255)
     gtactgtactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCCCCCGATTTGGGAGGcc                                                                                                                                            <  4:290114/149‑1 (MQ=255)
           tactAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCtttttt                                                                                                                                      >  3:869071/1‑149 (MQ=255)
             ctAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAt                                                                                                                                    >  4:417167/1‑149 (MQ=255)
              tAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgggggg                                                                                                                                              >  2:43182/1‑138 (MQ=255)
              tAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCgggggg                                                                                                                                              <  1:43182/138‑1 (MQ=255)
              tAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAtt                                                                                                                                   <  3:227545/149‑1 (MQ=255)
              tAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAtt                                                                                                                                   <  4:869071/149‑1 (MQ=255)
               aaaGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAtt                                                                                                                                   >  1:407064/1‑148 (MQ=255)
                   tCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTTCCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATaa                                                                                                                              <  3:215718/149‑1 (MQ=255)
                   tCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATaa                                                                                                                              <  3:843160/149‑1 (MQ=255)
                        tAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaa                                                                                                                          >  1:618794/1‑148 (MQ=255)
                        tAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCCCCCGATTTGGGAGGc                                                                                                                                             >  2:256010/1‑129 (MQ=255)
                         aaGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaaaa                                                                                                                        <  4:283276/149‑1 (MQ=255)
                         aaGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaaaa                                                                                                                        >  3:400053/1‑149 (MQ=255)
                          aGGAAACAAACATGAAACACATACCGTTTTTCTTCGCGTTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAg                                                                                                                       >  3:536164/1‑149 (MQ=255)
                          aGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAg                                                                                                                       >  4:772948/1‑149 (MQ=255)
                           ggAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAg                                                                                                                       >  2:363020/1‑148 (MQ=255)
                           ggAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAAAAAAAAg                                                                                                                       <  2:494275/148‑1 (MQ=255)
                            gAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAt                                                                                                                                    <  1:391050/134‑1 (MQ=255)
                            gAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAt                                                                                                                                    >  2:391050/1‑134 (MQ=255)
                             aaacaaacATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGc                                                                                                                     >  1:410770/1‑148 (MQ=255)
                             aaacaaacATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCa                                                                                                                    >  3:518905/1‑149 (MQ=255)
                              aacaaacaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCa                                                                                                                    <  2:407064/148‑1 (MQ=255)
                                  aacaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaaaa                                                                                                                        <  2:570330/140‑1 (MQ=255)
                                  aacaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACaaaaa                                                                                                                        >  1:570330/1‑140 (MQ=255)
                                   acaTGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTa                                                                                                              >  3:237397/1‑149 (MQ=255)
                                        aaaCACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAt                                                                                                                                    >  2:35746/1‑122 (MQ=255)
                                        aaaCACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTAt                                                                                                                                    <  1:35746/122‑1 (MQ=255)
                                              aTACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGa                                                                                                    <  1:202266/148‑1 (MQ=255)
                                                  cGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAgg                                                                                                                      <  2:150558/126‑1 (MQ=255)
                                                  cGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAgg                                                                                                                      >  1:150558/1‑126 (MQ=255)
                                                          ttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACAt                                                                                                        <  3:9612/132‑1 (MQ=255)
                                                          ttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACAt                                                                                                        >  4:9612/1‑132 (MQ=255)
                                                          ttcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTg                                                                                       <  3:417167/149‑1 (MQ=255)
                                                           tcGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTg                                                                                       >  1:632514/1‑148 (MQ=255)
                                                              cATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCCCCCGATTTGGGAGGCCt                                                                                                                                           >  3:583824/1‑93 (MQ=38)
                                                              cATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCCCCCGATTTGGGAGGCCt                                                                                                                                           <  4:583824/93‑1 (MQ=38)
                                                                ttCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCt                                                                                  >  1:432569/1‑148 (MQ=255)
                                                                 tCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTg                                                                                 >  2:706654/1‑148 (MQ=255)
                                                                 tCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGc                                                                                <  4:400053/149‑1 (MQ=255)
                                                                   ttttttACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCga                                                                              <  4:536164/149‑1 (MQ=255)
                                                                          ccTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGaaa                                                                        <  1:706654/148‑1 (MQ=255)
                                                                           cTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGAGATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt                                                                      >  4:611283/1‑149 (MQ=255)
                                                                            ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt                                                                      >  2:108958/1‑148 (MQ=255)
                                                                            ttCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAAt                                                                      >  1:593225/1‑148 (MQ=255)
                                                                               cccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACTCTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAgc                                                                  <  3:491646/149‑1 (MQ=255)
                                                                               cccTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAg                                                                   <  2:410770/148‑1 (MQ=255)
                                                                                      tGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTg                                                            >  1:493028/1‑148 (MQ=255)
                                                                                          aGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGa                                                       >  4:771353/1‑149 (MQ=255)
                                                                                                  cgtcgtGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAatta                                                <  1:108958/148‑1 (MQ=255)
                                                                                                       tgtgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGc                                           >  2:38195/1‑148 (MQ=255)
                                                                                                       tgtgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCg                                          >  4:770976/1‑149 (MQ=255)
                                                                                                         tgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGtt                                        >  4:266162/1‑149 (MQ=255)
                                                                                                         tgAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGtt                                        >  3:319070/1‑149 (MQ=255)
                                                                                                           aaaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTAc                                      >  3:217175/1‑149 (MQ=255)
                                                                                                           aaaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTAc                                      >  3:112290/1‑149 (MQ=255)
                                                                                                            aaCAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTAc                                      >  1:87549/1‑148 (MQ=255)
                                                                                                              cAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTg                                    >  1:460673/1‑148 (MQ=255)
                                                                                                                       aaGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACgg                           >  2:513019/1‑148 (MQ=255)
                                                                                                                       aaGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACgg                           <  2:460673/148‑1 (MQ=255)
                                                                                                                        aGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGCTAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCt                                                             >  2:550874/1‑113 (MQ=255)
                                                                                                                        aGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCt                                                             <  1:550874/113‑1 (MQ=255)
                                                                                                                        aGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGgcg                         >  4:85628/1‑149 (MQ=255)
                                                                                                                          aCGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACGATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGcgcg                       <  4:93959/149‑1 (MQ=255)
                                                                                                                          aCGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACGATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGcgcg                       >  3:93959/1‑149 (MQ=255)
                                                                                                                            gAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGaa                     <  3:772948/149‑1 (MQ=255)
                                                                                                                            gAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGaa                     >  4:660095/1‑149 (MQ=255)
                                                                                                                             ataaaaaaGGCCCCCCGATTTGGGAGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAAc                    >  3:279694/7‑149 (MQ=255)
                                                                                                                                 aTAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTgg                 >  1:270385/1‑148 (MQ=255)
                                                                                                                                       ccTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCg           <  2:618794/148‑1 (MQ=255)
                                                                                                                                       ccTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCg           >  2:689997/1‑148 (MQ=255)
                                                                                                                                        cTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAg         >  4:526557/1‑149 (MQ=255)
                                                                                                                                        cTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGACGTCGa          <  2:493028/148‑1 (MQ=255)
                                                                                                                                         tCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAatta                                                >  4:137882/1‑109 (MQ=255)
                                                                                                                                         tCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAatta                                                <  3:137882/109‑1 (MQ=255)
                                                                                                                                          cccAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAgg        >  2:900226/1‑148 (MQ=255)
                                                                                                                                          cccAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAgg        >  1:696952/1‑148 (MQ=255)
                                                                                                                                           ccAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTg      >  3:518054/1‑149 (MQ=255)
                                                                                                                                           ccAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTg      <  4:112290/149‑1 (MQ=255)
                                                                                                                                            cAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTgg     <  3:660095/149‑1 (MQ=255)
                                                                                                                                            cAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTgg     <  3:771353/149‑1 (MQ=255)
                                                                                                                                            cAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTg      >  1:411499/1‑148 (MQ=255)
                                                                                                                                            cAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTg      >  2:830077/1‑148 (MQ=255)
                                                                                                                                             aaaTCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTggg    <  3:85628/149‑1 (MQ=255)
                                                                                                                                                tCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGa          >  2:250634/1‑140 (MQ=255)
                                                                                                                                                tCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGa          <  1:250634/140‑1 (MQ=255)
                                                                                                                                                tCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGaa  >  2:637381/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TACGGGTACTGTACTAAAGTCACTTAAGGAAACAAACATGAAACACATACCGTTTTTCTTCGCATTCTTTTTTACCTTCCCCTGAATGGGAGGCGTTTCGTCGTGTGAAACAGAATGCGAAGACGAACAATAAGGCCTCCCAAATCGGGGGGCCTTTTTTATTGATAACAAAAAGGCAACACTATGACATCGGAAAACCCGTTACTGGCGCTGCGAGAGAAAATCAGCGCGCTGGATGAAAAATTATTAGCGTTACTGGCAGAACGGCGCGAACTGGCCGTCGAGGTGGGAA  >  NC_000913/2737562‑2737853

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