Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,500,180 A→G 8.3% H39R (CAC→CGC) ‡ PP_2193 → ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,500,1800AG 8.3% 432.7 / 19.7 145H39R (CAC→CGC) ‡PP_2193ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base A (73/60);  new base G (5/7);  total (78/67)
Fisher's exact test for biased strand distribution p-value = 5.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 146 total aligned reads displayed.

GAGACTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGATGTCGCCG  >  NC_002947/2500032‑2500329
                                                                                                                                                    |                                                                                                                                                      
gagaCTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCAc                                                                                                                                                       >  1:2819730/1‑150 (MQ=255)
   aCTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCa                                                                                                                                                    >  1:4162128/1‑150 (MQ=255)
     tGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAgg                                                                                                                                                 <  1:2042730/150‑1 (MQ=255)
      gATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAggg                                                                                                                                                <  1:3023677/150‑1 (MQ=255)
      gATCAGCAACCCATTCTCATTTACGACGTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAggg                                                                                                                                                <  2:4428513/150‑1 (MQ=255)
      gATCAGCAACCCATGCACATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCAGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAggg                                                                                                                                                <  2:5042739/150‑1 (MQ=255)
       aTCAGCGACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCTGCGTGCAggggg                                                                                                                                              <  2:299638/150‑2 (MQ=255)
       aTCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAgggg                                                                                                                                               >  1:4027624/1‑150 (MQ=255)
       aTCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCTGCGTGCAggggg                                                                                                                                              <  2:3512816/150‑2 (MQ=255)
       aTCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCTGCGTGCAggggg                                                                                                                                              <  2:299460/150‑2 (MQ=255)
        tCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGc                                                                                                                                              >  2:4473186/1‑150 (MQ=255)
        tCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGc                                                                                                                                              <  1:5063304/150‑1 (MQ=255)
         cAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCTGCGTGCAgggggcc                                                                                                                                            <  2:4512098/150‑4 (MQ=255)
           gCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCcgc                                                                                                                                           >  1:3164754/1‑150 (MQ=255)
           gCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCcgc                                                                                                                                           >  1:3158635/1‑150 (MQ=255)
           gCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCTGCGTGCAgggggccga                                                                                                                                          <  1:1475770/150‑6 (MQ=255)
           gCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCTGCGTGCAgggggccga                                                                                                                                          <  1:1133/150‑6 (MQ=255)
            cAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCcgcg                                                                                                                                          >  1:2922445/1‑150 (MQ=255)
             accccATTCTCATTTACGACTGGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCGGTCGGGGTTCTTTTCACAGGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGAc                                                                                                                                        <  2:4300448/148‑1 (MQ=255)
             aaCCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGa                                                                                                                                         <  1:2596983/150‑1 (MQ=255)
              aCCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGAc                                                                                                                                        >  1:775027/1‑150 (MQ=255)
              aCCCATTCTCATTTACGACATGTCAGGAAGCTCATTCGTGCCTCTGATGCACAATCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGAc                                                                                                                                        >  1:3492873/1‑150 (MQ=255)
               cccATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACg                                                                                                                                       <  2:3165000/150‑1 (MQ=255)
               cccATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACg                                                                                                                                       >  2:2135062/1‑150 (MQ=255)
                ccATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACgg                                                                                                                                      >  1:2703958/1‑150 (MQ=255)
                  aTTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCt                                                                                                                                    <  1:2153614/150‑1 (MQ=255)
                   ttCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTg                                                                                                                                   >  1:541782/1‑150 (MQ=255)
                   ttCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTg                                                                                                                                   <  1:1013428/150‑1 (MQ=255)
                        aTTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGaa                                                                                                                              <  1:718097/150‑1 (MQ=255)
                          ttACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCGTGCAg                                                                                                                                                  <  1:775765/128‑1 (MQ=255)
                             cGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAc                                                                                                                         >  1:5311903/1‑150 (MQ=255)
                               aCTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAccc                                                                                                                       <  2:3625087/150‑1 (MQ=255)
                               aCTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAccc                                                                                                                       >  1:1455976/1‑150 (MQ=255)
                                 ttGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAg                                                                                                                     >  2:67543/1‑150 (MQ=255)
                                    tCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGACGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGAAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATg                                                                                                                  >  2:2583991/1‑150 (MQ=255)
                                      aGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACa                                                                                                                            <  2:2029231/138‑1 (MQ=255)
                                      aGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACa                                                                                                                            >  1:2029231/1‑138 (MQ=255)
                                      aGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAccc                                                                                                                       <  1:2403120/143‑1 (MQ=255)
                                      aGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGAccc                                                                                                                       >  2:2403120/1‑143 (MQ=255)
                                       ggAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAg                                                                                                               <  1:519937/150‑1 (MQ=255)
                                       ggAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAg                                                                                                               <  2:2819730/150‑1 (MQ=255)
                                        gAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGACCAGACCCAGATGCAgg                                                                                                              >  2:3021588/1‑150 (MQ=255)
                                         aaGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGc                                                                                                             <  1:4381701/150‑1 (MQ=255)
                                           gCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGccc                                                                                                           >  1:4815469/1‑150 (MQ=255)
                                             tCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGc                                                                                                         >  1:2487922/1‑150 (MQ=255)
                                               aTTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCa                                                                                                       >  1:1465592/1‑150 (MQ=255)
                                                ttCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAg                                                                                                      <  2:3327885/150‑1 (MQ=255)
                                                   gTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTc                                                                                                   >  2:4145046/1‑150 (MQ=255)
                                                   gTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAAGCCCGCCAGGTc                                                                                                   <  2:1627590/150‑1 (MQ=255)
                                                     gCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCa                                                                                                 <  2:3974035/150‑1 (MQ=255)
                                                     gCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCACGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCa                                                                                                 >  2:5250393/1‑150 (MQ=255)
                                                      ccTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAg                                                                                                <  2:2115225/150‑1 (MQ=255)
                                                          tGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCg                                                                                            >  2:2789575/1‑150 (MQ=255)
                                                            aTGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAc                                                                                          <  2:2487922/150‑1 (MQ=255)
                                                            aTGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAc                                                                                          >  2:252461/1‑150 (MQ=255)
                                                            aTGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGAc                                                                                          <  1:302020/150‑1 (MQ=255)
                                                                tCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGTTGCAg                                                                                                               <  1:4648094/125‑1 (MQ=255)
                                                                tCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAg                                                                                                               >  2:4648094/1‑125 (MQ=255)
                                                                tCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGc                                                                                      >  2:3594796/1‑150 (MQ=255)
                                                                  aaGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGTCa                                                                                    >  1:4343989/1‑150 (MQ=255)
                                                                  aaGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCa                                                                                    >  2:2044104/1‑150 (MQ=255)
                                                                                 cgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGAc                                                                     <  1:4473186/150‑1 (MQ=255)
                                                                                 cgTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGAc                                                                     <  1:4145598/150‑1 (MQ=255)
                                                                                      cgGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgc                                                                >  2:5327899/1‑150 (MQ=255)
                                                                                        gTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGgcgc                                                              <  1:252461/150‑1 (MQ=255)
                                                                                            ttttCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGAc                                                          <  1:4087374/150‑1 (MQ=255)
                                                                                             tttCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGc                                                                                                                                                     <  2:753710/59‑1 (MQ=255)
                                                                                             tttCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGc                                                                                                                                                     >  1:753710/1‑59 (MQ=255)
                                                                                               tCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCAc                                                       >  1:4008313/1‑150 (MQ=255)
                                                                                                   cATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTc                                                   <  2:3158635/150‑1 (MQ=255)
                                                                                                   cATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTc                                                   <  2:5227615/150‑1 (MQ=255)
                                                                                                   cATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCAAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTc                                                   <  2:2463786/150‑1 (MQ=255)
                                                                                                        cTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGa                                              >  2:4797946/1‑150 (MQ=255)
                                                                                                        cTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGa                                              >  2:1122738/1‑150 (MQ=255)
                                                                                                         ttGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGAc                                             >  2:658439/1‑150 (MQ=255)
                                                                                                              gTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGa                                                           >  1:4541375/1‑131 (MQ=255)
                                                                                                                 ttGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccgg                                     >  1:87654/1‑150 (MQ=255)
                                                                                                                 ttGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccgg                                     >  1:5256160/1‑150 (MQ=255)
                                                                                                                 ttGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccgg                                     >  2:4321820/1‑150 (MQ=255)
                                                                                                                 ttGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccgg                                     >  2:3895298/1‑150 (MQ=255)
                                                                                                                   gCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAggccggcc                                   >  1:3720973/1‑150 (MQ=255)
                                                                                                                    cAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCt                                  >  1:4691574/1‑150 (MQ=255)
                                                                                                                         ccGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCAcc                             >  1:116351/1‑150 (MQ=255)
                                                                                                                          cGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACct                            >  1:3661344/1‑150 (MQ=255)
                                                                                                                            tGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgc                          <  2:4651073/150‑1 (MQ=255)
                                                                                                                            tGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgc                          >  2:1226221/1‑150 (MQ=255)
                                                                                                                             ggCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGTCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct                         <  1:4398993/150‑1 (MQ=255)
                                                                                                                             ggCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACctgct                         <  2:3164754/150‑1 (MQ=255)
                                                                                                                                ccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTa                      >  1:1534772/1‑150 (MQ=255)
                                                                                                                                ccTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTa                      <  1:2135062/150‑1 (MQ=255)
                                                                                                                                      gCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTAcaccac                >  1:5263654/1‑150 (MQ=255)
                                                                                                                                        tCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAg              >  2:1818878/1‑150 (MQ=255)
                                                                                                                                        tCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAg              >  2:4667047/1‑150 (MQ=255)
                                                                                                                                            gccccTGCGTGCAGGGGGCCGCGACGGCTGCCGAACAGACCCAGATCCAGGCCCGCCAGGTCCAGCTCGAGCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACCGGCg           >  1:169943/2‑150 (MQ=255)
                                                                                                                                            cccccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGa          <  2:541782/150‑1 (MQ=255)
                                                                                                                                             ccccTTCACGCAGGGGCCGCGACGTCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCTTCCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCTCAGACCATGCCGGCCTGCACCTGCTGTACACCACAGGCGAt         <  1:3952581/150‑1 (MQ=255)
                                                                                                                                             ccccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGAt         <  2:775027/150‑1 (MQ=255)
                                                                                                                                             ccccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGAt         >  1:1637737/1‑150 (MQ=255)
                                                                                                                                               ccTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGATGt       <  2:3661344/150‑1 (MQ=255)
                                                                                                                                                cTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGATGTc      >  1:3886696/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
GAGACTGATCAGCAACCCATTCTCATTTACGACTTGTCAGGAAGCTCATTCGTGCCTCTGATGCTCAAGCCCTTGCACCGCCGTCGCGGTTCTTTTCACCATGCCTTGATGTCTTGCAGTGCCGTGGCCCTGTTGCTCGGCCCCCTGCACGCAGGGGCCGCGACGGCTGCCGAACAGACCCAGATGCAGGCCCGCCAGGTCCAGCTCGACCTGCCAGCCCAACCGCTGGACCAGGCGCTGACCACTTTCGCCGACCAGGCCGGCCTGCACCTGCTGTACACCACAGGCGATGTCGCCG  >  NC_002947/2500032‑2500329

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