Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,352,025 G→T 10.8% intergenic (‑230/+34) PP_4707 ← / ← pnp OsmY‑like protein/polyribonucleotide nucleotidyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,352,0250GT10.8% 268.9 / 16.8 93intergenic (‑230/+34)PP_4707/pnpOsmY‑like protein/polyribonucleotide nucleotidyltransferase
Reads supporting (aligned to +/- strand):  ref base G (47/35);  new base T (5/5);  total (52/41)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

TTTGCTGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCG  >  NC_002947/5351878‑5352173
                                                                                                                                                   |                                                                                                                                                    
tttGCTGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCg                                                                                                                                                    >  2:2329310/1‑150 (MQ=255)
tttGCTGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCg                                                                                                                                                    >  1:1512275/1‑150 (MQ=255)
   gctgctTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAg                                                                                                                                                 <  1:514937/150‑1 (MQ=255)
    ctgntTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAgg                                                                                                                                                >  2:1501696/1‑150 (MQ=255)
    ctgctTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGNAgg                                                                                                                                                <  2:1512275/150‑1 (MQ=255)
         tCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTTCGCAGGGCCCt                                                                                                                                           <  1:1158810/150‑1 (MQ=255)
                 ttCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTTCGCAGGGCCCttt                                                                                                                                         <  2:2120916/144‑1 (MQ=255)
                 ttCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTTCGCAGGGCCCttt                                                                                                                                         >  1:2120916/1‑144 (MQ=255)
                  tctttttAAATCAGTAAGTTAGCGAATATCTATCACTTCGGTGGTCGTGAATTTTGCAATCAGCGCGGTAAACCTGCATGAAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTAtt                                                                                                                                  <  2:945318/147‑1 (MQ=255)
                        tAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGa                                                                                                                            <  1:801877/150‑1 (MQ=255)
                              aaTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTTCGCAGGGCCCtttt                                                                                                                                        <  1:2305570/132‑2 (MQ=255)
                              aaTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTTCGCAGGGCCCtttt                                                                                                                                        >  2:2305570/1‑131 (MQ=255)
                               aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                     <  1:1896342/150‑1 (MQ=255)
                               aTAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGc                                                                                                                     >  2:1420967/1‑150 (MQ=255)
                                tAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTNGCAATCNGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCt                                                                                                                    >  1:2162951/1‑150 (MQ=255)
                                   gTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAg                                                                                                                 >  2:1491127/1‑150 (MQ=255)
                                       gCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAg                                                                                                                                                 >  1:404376/1‑114 (MQ=255)
                                       gCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAg                                                                                                                                                 <  2:404376/114‑1 (MQ=255)
                                       gCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGc                                                                                                             >  1:685205/1‑150 (MQ=255)
                                         aaaTATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGccc                                                                                                           >  2:784480/1‑150 (MQ=255)
                                           atatCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTATCTATAGAGTTTCTTGCGCTTAGACGCCCGa                                                                                                         >  1:722195/1‑150 (MQ=255)
                                            tatctatcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAg                                                                                                        <  1:2329310/150‑1 (MQ=255)
                                              tctatcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGc                                                                                                      <  1:1648673/150‑1 (MQ=255)
                                                tatcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCt                                                                                                    >  2:2136957/1‑150 (MQ=255)
                                                tatcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCt                                                                                                    >  1:289050/1‑150 (MQ=255)
                                                tatcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCt                                                                                                    >  1:1755474/1‑150 (MQ=255)
                                                 atcACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCACAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTTACGTTCGCTTAGACGCCCGAGGCCtt                                                                                                   >  1:661089/1‑150 (MQ=255)
                                                    aCTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGGCCGAGGGCTTGGc                                                                                                >  1:405969/1‑150 (MQ=255)
                                                      ttCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTgg                                                                                                 >  2:1069701/1‑147 (MQ=255)
                                                      ttCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTgg                                                                                                 <  1:1069701/147‑1 (MQ=255)
                                                      ttCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCg                                                                                              >  1:52547/1‑150 (MQ=255)
                                                                 tGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCt                                                                                   <  2:685205/150‑1 (MQ=255)
                                                                    aTTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTg                                                                                >  1:597435/1‑150 (MQ=255)
                                                                     ttttGCAATCAGCGCGGTAAACCTGCATGCAAAATGCAAGTCCGCAAAAGGTTCCAGCTCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGATGCCTTGGCCGCATGGACGTCTTTGa                                                                               <  1:982763/150‑1 (MQ=255)
                                                                       ttGCAATCAGCGCGGTAACCCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTCCGAGGGGCCCCTTCGCCGCTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATc                                                                             <  1:292627/150‑1 (MQ=255)
                                                                         gCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGa                                                                           >  1:1069186/1‑150 (MQ=255)
                                                                           aaTCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACa                                                                         >  1:581761/1‑150 (MQ=255)
                                                                           aaTCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACa                                                                         >  1:2360163/1‑150 (MQ=255)
                                                                               aGCGCGGTAAACCGGCATGCAAGATGTAAGTCCGCAAAAGGTTCCCGTTCCCATGAAAAAAGGGCCCTGCGAAGGGCCTTGTCTCGGTTGCGTGACGTTCGCTTAGACGCCTGAGGGCTTGGAGGCAGGGAAGTCTTTAATCGACAGCtt                                                                     <  2:597435/150‑1 (MQ=255)
                                                                                gcgcGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTtctc                                                                                                                                      <  2:1020203/84‑1 (MQ=255)
                                                                                gcgcGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTtctc                                                                                                                                      >  1:1020203/1‑84 (MQ=255)
                                                                                gcgcGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTg                                                                    >  1:370854/1‑150 (MQ=255)
                                                                                 cgcgGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGa                                                                   <  1:2136957/150‑1 (MQ=255)
                                                                                    ggTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAgg                                                                                          >  1:1794000/1‑124 (MQ=255)
                                                                                    ggTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAgg                                                                                          <  2:1794000/124‑1 (MQ=255)
                                                                                    ggTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATAc                                                                >  1:1504425/1‑150 (MQ=255)
                                                                                    ggTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATAc                                                                <  2:1069186/150‑1 (MQ=255)
                                                                                      tAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACgg                                                              <  2:1435581/150‑1 (MQ=255)
                                                                                      tAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACgg                                                              >  1:2202454/1‑150 (MQ=255)
                                                                                        aaCCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGcc                                                            >  1:1438443/1‑150 (MQ=255)
                                                                                           cTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGa                                                                                                                <  1:925730/95‑1 (MQ=255)
                                                                                           cTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGa                                                                                                                >  2:925730/1‑95 (MQ=255)
                                                                                            tgcatgcaAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCgg                                                        >  2:1232725/1‑150 (MQ=255)
                                                                                                 gcaagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTgg                                                                                                 >  1:871155/1‑104 (MQ=255)
                                                                                                 gcaagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTgg                                                                                                 <  2:871155/104‑1 (MQ=255)
                                                                                                 gcaagatgcaagTCCGAAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTAGTGCGTGACGTTCGCTTAGACGCCCGAGGCCATGGCCGCAGGGACGTCTTTGATCGACAGCTTGCTACGGCCGAGGTTGTc                                                   >  2:1225024/1‑150 (MQ=255)
                                                                                                   aagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGAc                                                                                       <  2:128382/112‑1 (MQ=255)
                                                                                                   aagatgcaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGAc                                                                                       >  1:128382/1‑112 (MQ=255)
                                                                                                       tgaaagTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGTCTTCATCGCATTGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTc                                             <  2:1797486/147‑1 (MQ=255)
                                                                                                            gTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCt                                                                                                                    >  2:191830/1‑74 (MQ=255)
                                                                                                            gTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCt                                                                                                                    <  1:191830/74‑1 (MQ=255)
                                                                                                             tCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATa                                                                 <  2:581627/124‑1 (MQ=255)
                                                                                                             tCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATa                                                                 >  1:581627/1‑124 (MQ=255)
                                                                                                               cGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACaa                                     <  1:784480/150‑1 (MQ=255)
                                                                                                                gCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAg                                    >  2:241657/1‑150 (MQ=255)
                                                                                                                gCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAg                                    >  1:2232563/1‑150 (MQ=255)
                                                                                                                      ggTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGACGCGGTTGTCCACGTCCAGTACAAGCACTTc                              <  2:2162951/150‑1 (MQ=255)
                                                                                                                       gTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTGCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCa                             >  2:1235344/1‑150 (MQ=255)
                                                                                                                        ttCCCGTCCCCATGAACGGAGCGCCCGGCGAAGTGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCaa                            <  1:122637/150‑1 (MQ=255)
                                                                                                                          cccGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCtt                                                                                  >  1:2318592/1‑94 (MQ=255)
                                                                                                                          cccGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCtt                                                                                  <  2:2318592/94‑1 (MQ=255)
                                                                                                                                cccATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTgg                    >  1:331802/1‑150 (MQ=255)
                                                                                                                                cacATGAAAAAATGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGATTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTgg                    >  2:1912931/3‑150 (MQ=255)
                                                                                                                                 ccATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGc                   <  2:331802/150‑1 (MQ=255)
                                                                                                                                 ccATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGc                   >  1:2011865/1‑150 (MQ=255)
                                                                                                                                  cATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGcc                  >  1:2249552/1‑150 (MQ=255)
                                                                                                                                  cATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGcc                  >  2:2243567/1‑150 (MQ=255)
                                                                                                                                    tagagaaaGGGCCCTTCGCAGGGCCCTTTTTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCaa                            <  1:582771/133‑1 (MQ=255)
                                                                                                                                    tagagaaaGGGCCCTTCGCAGGGCCCTTTTTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCaa                            >  2:582771/6‑138 (MQ=255)
                                                                                                                                     gAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTc                              >  2:1850321/1‑135 (MQ=255)
                                                                                                                                     gAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTc                              <  1:1850321/135‑1 (MQ=255)
                                                                                                                                        aaaaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGccttcctt            <  2:1755474/150‑1 (MQ=255)
                                                                                                                                         aaaGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACtt                               >  2:1777255/1‑130 (MQ=255)
                                                                                                                                         aaaGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACtt                               <  1:1777255/130‑1 (MQ=255)
                                                                                                                                          aaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGa          <  2:289050/150‑1 (MQ=255)
                                                                                                                                          aaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGa          <  2:370854/150‑1 (MQ=255)
                                                                                                                                          aaGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGa          <  1:1232725/150‑1 (MQ=255)
                                                                                                                                           aGGGCCCTTCGCAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAg         >  1:2442304/1‑150 (MQ=255)
                                                                                                                                                ccTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTa                                        >  2:551846/1‑114 (MQ=255)
                                                                                                                                                ccTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTa                                        <  1:551846/114‑1 (MQ=255)
                                                                                                                                                  tGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCa                                           <  1:1363212/109‑1 (MQ=255)
                                                                                                                                                  tGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCa                                           >  2:1363212/1‑109 (MQ=255)
                                                                                                                                                  tGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCg  >  2:1915585/1‑150 (MQ=255)
                                                                                                                                                   gCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCTGAg                                                                                                        >  2:1901818/1‑47 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTGCTGCTTCTGTTATTTCCTTTTAAATCAATAAGTTAGCAAATATCTATCACTTCGGTGGTCGTGCATTTTGCAATCAGCGCGGTAAACCTGCATGCAAGATGCAAGTCCGCAAAAGGTTCCCGTCCCCATGAAAAAAGGGCCCTGCGAAGGGCCCTTTCTCTATTGCGTGACGTTCGCTTAGACGCCCGAGGCCTTGGCCGCAGGGACGTCTTTGATCGACAGCTTGATACGGCCGCGGTTGTCCACGTCCAGTACAAGCACTTCAACTTCCTGGCCTTCCTTGAGGATGTCG  >  NC_002947/5351878‑5352173

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