Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,416,555 A→T 13.5% intergenic (+88/+27) pdxB → / ← PP_2118 erythronate‑4‑phosphate dehydrogenase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,416,5550AT13.5% 223.6 / 21.4 88intergenic (+88/+27)pdxB/PP_2118erythronate‑4‑phosphate dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (40/36);  new base T (6/6);  total (46/42)
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.90e-01

GGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTGCG  >  NC_002947/2416411‑2416703
                                                                                                                                                 |                                                                                                                                                    
ggggAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTgcgc                                                                                                                                                  >  1:328081/1‑149 (MQ=255)
   gAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCgg                                                                                                                                               <  5:410776/149‑1 (MQ=255)
   gAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCgg                                                                                                                                               <  7:146260/149‑1 (MQ=255)
    aGGCTCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCggg                                                                                                                                              <  5:413677/149‑1 (MQ=255)
    aGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCggg                                                                                                                                              <  2:328081/149‑1 (MQ=255)
      gcgcCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGtt                                                                                                                                            <  6:154172/149‑1 (MQ=255)
       cgccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGttt                                                                                                                                           <  7:219857/149‑2 (MQ=255)
       cgccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGttt                                                                                                                                           <  7:342535/149‑2 (MQ=255)
        gccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGc                                                                                                                                          <  4:19788/149‑1 (MQ=255)
         ccgcAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCa                                                                                                                                         <  4:108650/149‑1 (MQ=255)
           gcagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGtttttat                                                                                                                                       >  4:341949/1‑144 (MQ=255)
            cagttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTTGCGCCGGGtttttatt                                                                                                                                      <  3:341949/149‑7 (MQ=255)
             agttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCagag                                                                                                                                      <  3:253947/148‑1 (MQ=255)
             agttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCagag                                                                                                                                      >  4:253947/1‑148 (MQ=255)
              gttggcgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGaa                                                                                                                                    <  2:208688/149‑1 (MQ=255)
                   cgcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGc                                                                                                                               >  2:193885/1‑149 (MQ=255)
                    gcagttggTCAGTGCCGTGGGTGGGGTGTTGGTCTTAGCACGCGTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGcc                                                                                                                              <  7:24485/149‑1 (MQ=255)
                    gcagttggTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGcc                                                                                                                              >  1:180867/1‑149 (MQ=255)
                    gcagttggTCAGTGCACTGGGTGGGGTGTGGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGcc                                                                                                                              <  5:66628/149‑1 (MQ=255)
                            tCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAAcc                                                                                                                      >  5:269989/1‑149 (MQ=255)
                            tCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGACCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAAcc                                                                                                                      >  7:15442/1‑149 (MQ=255)
                                       ggtggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTAcgc                                                                                                           <  4:388496/149‑1 (MQ=255)
                                       ggtggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTAcgc                                                                                                           >  8:145258/1‑149 (MQ=255)
                                       ggtggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTAcg                                                                                                            <  3:194188/148‑1 (MQ=255)
                                       ggtggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTAc                                                                                                             <  2:180867/147‑1 (MQ=255)
                                         tggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACgcgc                                                                                                         >  6:185828/1‑149 (MQ=255)
                                         tggggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACgcgc                                                                                                         <  2:420593/149‑1 (MQ=255)
                                           gggtgTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTc                                                                                                       <  2:303331/149‑1 (MQ=255)
                                             gtgtTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGa                                                                                                     >  4:302674/1‑149 (MQ=255)
                                              tgtTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGac                                                                                                    <  6:269989/149‑1 (MQ=255)
                                              tgtTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGa                                                                                                     >  7:46521/1‑148 (MQ=255)
                                                    tCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCt                                                                                              >  7:386665/1‑149 (MQ=255)
                                                    tCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCt                                                                                              <  2:405085/149‑1 (MQ=255)
                                                       aaGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCttt                                                                                            <  7:145258/148‑1 (MQ=255)
                                                        aGCACGCGTTGTCTTGACCTTGGAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGTGCGAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCttt                                                                                            <  1:378624/147‑1 (MQ=255)
                                                            cGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTCGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCa                                                                                      <  1:193885/149‑1 (MQ=255)
                                                              ccTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTAGCCGATCAACATTTGCGTACCCGCTCGACACGGCTTTCTTCCAgg                                                                                    >  8:117980/1‑149 (MQ=255)
                                                                   tCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGcc                                                                                 >  1:156837/1‑147 (MQ=255)
                                                                    tctgtccTTGCACGTAAGGTGGGGTAACTGGTGGGGGGGGCAACGCCAATGGGCGGGCCAAAAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCtt                                                                               <  6:139503/143‑1 (MQ=255)
                                                                       gACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAg                                                                           >  3:125316/1‑149 (MQ=255)
                                                                        aCCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGACCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAg                                                                           <  8:15442/148‑1 (MQ=255)
                                                                                      gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTgg                                                                  <  2:298224/144‑1 (MQ=255)
                                                                                      gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTgg                                                                  >  1:298224/1‑144 (MQ=255)
                                                                                      gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcat                                                             <  5:169036/149‑1 (MQ=255)
                                                                                      gtgtGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatca                                                              >  4:244233/1‑148 (MQ=255)
                                                                                       tgtgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGatcatc                                                            <  3:302674/149‑1 (MQ=255)
                                                                                        gtgGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCt                                                           >  6:135365/1‑149 (MQ=255)
                                                                                          ggTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGc                                                         >  5:161507/1‑149 (MQ=255)
                                                                                           gTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCt                                                        >  1:179639/1‑149 (MQ=255)
                                                                                            ttACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTc                                                       >  1:8704/1‑149 (MQ=255)
                                                                                              aCCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCgg                                                     >  1:159651/1‑149 (MQ=255)
                                                                                                cTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGt                                                    >  3:120413/1‑148 (MQ=255)
                                                                                                cTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGt                                                    >  6:4275/1‑148 (MQ=255)
                                                                                                  gTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGAt                                                 >  1:377808/1‑149 (MQ=255)
                                                                                                   ttGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATc                                                <  8:386665/149‑1 (MQ=255)
                                                                                                   ttGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATc                                                >  2:166994/1‑149 (MQ=255)
                                                                                                     ggggCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCgg                                              >  7:392546/1‑149 (MQ=255)
                                                                                                      gggCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCgg                                              <  3:244233/148‑1 (MQ=255)
                                                                                                        gcggcCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCTAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCAc                                           >  4:327307/1‑149 (MQ=255)
                                                                                                            ccAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATc                                                >  2:376992/1‑140 (MQ=255)
                                                                                                            ccAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATc                                                <  1:376992/140‑1 (MQ=255)
                                                                                                              aaGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTcgcg                                     <  5:4275/149‑1 (MQ=255)
                                                                                                                gCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTTCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcc                                   >  4:73244/1‑149 (MQ=255)
                                                                                                                gCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGTATCATCTGCTCGGTGATCGGCACTTCGCGcc                                   >  5:116138/1‑149 (MQ=255)
                                                                                                                 cccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGccc                                  >  4:111816/1‑149 (MQ=255)
                                                                                                                  ccACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGcccc                                 >  3:7644/1‑149 (MQ=255)
                                                                                                                   cACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCt                                >  6:88094/1‑149 (MQ=255)
                                                                                                                    aCTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCtt                               <  1:166994/149‑1 (MQ=255)
                                                                                                                            ggCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTc                         <  2:159651/147‑1 (MQ=255)
                                                                                                                            ggCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCg                        <  2:156837/148‑1 (MQ=255)
                                                                                                                              ccAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgatg                     <  2:179639/149‑1 (MQ=255)
                                                                                                                               cAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGTATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgatga                    <  6:116138/149‑1 (MQ=255)
                                                                                                                               cAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgttgg                    >  1:156399/1‑148 (MQ=255)
                                                                                                                                aaTAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCgatgat                   >  3:85016/1‑149 (MQ=255)
                                                                                                                                 tctgcaaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                  >  4:333741/6‑149 (MQ=255)
                                                                                                                                 tctgcaaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                  >  8:253043/6‑149 (MQ=255)
                                                                                                                                  tAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATc                  >  5:326861/1‑148 (MQ=255)
                                                                                                                                     aaaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCt                                                        >  3:177421/1‑107 (MQ=255)
                                                                                                                                     aaaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCt                                                        <  4:177421/107‑1 (MQ=255)
                                                                                                                                     aaaCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAAc              <  1:54366/149‑1 (MQ=255)
                                                                                                                                      aaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAAcc             >  8:355089/1‑149 (MQ=255)
                                                                                                                                      aaCCCGGCGCTTGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTGGGCGTCGATGATCGAAcc             >  8:18627/1‑149 (MQ=255)
                                                                                                                                       aCCCGGCGCAAGCGCCGGGTTGCAGAGAATTCCCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCg            <  2:412938/149‑1 (MQ=255)
                                                                                                                                          cGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCa          >  3:115491/1‑148 (MQ=255)
                                                                                                                                              gcAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGATACGGCTTTCTTCCAGGGCCTTGCAGGCCTTCTGGATAATCTGATCTGTGATCGTCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACtt     >  4:236706/1‑149 (MQ=255)
                                                                                                                                              gcAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTc                                        >  7:129311/1‑114 (MQ=255)
                                                                                                                                              gcAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTc                                        <  8:129311/114‑1 (MQ=255)
                                                                                                                                              gcAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACtt     >  5:70429/1‑149 (MQ=255)
                                                                                                                                                 aGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTGCg  <  3:111816/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGGGAGGCGCCGCAGTTGGCGCAGTTGGTCAGTGCCCTGGGTGGGGTGTTGGTCTAAGCACGCCTTGTCTTGACCTTGCAGGTAGGTGTGGTTACCTGTTGGGGCGGCCAAGCCCACTGGCTGGGCCAATAAAAACCCGGCGCAAGCGCCGGGTTGCAGAGAATTCGCCGATCAACCTTTGCGTACGCGCTCGACACGGCTTTCTTCCAGGGCCTTGCAGGCCTGCTGGATCATCTGCTCGGTGATCGGCACTTCGCGCCCCTTGGCGTCGATGATCGAACCGCATACTTGCG  >  NC_002947/2416411‑2416703

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