Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 182,034 G→A 10.2% intergenic (+23/+42) 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,0340GA10.2% 176.6 / 8.7 69intergenic (+23/+42)PP_5SB/PP_01605S ribosomal RNA/ferrioxamine receptor
Reads supporting (aligned to +/- strand):  ref base G (31/30);  new base A (4/3);  total (35/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.87e-01

AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGT  >  NC_002947/181887‑182171
                                                                                                                                                      |                                                                                                                                         
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                                                                                                                              >  7:222435/1‑127 (MQ=255)
                                  ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                                                                                                              <  8:222435/127‑1 (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                                                                                                            <  1:319569/90‑1 (MQ=255)
                                                                                          ggTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCCGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                            >  2:319569/1‑90 (MQ=255)
                                                                                                cccATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                            <  2:311120/148‑1 (MQ=255)
                                                                                                cccATGTGAGAGT‑‑AGGTCATAGTCAAGATTCAAACCCAAATCCCCTGTCTGGGATGCAGAAAGGGGCTTTGTGTTTTACGGGTTTCAATTTATGTTAGTACTGATAGATTAACGTGGCGCTCGCTTTGCGGTCTTCACGCATGTAGGag                                           >  8:151717/1‑149 (MQ=255)
                                                                                                 ccATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg                                              <  1:256488/145‑1 (MQ=255)
                                                                                                 ccATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg                                              >  2:256488/1‑145 (MQ=255)
                                                                                                      tgAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                   >  2:326782/1‑87 (MQ=255)
                                                                                                      tgAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                   <  1:326782/87‑1 (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)
                                                                                                                                                      |                                                                                                                                         
AGAATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGT  >  NC_002947/181887‑182171

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: