Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 182,035 A→T 11.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,0350AT11.3% 198.5 / 13.5 78intergenic (+24/+41)PP_5SB/PP_01605S ribosomal RNA/ferrioxamine receptor
Reads supporting (aligned to +/- strand):  ref base A (35/33);  new base T (5/4);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTG  >  NC_002947/181889‑182181
                                                                                                                                                  |                                                                                                                                                   
aaTTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                                                                                                                   >  2:138029/1‑149 (MQ=255)
         acgacCATAGAGCCTTTGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGGGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCATGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgg                                                                                                                                           <  6:195740/148‑1 (MQ=255)
         acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                                                                                                                           <  6:283491/148‑1 (MQ=255)
          agacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                          <  4:69149/147‑1 (MQ=255)
           gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                          >  5:186195/1‑147 (MQ=255)
           gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                          <  6:186195/147‑1 (MQ=255)
           gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGc                                                                                                                                        >  8:142272/1‑149 (MQ=255)
              cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                          >  3:294450/1‑144 (MQ=255)
              cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                          <  4:294450/144‑1 (MQ=255)
                tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGt                                                                                                                                   <  6:211837/149‑1 (MQ=255)
                  gagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgggt                                                                                                                                 <  2:126140/149‑4 (MQ=255)
                  gagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGCTTTgggt                                                                                                                                 <  7:142272/149‑4 (MQ=255)
                     cATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtttt                                                                                                                              <  4:231383/149‑1 (MQ=255)
                      aTTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTc                                                                                                                             >  1:207665/1‑149 (MQ=255)
                           aaCCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGt                                                                                                                        >  6:233093/1‑149 (MQ=255)
                                ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCaa                                                                                                                   <  2:38365/149‑1 (MQ=255)
                                   gATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAata                                                                                                                 >  3:275086/1‑148 (MQ=255)
                                       ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAatat                                                                                                                <  4:252498/145‑1 (MQ=255)
                                       ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAatat                                                                                                                >  3:252498/1‑145 (MQ=255)
                                            ccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGaa                                                                                                       <  4:301306/149‑1 (MQ=255)
                                             cGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTg                                                                                                                                    >  3:130047/1‑119 (MQ=255)
                                             cGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTg                                                                                                                                    <  4:130047/119‑1 (MQ=255)
                                             cGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAAc                                                                                                      >  5:15283/1‑149 (MQ=255)
                                             cGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCACCTGTCTTCGATGCAGACAGGGGCTTTGTCTTTTCAGGGTTTCAAAAAAGGATAGAAc                                                                                                      >  3:13051/1‑149 (MQ=21)
                                                 cTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGAt                                                                                                  <  2:207665/149‑1 (MQ=255)
                                                  tCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATa                                                                                                 >  6:204715/1‑149 (MQ=255)
                                                      cagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt                                                                                              >  6:331513/1‑148 (MQ=255)
                                                       agTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGa                                                                                            <  1:138029/149‑1 (MQ=255)
                                                       agTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGa                                                                                            >  7:61326/1‑149 (MQ=255)
                                                        gTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGAc                                                                                           >  6:176884/1‑149 (MQ=255)
                                                        gTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGAc                                                                                           >  1:158474/1‑149 (MQ=255)
                                                         tGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGAcc                                                                                          <  5:331513/149‑1 (MQ=255)
                                                            aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTg                                                                                       >  3:272065/1‑149 (MQ=255)
                                                            aaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTg                                                                                       <  4:272065/149‑1 (MQ=255)
                                                             aCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTgg                                                                                      <  2:167824/149‑1 (MQ=255)
                                                                 tGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCt                                                                                  <  5:173653/149‑1 (MQ=255)
                                                                  gCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCt                                                                                  >  8:102109/1‑148 (MQ=255)
                                                                    aTCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTca                                                                                >  2:63003/1‑148 (MQ=255)
                                                                     tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTcaca                                                                              <  6:15908/149‑1 (MQ=255)
                                                                       gCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACAtt                                                                            <  5:182593/149‑1 (MQ=255)
                                                                         cGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgc                                                                          <  5:216103/149‑1 (MQ=255)
                                                                           aTGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgcgc                                                                        <  4:275086/149‑1 (MQ=255)
                                                                               tAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGctct                                                                     >  7:118055/1‑148 (MQ=255)
                                                                                 gtgtGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCa                                                                  <  5:176884/149‑1 (MQ=255)
                                                                                       ggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACAtt                                                                            >  2:169814/1‑133 (MQ=255)
                                                                                       ggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACAtt                                                                            <  1:169814/133‑1 (MQ=255)
                                                                                            ccccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGCTGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAg                                                        >  6:45263/1‑147 (MQ=255)
                                                                                             cccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCag                                                     >  6:283834/1‑149 (MQ=255)
                                                                                             cccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                                      <  5:204715/148‑1 (MQ=255)
                                                                                               cATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCTCATGtgccagc                                                    >  6:174688/1‑143 (MQ=255)
                                                                                               cATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCag                                                     >  4:78973/1‑147 (MQ=255)
                                                                                                   tgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGtt                                                >  1:65482/1‑148 (MQ=255)
                                                                                                     agagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCa                                              <  6:15283/148‑1 (MQ=255)
                                                                                                      gagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAgg                                            <  5:283834/149‑1 (MQ=255)
                                                                                                         tAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTg                                         >  3:94661/1‑149 (MQ=255)
                                                                                                          aGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAAATGATAGCTGACCGTGGCGCTCACATTGCGATCTTAACCCATGTAGAAGTAGTTCatgatg                                         >  7:221942/1‑143 (MQ=255)
                                                                                                            gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAACCATGTAGCAGTAGTTCAGGCTGGca                                      >  2:248775/1‑149 (MQ=255)
                                                                                                            gTCATCGTCAAGATTCAAACCCAAAGCCACTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCCATATCGGGTAGAACTGAAAGCTGACCATGGCGATCACATTGCGCTCTTCAAGCGTGTAGCAGTAGTTCAGGCTGGc                                       >  7:120328/1‑148 (MQ=255)
                                                                                                             tCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGca                                      >  2:169538/1‑148 (MQ=255)
                                                                                                                      aGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGa                            <  1:63003/149‑1 (MQ=255)
                                                                                                                        aTTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTCCAGGCTGGcac                                     <  3:197394/138‑1 (MQ=255)
                                                                                                                        aTTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTCCAGGCTGGcac                                     >  4:197394/1‑138 (MQ=255)
                                                                                                                          tCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATg                          >  3:292232/1‑147 (MQ=255)
                                                                                                                           cAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTa                        <  4:316425/148‑1 (MQ=255)
                                                                                                                            aaaCCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGc                      >  3:204357/1‑149 (MQ=255)
                                                                                                                                 cAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTg                                                                                       >  3:107184/1‑80 (MQ=255)
                                                                                                                                 cAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTg                                                                                       <  4:107184/80‑1 (MQ=255)
                                                                                                                                    aGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg              <  8:61326/149‑1 (MQ=255)
                                                                                                                                    aGCCCCTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg              >  4:149697/1‑149 (MQ=255)
                                                                                                                                       cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt               <  4:256492/145‑1 (MQ=255)
                                                                                                                                       cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAtt               >  3:256492/1‑145 (MQ=255)
                                                                                                                                       cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGt            >  4:118954/1‑148 (MQ=255)
                                                                                                                                        ccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGt            >  1:13328/1‑147 (MQ=255)
                                                                                                                                               gCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGGAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGgttgtt   >  4:2030/1‑149 (MQ=255)
                                                                                                                                                cGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGttgttg  <  3:78973/149‑1 (MQ=255)
                                                                                                                                                cGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGttgttg  <  2:158474/149‑1 (MQ=255)
                                                                                                                                                cGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGttgttg  <  1:169538/149‑1 (MQ=255)
                                                                                                                                                cGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGGTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGttgttg  >  7:125685/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AATTTCTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTG  >  NC_002947/181889‑182181

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