Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 182,034 G→A 9.0% 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,0340GA 9.0% 211.0 / 8.9 79intergenic (+23/+42)PP_5SB/PP_01605S ribosomal RNA/ferrioxamine receptor
Reads supporting (aligned to +/- strand):  ref base G (38/32);  new base A (3/4);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

CTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGG  >  NC_002947/181894‑182182
                                                                                                                                            |                                                                                                                                                     
cTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGa                                                                                                                                                >  4:110621/1‑148 (MQ=255)
   gacgacGATAGAGCATTGGATCCAGCTGATCGCTTCCCGATCTGAGCAGTGATACGATGGATCGGCGATGGTAGTGTGGGGTTTCCGCTTGTGAGTGTAGGTCATCGTCTAGATTCTAACCCAAAGCCCCTGTCTGGGATGCAGACAg                                                                                                                                             <  2:176198/148‑1 (MQ=17)
    acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAggg                                                                                                                                           <  2:300364/149‑1 (MQ=255)
    acgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgg                                                                                                                                            <  1:84580/148‑1 (MQ=255)
     cgacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAgggg                                                                                                                                          <  3:110621/149‑1 (MQ=255)
      gacCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAGGGGc                                                                                                                                         <  6:308068/149‑1 (MQ=255)
        ccATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTAAAGATTCAAACCCAAAGCCCCTGTCTGTGATGCAGACAGGGGCtt                                                                                                                                       >  8:282394/1‑149 (MQ=255)
         cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCttt                                                                                                                                      <  4:290656/149‑1 (MQ=255)
         cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                                                                                                                            >  1:215052/1‑143 (MQ=255)
         cATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAgg                                                                                                                                            <  2:215052/143‑1 (MQ=255)
           tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCg                                                                                                                                                       <  6:322554/130‑1 (MQ=255)
           tAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCg                                                                                                                                                       >  5:322554/1‑130 (MQ=255)
            agagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCa                                                                                                                                                  >  2:281051/1‑134 (MQ=255)
            agagCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTc                                                                                                                                   >  4:287955/1‑149 (MQ=255)
              agCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCtt                                                                                                                                 >  4:256605/1‑149 (MQ=255)
                 aatGGAACCACCTCAACCCATCCCGAACTCAGCAGTGAAACGAAGCATCGCCCATGGTAGTGTGGGGTTTCCCCAAGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCAGACTGCGATCCAGACAGGGGCTATGTCTTTTc                                                                                                                              >  6:350359/3‑149 (MQ=18)
                 aTTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTc                                                                                                                              >  2:339334/1‑149 (MQ=255)
                 aTTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTc                                                                                                                              <  3:287955/149‑1 (MQ=255)
                    ggAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCg                                                                                                                            <  6:28609/148‑1 (MQ=255)
                        ccaccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGtt                                                                                                                        >  4:29554/1‑148 (MQ=255)
                          accTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGCGCttt                                                                                                                                      >  6:64763/1‑132 (MQ=255)
                           ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCaa                                                                                                                    >  2:142350/1‑149 (MQ=255)
                           ccTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCaa                                                                                                                    >  5:206089/1‑149 (MQ=255)
                              gATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAatat                                                                                                                 >  1:250391/1‑149 (MQ=255)
                                  ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCgg                                                                                                              <  4:283447/148‑1 (MQ=255)
                                  ccatcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                             >  7:323113/1‑149 (MQ=255)
                                   catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                             >  5:233133/1‑148 (MQ=255)
                                   catcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGt                                                                                                             >  6:64507/1‑148 (MQ=255)
                                    atcccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa                                                                                                           <  1:142350/149‑1 (MQ=255)
                                      cccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAg                                                                                                                                             >  2:309482/1‑113 (MQ=255)
                                      cccGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCATCGCAGACAg                                                                                                                                             <  1:309482/113‑1 (MQ=255)
                                           aCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa                                                                                                   <  5:64507/149‑1 (MQ=255)
                                           aCTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGa                                                                                                   <  3:256605/149‑1 (MQ=255)
                                               agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                                                                                                                    <  8:158481/97‑1 (MQ=255)
                                               agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATg                                                                                                                                                    >  7:158481/1‑97 (MQ=255)
                                               agcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGc                                                                                               >  6:310548/1‑149 (MQ=255)
                                                gcagTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCt                                                                                              <  6:233133/149‑1 (MQ=255)
                                                  agTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGa                                                                                            <  6:206089/149‑1 (MQ=255)
                                                   gTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTa                                                                                                           >  4:150093/1‑134 (MQ=255)
                                                      aaaCGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGt                                                                                        <  1:339334/149‑1 (MQ=255)
                                                         cGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGgc                                                                                     <  3:29554/149‑1 (MQ=255)
                                                          gATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTTGAACTGATAGCTGACCGTGgcg                                                                                    >  7:144565/1‑149 (MQ=255)
                                                                tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTc                                                                                                                      <  3:7090/110‑1 (MQ=255)
                                                                tCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTc                                                                                                                      >  4:7090/1‑110 (MQ=255)
                                                                 cGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                               >  2:245358/1‑116 (MQ=255)
                                                                 cGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCg                                                                                                               <  1:245358/116‑1 (MQ=255)
                                                                     gATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTgcg                                                                         <  8:296012/149‑1 (MQ=255)
                                                                               tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCAcc                                                                >  2:248510/1‑148 (MQ=255)
                                                                               tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTca                                                                                >  8:5147/1‑131 (MQ=255)
                                                                               tgGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTca                                                                                <  7:5147/131‑1 (MQ=255)
                                                                                 gggTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCa                                                              <  5:310548/148‑1 (MQ=255)
                                                                                   ggttCCCCATGTGAGAGTAGGTCATCGTAAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACATGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGTCACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATg                                                            <  1:248510/147‑1 (MQ=39)
                                                                                   gTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGt                                                           <  8:323113/149‑1 (MQ=255)
                                                                                    tttCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATg                                                            >  8:315489/1‑147 (MQ=255)
                                                                                       ccccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGCCAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCCTGGCGCTCACATTGCGCGATTCACCCATGtcg                                                         >  8:240735/1‑145 (MQ=255)
                                                                                        cccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                                       <  2:250391/148‑1 (MQ=255)
                                                                                         ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCag                                                      >  8:326762/1‑148 (MQ=255)
                                                                                         ccATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCag                                                      >  3:292665/1‑148 (MQ=255)
                                                                                             gtgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                                       >  4:245520/1‑143 (MQ=255)
                                                                                             gtgAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                                       <  3:245520/143‑1 (MQ=255)
                                                                                                agagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAg                                              >  1:189174/1‑149 (MQ=255)
                                                                                                 gagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAgg                                             >  7:29241/1‑149 (MQ=255)
                                                                                                 gagTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAgg                                             >  5:142335/1‑149 (MQ=255)
                                                                                                       gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCTATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGca                                       <  7:315489/149‑1 (MQ=255)
                                                                                                       gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCAt                                                             >  4:5264/1‑127 (MQ=255)
                                                                                                       gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCAt                                                             <  3:5264/127‑1 (MQ=255)
                                                                                                       gTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGca                                       >  8:301427/1‑149 (MQ=255)
                                                                                                             gTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAg                                 >  4:131282/1‑149 (MQ=255)
                                                                                                              tCAAGATTCAAACCCAAAGCCCCTGTCTGCCATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAgg                                >  1:32194/1‑149 (MQ=255)
                                                                                                               cAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGt                               >  7:334747/1‑149 (MQ=255)
                                                                                                                   aTTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATg                           >  4:252513/1‑149 (MQ=255)
                                                                                                                    ttCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGt                          <  8:29241/149‑1 (MQ=255)
                                                                                                                             aaaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                                       <  7:72175/111‑1 (MQ=255)
                                                                                                                             aaaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCa                                                       >  8:72175/1‑111 (MQ=255)
                                                                                                                             aaaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCAt                 >  4:324591/1‑149 (MQ=255)
                                                                                                                              aaGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGTTGATGTAGCTTTCAtt                <  8:334747/149‑1 (MQ=255)
                                                                                                                                gCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTgg              <  3:131282/149‑1 (MQ=255)
                                                                                                                                 ccccTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGt             <  7:301427/149‑1 (MQ=255)
                                                                                                                                  cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATATAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg               <  8:82444/146‑1 (MQ=255)
                                                                                                                                  cccTGTCTGCATCGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATATAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTg               >  7:82444/1‑146 (MQ=255)
                                                                                                                                   ccTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTaa           <  4:292665/149‑1 (MQ=255)
                                                                                                                                            gATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTGATTGGTAAGGTTGTTgg  >  4:328926/1‑149 (MQ=255)
                                                                                                                                            gATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTgg  >  6:324413/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                     
CTTGACGACCATAGAGCATTGGAACCACCTGATCCCATCCCGAACTCAGCAGTGAAACGATGCATCGCCGATGGTAGTGTGGGGTTTCCCCATGTGAGAGTAGGTCATCGTCAAGATTCAAACCCAAAGCCCCTGTCTGCGATGCAGACAGGGGCTTTGTCTTTTCCGGGTTTCAATATCGGTTAGAACTGATAGCTGACCGTGGCGCTCACATTGCGCTCTTCACCCATGTAGCAGTAGTTCAGGCTGGCACATGAGGTGATGTAGCTTTCATTGGTAAGGTTGTTGG  >  NC_002947/181894‑182182

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