Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 182,035 A→T 10.3% intergenic (+24/+41) PP_5SB → / ← PP_0160 5S ribosomal RNA/ferrioxamine receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947182,0350AT10.3% 164.5 / 10.4 68intergenic (+24/+41)PP_5SB/PP_01605S ribosomal RNA/ferrioxamine receptor
Reads supporting (aligned to +/- strand):  ref base A (30/31);  new base T (4/3);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.65e-01

AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGC  >  NC_002947/181887‑182183
                                                                                                                                                       |                                                                                                                                                    
aGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGa                                                                                                                                                      >  2:293455/1‑149 (MQ=255)
  aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                                                                                                                    >  3:362874/1‑149 (MQ=255)
          gacgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgg                                                                                                                                            >  8:247462/1‑149 (MQ=255)
           acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAggg                                                                                                                                           <  5:98833/149‑1 (MQ=255)
           acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAggg                                                                                                                                           <  2:53366/149‑1 (MQ=255)
            cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                                                                                                                          <  1:273795/149‑1 (MQ=255)
             gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGAGATGTAGTTTGGGGTTTCCCCATGTGAGGTT‑‑AGGTAATAGTCAAGATTAAAACCCAAAACCCCTTTGAGGGATGCAAACAgggg                                                                                                                                          >  1:311120/1‑149 (MQ=14)
               ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCtt                                                                                                                                       <  4:362874/149‑1 (MQ=255)
                cATAGCGCATTGGAACCACCTGATCCCCTCCCGAACTCAGCAGTGAAACGCTGCATCGCCGATGGGAGTGGGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                      <  7:16972/149‑1 (MQ=18)
                cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                      <  4:350983/149‑1 (MQ=255)
                cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCt                                                                                                                                        >  1:215454/1‑147 (MQ=255)
                cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCt                                                                                                                                        <  2:215454/147‑1 (MQ=255)
                cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGCttt                                                                                                                                      <  7:6919/149‑1 (MQ=255)
                       ccttcgAACCACCTGAACCCAACCCCAACTCACCACTTAAACCAAGCATCGCCGATGGTAGGGGGGGGTTTCCCCATGTGAGAGTCCAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtt                                                                                                                                 <  1:111643/144‑1 (MQ=18)
                       cATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                      <  3:185666/142‑1 (MQ=255)
                       cATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                      >  4:185666/1‑142 (MQ=255)
                       cATTGGAACCACCTGATCCCATCCCGAACTAAGCAGTTAAACGATGCATCGCCGATGGTAGTGTGGGGTTTACCCATGTGAGAGT‑‑AGTTCATCGTCAAGATTCAAACCCAAAGCCCCTTTCTGCGATCCATAAAGGGGCTATGtaatt                                                                                                                                >  7:324140/1‑144 (MQ=17)
                         ttGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTcc                                                                                                                             >  6:18548/1‑149 (MQ=255)
                            gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCggg                                                                                                                          >  7:91480/1‑149 (MQ=255)
                            gAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCgg                                                                                                                           <  1:293455/148‑1 (MQ=255)
                              accaccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                                        <  3:261130/149‑1 (MQ=255)
                               ccaccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                                        >  8:193891/1‑148 (MQ=255)
                                 accTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCa                                                                                                                     <  5:267541/149‑1 (MQ=255)
                                  ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                                                                                                                          <  8:222435/127‑1 (MQ=255)
                                  ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                                                                                                                          >  7:222435/1‑127 (MQ=255)
                                         ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                             >  8:189354/1‑149 (MQ=255)
                                         ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                             <  4:115726/149‑1 (MQ=255)
                                           atcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa                                                                                                           <  5:39679/149‑1 (MQ=255)
                                           atcccGAACTCAGCAGTGAAACGATGCATCGCAGATGGTAGTGTGCGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAAAACAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTATTTTCCGGGTTTCAATATAGGt                                                                                                             >  3:108865/1‑147 (MQ=21)
                                                      agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGc                                                                                                >  8:100822/1‑149 (MQ=255)
                                                      agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCAGCGAGGCAGACac                                                                                                                                             >  6:170507/1‑103 (MQ=255)
                                                         agTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTGTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGa                                                                                             >  7:204095/1‑149 (MQ=255)
                                                             aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGt                                                                                         >  7:256496/1‑149 (MQ=255)
                                                               aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTgg                                                                                       >  8:140383/1‑149 (MQ=255)
                                                               aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTgg                                                                                       <  2:328974/149‑1 (MQ=255)
                                                                         gCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACAtt                                                                             <  7:100822/149‑1 (MQ=255)
                                                                               tGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGct                                                                        <  8:91480/149‑1 (MQ=255)
                                                                                 gTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGGCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGctc                                                                       <  7:294695/148‑1 (MQ=255)
                                                                                 gTAGTGTGGGGTTTCCACATGTGATAGT‑‑AGTTCATCGTCAAGATTAAAACCCAAATCCCCTGTCTTCGATGCAGACAGGTGCTTTGTCTTTTACGGGTTTCAATATCGGTTAGAACTGATAGCTGACCTTGGCGCTCACATTGCGctc                                                                       <  8:324140/148‑1 (MQ=255)
                                                                                    gtgtGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTc                                                                    <  8:204095/148‑1 (MQ=255)
                                                                                       tgGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAcc                                                                 >  3:47197/1‑148 (MQ=255)
                                                                                        ggggTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCa                                                               <  7:247462/149‑1 (MQ=255)
                                                                                         gggTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCAt                                                              <  5:18548/149‑1 (MQ=255)
                                                                                          ggTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCCGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                                        >  2:319569/1‑90 (MQ=255)
                                                                                          ggTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCCGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                                        <  1:319569/90‑1 (MQ=255)
                                                                                                cccATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                                        <  2:311120/148‑1 (MQ=255)
                                                                                                cccATGTGAGAGT‑‑AGGTCATAGTCAAGATTCAAACCCAAATCCCCTGTCTGGGATGCAGAAAGGGGCTTTGTGTTTTACGGGTTTCAATTTATGTTAGTACTGATAGATTAACGTGGCGCTCGCTTTGCGGTCTTCACGCATGTAGGag                                                       >  8:151717/1‑149 (MQ=255)
                                                                                                 ccATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg                                                          >  2:256488/1‑145 (MQ=255)
                                                                                                 ccATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg                                                          <  1:256488/145‑1 (MQ=255)
                                                                                                      tgAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                               <  1:326782/87‑1 (MQ=255)
                                                                                                      tgAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                               >  2:326782/1‑87 (MQ=255)
                                                                                                        agtgT‑‑AGGTCTTCGTCAAGATTCAAACCCAAAGGCCCTGTCTGGTATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                               <  6:214502/82‑1 (MQ=255)
                                                                                                        agagT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATCCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAAAATCg                                                                                                               >  5:214502/1‑85 (MQ=255)
                                                                                                                ggTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTgg                                          <  7:140383/148‑1 (MQ=255)
                                                                                                                ggTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGc                                         <  7:193891/149‑1 (MQ=255)
                                                                                                                  tCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGca                                        >  4:173808/1‑148 (MQ=255)
                                                                                                                        tCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAgg                                 >  1:195499/1‑149 (MQ=255)
                                                                                                                        tCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACACATGTAGCAGTAGTTCAGGCTGGCACATGAgg                                 >  8:123780/1‑149 (MQ=255)
                                                                                                                            gATTCAAACACAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAAAATCGGTTAGAACTGATAGCTGAACGTGGCGCTCACATTGCGATCTTCACCAATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGAt                             >  4:222153/1‑149 (MQ=255)
                                                                                                                             aTTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATg                            >  6:227158/1‑149 (MQ=255)
                                                                                                                                cAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTa                          <  3:173808/148‑1 (MQ=255)
                                                                                                                                    cccAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGGCTTGTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGCCCGTGGCGCTCACAGTGCGCTCTTCCCCCATGTAGCAGTAGTTCAGGCTGGGACATGAGGGGATGTAGCttt                     >  5:243899/1‑149 (MQ=255)
                                                                                                                                     ccAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCttt                     >  6:200513/1‑148 (MQ=255)
                                                                                                                                     acAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCttt                     >  7:293141/2‑148 (MQ=255)
                                                                                                                                        aaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt                 <  7:189354/149‑1 (MQ=255)
                                                                                                                                         aGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg                <  5:227158/149‑1 (MQ=255)
                                                                                                                                            cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGt              >  8:213893/1‑148 (MQ=255)
                                                                                                                                            cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGt              >  8:338301/1‑148 (MQ=255)
                                                                                                                                            cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGt              >  1:114244/1‑148 (MQ=255)
                                                                                                                                                       aTGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGc  <  2:195499/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGC  >  NC_002947/181887‑182183

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