Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,782 G→A 22.0% H248H (CAC→CAT) ‡ ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,7820GA22.0% 422.1 / 117.3 214H248H (CAC→CAT) ‡ftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (86/81);  new base A (22/25);  total (108/106)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

CGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGG  >  NC_002947/5367634‑5367930
                                                                                                                                                    |                                                                                                                                                    
cgatgaTGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGt                                                                                                                                                     <  2:1810279/150‑1 (MQ=255)
 gatgatGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTg                                                                                                                                                    >  1:2993755/1‑150 (MQ=255)
 gatgatGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTg                                                                                                                                                    >  1:4982645/1‑150 (MQ=255)
    gatgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGctt                                                                                                                                                 >  1:2742235/1‑150 (MQ=255)
      tgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttct                                                                                                                                               >  1:5036132/1‑150 (MQ=255)
      tgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttct                                                                                                                                               <  2:3994282/150‑1 (MQ=255)
      tgCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttct                                                                                                                                               <  1:2330399/150‑1 (MQ=255)
       gCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttctt                                                                                                                                              <  2:2741428/150‑1 (MQ=255)
       gCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttctt                                                                                                                                              <  2:3730254/150‑1 (MQ=255)
        ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                                                                                             >  2:2569921/1‑150 (MQ=255)
         catcatTCATCTCGAAGCAGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACTCTCGTCGTGACCTCCGCCCATGCCTGCGCCACTATGGCTACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGTGCGTGCTTCTTgg                                                                                                                                            <  2:1584039/150‑1 (MQ=255)
          atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGc                                                                                                                                           >  2:4393014/1‑150 (MQ=255)
          atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTGCTTGGc                                                                                                                                           >  2:4946581/1‑150 (MQ=255)
            cattcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCt                                                                                                                                         <  1:4053126/150‑1 (MQ=255)
             attcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTg                                                                                                                                        >  1:3672688/1‑150 (MQ=255)
              ttcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttcttcttcatc                                                                                                                                       <  1:4487547/150‑8 (MQ=255)
              ttcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGc                                                                                                                                       <  2:4792486/150‑1 (MQ=255)
               tcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCt                                                                                                                                      >  1:3718315/1‑150 (MQ=255)
                  tctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGa                                                                                                                                   >  2:4522051/1‑150 (MQ=255)
                    tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                                                                                                                 >  2:3390000/1‑150 (MQ=255)
                    tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                                                                                                                 >  2:3391952/1‑150 (MQ=255)
                    tcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAAc                                                                                                                                 >  1:3915792/1‑150 (MQ=255)
                         gCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTc                                                                                                                            <  1:2097925/150‑1 (MQ=255)
                           cGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgc                                                                                                                          <  1:1884641/150‑1 (MQ=255)
                            gTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTcgcg                                                                                                                         >  1:1838527/1‑150 (MQ=255)
                            atCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacg                                                                                                                         >  2:101391/2‑150 (MQ=255)
                             tCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacgg                                                                                                                        >  2:1757060/1‑149 (MQ=255)
                              ccATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGca                                                                                                                       >  2:3848818/1‑150 (MQ=255)
                                 tctcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCTCGATGGCGACCAACGGCGTCGTTTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCacac                                                                                                                    <  1:3802804/150‑1 (MQ=255)
                                  ctcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACg                                                                                                                   >  1:2623926/1‑150 (MQ=255)
                                   tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACgg                                                                                                                  <  1:3605761/150‑1 (MQ=255)
                                    cgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGc                                                                                                                 <  2:531864/148‑1 (MQ=255)
                                      accaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTg                                                                                                               <  1:1454551/150‑1 (MQ=255)
                                      accaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCtt                                                                                                               <  1:1107363/150‑2 (MQ=255)
                                      accaGAAGCTGGTACAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACAACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTACTAGGCCTGTTCGAACATGTCACGGACGCGGCta                                                                                                               >  1:4186854/1‑149 (MQ=255)
                                       ccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTgg                                                                                                              <  1:3273860/150‑1 (MQ=255)
                                        caGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGc                                                                                                             >  1:1986071/1‑150 (MQ=255)
                                        caGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGAAGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGc                                                                                                             >  1:4561798/1‑150 (MQ=255)
                                         agaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCa                                                                                                            >  2:1064883/6‑150 (MQ=255)
                                           cAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCAcc                                                                                                          >  2:1447651/1‑150 (MQ=255)
                                             aCTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGa                                                                                                        <  1:518542/150‑1 (MQ=255)
                                                      gggTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGa                                                                                               >  2:1137045/1‑150 (MQ=255)
                                                        gTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAAc                                                                                             <  2:1937235/150‑1 (MQ=255)
                                                             ttcgcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTc                                                                                        <  1:2891887/143‑1 (MQ=255)
                                                              tcgcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTca                                                                                       <  2:3275220/144‑2 (MQ=255)
                                                              tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTcc                                                                                       >  2:2391080/1‑150 (MQ=255)
                                                              tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTcc                                                                                       >  1:2953525/1‑150 (MQ=255)
                                                              tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTcc                                                                                       >  1:2866656/1‑150 (MQ=255)
                                                                aCGCTCGTCGTGATCGCCGCCCATGCCGGCGCCTCGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCAc                                                                                     <  1:1054362/150‑1 (MQ=255)
                                                                  gCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGa                                                                                   >  1:4767669/1‑150 (MQ=255)
                                                                   cTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGaa                                                                                  <  2:2942993/150‑1 (MQ=255)
                                                                    tcgtcgTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAg                                                                                 >  1:903924/1‑150 (MQ=255)
                                                                      gtcgtGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTc                                                                               >  2:4352706/1‑150 (MQ=255)
                                                                        cgttgccgccgCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTcag                                                                             >  1:4446497/6‑148 (MQ=255)
                                                                           gACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAc                                                                          <  1:3652638/150‑1 (MQ=255)
                                                                            aCCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAcc                                                                         >  1:1465568/1‑150 (MQ=255)
                                                                              cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAAc                                                                          >  1:820667/1‑147 (MQ=255)
                                                                              cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                       >  2:4950577/1‑150 (MQ=255)
                                                                              cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                       >  1:633411/1‑150 (MQ=255)
                                                                                   cccATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATg                                                                  <  1:1282158/150‑1 (MQ=255)
                                                                                    ccATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATgg                                                                 <  2:4324980/150‑1 (MQ=255)
                                                                                     aagaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGt                                                                <  2:793100/146‑1 (MQ=255)
                                                                                     aagaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGt                                                                >  1:3976061/5‑150 (MQ=255)
                                                                                      aTGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTg                                                               >  1:1022562/1‑150 (MQ=255)
                                                                                       tGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTga                                                              <  1:3391952/150‑1 (MQ=255)
                                                                                       gaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTga                                                              <  2:3314582/148‑1 (MQ=255)
                                                                                        gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                             <  1:2531581/150‑1 (MQ=255)
                                                                                        gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                             <  2:1730064/150‑1 (MQ=255)
                                                                                        accGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTgaa                                                             <  2:4182869/149‑1 (MQ=255)
                                                                                        accGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTgaa                                                             <  2:1053574/149‑1 (MQ=255)
                                                                                          cGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTgaaga                                                           >  2:3644994/1‑150 (MQ=255)
                                                                                            gcaccACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAAc                                                         >  1:2542015/4‑150 (MQ=255)
                                                                                            gcaccACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAAc                                                         <  1:1429764/147‑1 (MQ=255)
                                                                                             cgcCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACg                                                        <  2:1465568/150‑1 (MQ=255)
                                                                                             cgcCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACg                                                        <  2:2448675/150‑1 (MQ=255)
                                                                                                cACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTa                                                     >  2:2022813/1‑150 (MQ=255)
                                                                                                 aCGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTAc                                                    >  2:4478298/1‑150 (MQ=255)
                                                                                                 aCGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTAc                                                    >  2:5014065/1‑150 (MQ=255)
                                                                                                    aTGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACttt                                                 >  1:1791622/1‑150 (MQ=255)
                                                                                                      ggCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACattgg                                               <  2:651794/150‑3 (MQ=255)
                                                                                                      ggCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACattgg                                               <  2:171763/150‑3 (MQ=255)
                                                                                                       gCGACCGACGGCGTCGATTTCGTCGATGAAGATGACGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGTACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACAttggc                                              <  1:4787628/150‑1 (MQ=255)
                                                                                                            cAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgcttcgc                                         <  1:2155473/150‑1 (MQ=255)
                                                                                                             aaCGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGcc                                        >  1:23898/1‑150 (MQ=255)
                                                                                                              aCGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGccc                                       >  2:1204282/1‑150 (MQ=255)
                                                                                                               cGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCcg                                      >  1:1661353/1‑150 (MQ=255)
                                                                                                                   gTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGCCagcaac                                  <  2:2963859/150‑4 (MQ=255)
                                                                                                                   gTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGCCagcaac                                  >  1:452091/1‑147 (MQ=255)
                                                                                                                        tttCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCt                             <  1:3618804/150‑1 (MQ=255)
                                                                                                                            gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                         >  1:4129704/1‑150 (MQ=255)
                                                                                                                            gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                         <  2:2114336/150‑1 (MQ=255)
                                                                                                                            gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGACCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                         >  2:3502237/1‑150 (MQ=255)
                                                                                                                              cgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCca                       >  1:1376413/1‑150 (MQ=255)
                                                                                                                               gATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCcag                      <  2:1577126/150‑1 (MQ=255)
                                                                                                                                  taaGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGTACATGTCGAGCACACTGCTTGAAGCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAATAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCcaacag                   <  2:4582105/149‑1 (MQ=255)
                                                                                                                                     gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt                >  1:3512140/1‑150 (MQ=255)
                                                                                                                                     gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt                >  1:4715463/1‑150 (MQ=255)
                                                                                                                                          tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac           <  2:3006371/150‑1 (MQ=255)
                                                                                                                                             aCGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccgg        >  1:2398501/1‑150 (MQ=255)
                                                                                                                                                   cGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccggtaccgg  >  2:3323237/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACCGG  >  NC_002947/5367634‑5367930

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: