Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 182,036 T→C 10.2% intergenic (+25/+40) 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,0360TC10.2% 171.9 / 10.1 69intergenic (+25/+40)PP_5SB/PP_01605S ribosomal RNA/ferrioxamine receptor
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base C (4/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 7.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.02e-01

AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGCG  >  NC_002947/181889‑182184
                                                                                                                                                      |                                                                                                                                                    
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                                                                                                                                       >  4:185666/1‑142 (MQ=255)
                     cATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                       <  3:185666/142‑1 (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                                                                                        <  2:328974/149‑1 (MQ=255)
                                                             aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTgg                                                                                        >  8:140383/1‑149 (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)
                                                                                                                                                      tGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGCg  <  8:256496/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGT‑‑AGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGGCG  >  NC_002947/181889‑182184

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: